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SECTION  I 
INTRODUCTION 


The  present  report  gives  a  version  of  linearized  airfoil 
theory  for  subsonic  oscillatory  flows  which  avoids  the  frequently- 
used  idea  of  concentrating  the  pressures  into  lines  or  points.  As 
usual,  the  linearized  flow  differential  equation  is  solved  by 
means  of  fundamental  solutions.  The  boundary  conditions  at  the 
planform  then  are  expressed  by  an  integral  equation.  In  this 
report,  the  acceleration  potential  is  used.  It  has  the  advantage 
of  giving  directly,  without  differentiations,  the  pressure 
distributions.  Consequently,  the  Kutta  condition  expresses  itself 
very  simply  by  the  postulate  that  there  is  no  pressure  difference 
between  the  upper  and  lower  side  of  the  wing.  However,  the 
singularities  which  occur  in  the  governing  integral  equation  are 
very  strong  and  this  causes  conceptual  and  practical 
difficulties. 

For  numerical  purposes,  one  always  expresses  the  pressure 
distribution  in  terms  of  a  finite  number  of  parameters.  In  the 
vortex  lattice  method  (which  can  be  viewed  as  a  discretization  of 
the  integral  equation  formulation)  the  pressures  are  concentrated 
into  lines  along  which  the  force  per  unit  of  length  is  piece-wise 
constant.  For  such  lines  (and  the  pertinent  trailing  vortices), 
the  upwash  can  be  computed  by  the  Blot-Savart  law.  Dowell  and 
Ueda  (Refs.  1,2),  on  the  other  hand,  concentrate  the  pressures 
into  "pressure  points;"  although  for  points  of  the  wing  surface 
lying  in  the  wake  of  the  pressure  points,  this  concept  must  be 
modified.  Dowell  and  Ueda  do  this  without  further  explanation  by 
reference  to  Mangler’s  work  (Ref.  3).  (It  seems  to  me  that 
Mangler's  approach  is  applicable  only  if  for  the  upwash  at  points 
of  the  wake  one  replaces  the  point  force  by  a  line  force,  or 
alternatively  if  one  averages  the  upwash  along  a  line  in  the 
wake.)  To  express  the  boundary  conditions  at  the  plan  form,  one 
usually  matches  at  a  sufficient  number  of  "control"  points,  the 
upwash  expressed  in  terms  of  the  pressure  parameters  and  the 
upwash  given  by  che  boundary  conditions.  One  then  obtains  a 
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linear  system  (with  a  full  matrix)  from  which  the  pressure 
parameters  are  determined. 

The  upwash  field  obtained  from  the  pressure  field  is  by  no 
means  smooth.  The  choice  of  control  points  therefore  introduces 
an  element  of  arbitrariness.  Nevertheless,  these  methods  are 
successful,  although  perhaps  not  overly  accurate. 

In  the  present  study,  the  author  reduces  the  arbitrariness 
caused  by  the  concept  of  pressure  points  or  pressure  lines  and  by 
the  choice  of  control  points.  (Some  arbitrariness  is  inherent  in 
any  discretization.)  Within  surface  elements  (preferably 
triangles),  an  expression  closely  related  to  the  pressure  is 
approximated  by  linear  functions.  The  pressure  is  continuous  as 
one  proceeds  from  one  element  to  its  neighbors.  In  a  cruder  form, 
one  may  also  use  constant  pressure  elements.  (Even  in  the  latter 
case,  the  wake  has  only  logarithmic  singularities  if  no  side  of 
the  (polygonal)  element  is  parallel  to  the  wake  streamlines.)  The 
arbitrariness  due  to  choice  of  the  control  points  is  avoided  by 
matching  the  upwash  velocities  in  the  average  over  control  areas 
(which  include  lines  along  which  the  flow  field  is  singular). 

While  this  is  conceptually  satisfactory,  it  complicates  the 
procedure  since  it  requires  further  integrations.  For  control 
areas  at  a  distance  from  the  pressure  areas  and  their  wakes,  the 
functions  to  be  integrated  are  smooth.  The  integration  can  then 
be  carried  out  numerically,  for  instance,  by  (low  order)  Gaussian 
integration.  The  results  are  rather  close  to  those  obtained  with 
the  idea  of  pressure  and  control  points.  (This  is  the  reason  for 
the  success  of  the  method  of  Ueda  and  Dowell.)  In  the  vicinity  of 
the  pressure  element,  or  of  its  wake,  it  is  preferable  to 
determine  the  part  of  the  upwash  in  which  singularities  occur 
analytically.  One  then  uses  a  mixed  numerical  analytical  method. 
This  works  well  for  control  areas  close  to  the  wake  of  the 
pressure  element  but  at  some  distance  from  the  pressure  area. 

If  the  control  element  is  close  to  the  pressure  element,  a 
corresponding  splitting  of  the  upwash  is  possible  in  principle. 

But  then  even  the  "smooth"  function  to  be  integrated  proves  to  be 


rather  intractable  for  numerical  integration  methods. 

Fortunately,  all  necessary  integrations  can  be  carried  out 
analytically  in  terms  of  elementary  transcedental  functions.  This 
is  possible  because  of  several  factors:  (1)  Ueda  has  provided  a 
development  of  the  kernel  of  the  integral  equation  with  respect  to 
the  reduced  frequency.  If  the  elements  are  close  to  each  other, 
only  the  terms  of  the  lowest  order,  which  are  relatively  simple, 
are  of  importance.  (2)  If  the  wing  surface  is  subdivided  into 
triangles,  then  it  is  possible  to  represent  the  pressure 
distribution  (roughly  speaking)  by  linear  functions  which  are 
continuous,  and  at  the  same  time,  allow  one  to  carry  out  the 
necessary  Integrations.  (3)  The  choice  of  the  coordinate  systems 
is  important.  To  obtain  the  upwash  at  a  point  (x,y)  one  uses 
coordinates  which  have  this  point  as  origin.  One  of  the  two 
integrations  over  the  pressure  element  is  then  trivial.  For  the 
integration  over  the  control  elements,  one  uses  as  origins  the 
corner  points  of  the  pressure  elements  with  a  similar  effect. 

The  basic  equations  and  the  formulae  of  Ueda  (with  some 
extension)  are  derived  in  Appendix  A  and  B.  The  integral  equation 
directly  written  for  the  plane  of  the  wing  is  meaningless.  The 
author  found  it  preferable  to  go  back  to  the  original  meaning  and 
include  the  limiting  process  in  which  one  approaches  the  plane  of 
the  wing  from  above  or  below.  Mostly,  but  not  always,  one  obtains 
the  same  result  as  in  the  usual  less  cautious  approach.  Integra¬ 
tion  formulae  needed  for  the  determination  of  the  upwash  and  later 
for  the  averaging  over  the  control  area  are  first  derived  for 
distant  points  close  to  the  wake  and  later  for  cases  where 
pressure  elements  and  control  elements  are  close  to  each  other. 

The  complexity  of  the  resulting  formulae  has  given  the 
author  some  uneasiness  because  of  the  danger  of  errors.  The 
analytical  procedure  is  shown  in  sufficient  detail  so  that  the 
reader  will  be  able  to  check  the  results  himself. 

A  survey  of  essential  ideas  and  a  listing  of  the  formulae 
needed  in  the  computation  is  found  in  Section  VIII. 
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SECTION  II 
BASIC  EQUATIONS 


Let  x,y,z  be  a  system  of  Cartesian  coordinates.  We  consider 
the  subsonic  oscillatory  flow  over  a  wing  in  a  linearized 
approach.  The  plane  of  the  wing  is  given  by  z  =  0 .  Coordinates 
within  the  wing  planform  corresponding  to  x  and  y  are  5  and  n 
respectively.  The  free  stream  velocity  is  U,  the  free  stream  Mach 
number  M,  and  the  free  stream  density  p.  All  lengths  are  made 
dimensionless  with  a  characteristic  length  L.  The  time  dependence 
of  the  oscillatory  motion  is  given  by  the  factor  exp(ivt),  where  t 
is  the  time.  The  dimensionless  frequency  is  k  =  vL/U.  The 
deviation  of  the  pressure  from  the  free  stream  pressure  is  made 
dimensionless  with  p  U“/2  and  denoted  by  p.  The  dimensionless 
pressure  difference  between  the  lower  and  upper  sides  of  the  wing 
is  denoted  by  Ap.  The  dimensionless  upwash  at  a  point  x,y,z  is 
denoted  by  w.  The  upwash  at  a  point  x,y,z  expressed  in  terms  of 
the  pressure  difference  is  given  by 

w(x,y,z)  -  (8tr)“^  ~  [z  jj  Ap(5.n)K(^-x,  r,i<)d^dri]  (1) 

A 

The  region  of  integration  is  the  wing  area  A.  K  is  given  by  the 
classical  formula 


K(5-x,r,k)  -  exp(ik(C-x)) 


Mexp( ikV ) 


R(V^+r^)’ 


exp( ikv) 
(V^?2)3/2 


dv 


(2) 


with 

r^  =  (n-y)^  z^ 

•  [(x-U^  +  3^?^] 

V  -  (1 -M^)"^(-(^-x)-MR) 


i4 


(3) 


(3)  cont'd. 


A  derivation  is  found  in  Appendix  A.  We  introduce  the 
Prandtl  Glauert  transformation 


with 


Then 


-  -  1  -1 

w(x,6  y,6  z)  -  w(x,y,z) 


Ap(5.8"'n)  -  Ap(t,n) 


6'^K(^-x,b"^ r,k)  -  K(^-x,r,k) 


(n-y)^  + 


w(x,y,z)  -  (8ir)~^  Tz”^^  jj  . n)K( ^-x , r , k) d^dn)  (5) 

A 


The  power  series  development  with  respect  to  the  reduced  frequency 
k  of  K,  Eq.  (2),  has  been  derived  by  Ueda.  (The  result  is  by  no 
means  obvious,  because  the  integrals  which  arise  if  one  develops 
the  integrand  of  Eq.  (2)  with  respect  to  k  do  not  converge.)  In 
Appendix  B  these  formulae  have  been  rederived.  Simplifications 
which  are  rather  important  for  the  present  approach  arise  if  one 
partially  combines  the  two  terms  on  the  right  of  Eq.  (2).  The 
form  used  here  is  given  in  Eqs.  (B.63,  B.64,  and  B.65). 


The  pressure  difference  Ap  is  found  from  an  integral 
equation,  obtained  from  Eq.  (5)  by  making  the  limiting  process  z*0 
on  the  right  and  by  substituting  on  the  left  the  upwash 
which,  of  course,  is  known  over  the  wing.  Because  of  the 
singularity  of  K  the  evaluation  of  the  right-hand  side  involves  a 
limiting  process  which  is  not  trivial  except  for  points  (x,y)  at 
some  distance  from  the  point  (5,n)  and  its  wake.  It  is  true  one 

obtains  ultimately  lim  -r— (z  ff..)  -  lim  ff  ...,  except  for  one 

z*0  z-^0 

important  exception  which  justifies  the  cautious  approach  taken  in 
this  report.  A  reference  to  the  work  of  Mangier  (Ref.  4)  does  not 
seem  to  be  sufficient. 

In  the  approach  taken  here  the  wing  planform  is  divided 
into  elements  and  within  each  element  the  (unknown)  pressure  is 
approximated  by  a  linear  combination  of  shape  functions.  The 
coefficients  by  which  the  shape  functions  are  multiplied  are  the 
pressure  parameters  so  far  unknown.  The  upwash  distribution  is 
expressed  in  terms  of  these  parameters.  Usually  a  system  of 
equations  for  these  parameters  is  obtained  by  equating  the  upwash 
at  certain  points,  called  upwash  points  or  control  points,  with 
the  upwash  given  by  the  boundary  conditions.  The  upwash 
pertaining  to  the  chosen  pressure  distribution  is  by  no  means 
smooth.  At  the  element  boundaries  and  certain  lines  of  the  wake 
pertaining  to  a  pressure  element  it  goes  logarithmically  to 
infinity;  for  constant  pressure  elements  it  even  behaves  as  d  ^ 
where  d  is  the  distance  from  the  boundary.  One  chooses  the 
control  points  at  a  distance  from  these  singularities,  but  in  any 
case,  the  results  will  be  rather  Inaccurate  (except  of  course  if 
the  points  (x,y)  are  at  some  distance  from  the  point  (e.n)  and  its 
wake) . 

In  this  report  we  equate  the  integrals  over  the  given 
downwash  over  certain  areas  of  the  wing  with  the  same  integral 
over  the  downwash  expressed  in  terms  of  the  parameters  for  the 
pressures.  These  areas  will  be  called  control  areas. 

As  the  work  progressed  it  became  more  and  more  apparent  that 
analytical  integrations  would  play  an  important  role.  The  kernel 
of  the  integral  equation  is  accompanied  by  a  factor 
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exp  ik(^-x).  We  set 


Ap(5in)  =■  Ap(  5,  n)exp(  ikC)  (6) 

This  is  one  of  the  steps  which  make  an  analytical  integration 
feasible.  Moreover,  we  introduce  a  weight  function  exp(ikx)  in 
the  integrals  for  the  upwash.  Since  the  reduced  frequency  is  not 
large,  expdk^)  and  exp(ikx)  vary  only  by  a  small  amount  within  an 
element;  or  even  within  the  pressure  and  control  areas. 

During  the  course  of  the  work  the  author's  idea  about  the 
choice  of  the  pressure  and  of  the  control  elements  have  undergone 
changes.  First  he  had  a  subdivision  of  the  planform  in  mind  which 
is  suggested  by  the  vortex-lattice  method,  and  which  has  been  used 
in  the  work  of  Ueda  and  Dowell  (Fig.  1).  Then  he  realized  that 
for  such  trapezoidal  elements  it  is  not  possible  to  find  elemental 
pressure  shape  functions  which  satisfy  the  two  requirements  that 
the  pressures  be  continuous  as  one  passes  from  one  element  to  the 
next  and  that  the  resulting  integrals  can  be  integrated 
analytically.  Triangular  elements  with  linear  shape  functions  for 
Ap(^,Ti)  are  preferable.  Such  elements  can  be  obtained  by  drawing 
one  diagonal  into  the  trapezoids  of  Fig.  1  (see  Fig.  2). 

In  the  triangles  in  Fig.  2  one  side  is  parallel  to  the  x- 
axis  (to  the  direction  of  the  wake  streamlnes).  Consider  a  single 
element  and  assume  that  Ap(^,r))  is  constant  (Fig.  3a).  Then  one 
obtains  at  the  wake  boundary  n  •  a  singularity  in  the  upwash  as 
(n^'y)  »  element  boundary  and  at  n  »  a  singularity  as 

log(n'j“y).  In  contrast,  if  none  of  the  sides  of  the  triangle  is 
parallel  to  the  x-axis,  then  one  obtains  in  the  wake  three 
singular  lines  n  -  Pi  .  n  •  n2f  n  •  but  only  with 
logarithmic  singularities.  In  other  words,  one  obtains  a  smoother 
upwash  although  infinities  are  still  present.  This  is  of 
particular  interest  if  one  works  with  elements  of  constant  &p.  If 
the  pressure  elements  have  sides  parallel  to  the  x-axis,  then  the 
control  elements  must  overlap  lines  of  the  wake  pertaining  to  the 
x-axis,  otherwise  integrals  over  the  upwash  will  be  infinite.  For 
pressure  elements  with  no  side  parallel  to  the  x-axis  this  is  not 
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necessary,  as  the  logarithmic  infinities  give  a  contribution  which 
remains  finite  after  the  integration.  Incidentally,  for  constant 
pressure  elements  one  can  use  for  the  elemental  areas  quadrangles 
as  well  as  triangles.  Such  an  elemental  subdivision  is  shown  in 
Fig.  (Special  measures  may  be  needed  at  the  wing  tips  and  at 

the  center  line.)  A  subdivision  in  triangles  is  obtained  by 
drawing  into  the  quadrangles  the  diagonal  which  is  not  parallel  to 
the  x-axis.  Originally,  the  author  thought  it  desirable  to  choose 
the  quadrangles  or  the  triangles  so  that  their  corners  lie  on 
lines  of  constant  y  but  actually  the  integration  formulae  are  of  a 
nature  that  a  proliferation  of  singular  lines  in  the  wake  causes 
no  additional  computational  labor. 

The  numerical  work  turns  out  to  be  quite  complicated.  It 
would  be  greatly  simplified  if  one  could  use  a  subdivision  into 
elements  which  possess  a  repetitive  pattern.  For  wing  plan  forms 
with  parallel  leading  and  trailing  edges  this  can  easily  be 
achieved  (Fig.  5).  Here  one  has  only  two  types  of  elements  (see 
Fig.  5a)  and  it  suffices  if  one  establishes  only  once  for  each 
type  the  necessary  integrals  over  the  combinations  of  one  pressure 
element  and  close  upwash  elements.  If  the  trailing  edge  is  not 
parallel  to  the  leading  edge  and  one  uses  these  subdivisions,  then 
one  must  admit  exceptional  elements  at  the  trailing  edge  (see 
Fig.  6). 

For  straight  leading  and  trailing  edges  one  can  obtain  a 
subdivision  which  at  least  has  self-similarity  (Fig.  7).  Let  a 
and  b  be  the  chords  of  the  wing  at  its  root  and  at  the  tip, 
respectively.  The  net  of  element  boundaries  to  be  drawn  is 
self-similar  with  respect  to  the  point  of  intersection  of  leading 
and  trailing  edges  (point  0).  First  we  choose  the  cornerpoints  of 
the  net  that  lie  on  the  leading  edge.  The  distances  from  point  0 
are  chosen  in  such  a  manner  that  each  point  is  mapped  into  the 
next  one  closer  to  0  if  one  multiplies  the  scale  of  the  figure  by 
a  factor  r  <  1.  If  there  are  n  ♦  1  points  along  the  leading  edge 
(including  the  points  at  the  wing  root  and  at  the  wing  tip),  then 
one  needs  n  such  transformations  to  transform  the  wing  root  a  into 
the  wing  tip  b.  Therefore, 
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If  S  is  the  wing  span  (measured  in  the  y-direction)  then  point  0 
lies  at  a  distance 


S,  -  S  a(a-b) 

The  distance  of  the  point  along  the  leading  edge  from  point  0, 
(if  one  counts  from  the  root  outward  and  assigns  m  -  0  to  the 
root)  is  then 


It  gives  some  simplification  if  the  sides  of  the  elements  arrange 
themselves  in  straight  lines.  The  individual  triangular  elements 
are  then  embedded  in  larger  triangles.  Such  a  net  is  obtained  in 
the  following  manner.  First  one  chooses  one  of  these  larger 
triangles  with  two  corners  at  the  points  of  the  leading  edge 
determined  above  and  one  corner  at  the  trailing  edge.  Next  one 
draws  straight  lines  parallel  to  the  sides  of  these  triangles 
through  all  those  points  of  the  leading  edge  (including  some 
outside  of  the  wing).  This  divides  the  plan  form  into 
quadrangles.  Because  of  the  choice  of  the  initial  points  at  the 
leading  edge  one  set  of  diagonals  in  these  quadrangles  form 
straight  lines  through  point  0.  By  drawing  these  diagonals  one 
obtains  the  desires  self-similar  net.  In  general,  the  other 
diagonals  will  not  form  a  straight  line  and,  therefore,  the 
cornerpoints  lying  on  a  tract  of  these  diagonals  will  not  form  a 
line  y  -  const.  Of  course  if  one  of  the  sides  of  the  original 
large  triangle  is  line  y  -  const,  then  all  triangles  will  have 
such  a  side.  For  elements  of  constant  Ap  this  complicates  the 
choice  of  upwash  areas.  This  does  not  happen  for  elements  in 
which  Ap  is  linear. 

The  parameter  describing  the  pressure  distribution  are  the 
values  of  Ap(t.n)  -  Ap( 5, n)exp( ix^)  at  the  corners  of  the  grid 
formed  by  element  boundaries.  One  parameter  therefore  generates 
pressure  distributions  in  all  elements  that  contain  the  corner 


which  pertains  to  the  parameter  whose  influence  one  wants  to 
compute.  The  region  covered  by  the  elements  pertaining  to  a 
certain  pressure  parameter  will  be  called  the  pressure  area.  In 
the  interior  such  an  area  consists  of  six  triangular  elements  (see 
Fig.  8).  The  integration  of  the  upwash  (with  a  weight  function 
exp(ikx)  is  carried  out  over  "upwash  areas."  They  are  identical 
with  the  pressure  areas.  This  ensures  that  the  number  of 
equations  equals  the  number  of  unknowns.  The  elements  of  the 
matrix  of  the  system  for  the  determination  of  the  pressure 
distribution  are  given  by  the  upwash  due  to  the  pressure  in  an 
individual  pressure  area  (with  the  value  1  assigned  to  the 
pertinent  pressure  parameter)  integrated  over  the  upwash  areas. 
(This  is,  of  course,  nothing  new,  the  idea  is  inherent  in  all 
methods.)  Each  element  belongs  to  several  pressure  and  upwash 
areas.  To  each  triangular  element  belong  three  "elemental" 
pressure  distributions.  The  primary  task  is  the  determination  of 
the  upwash  due  to  the  elemental  pressure  distribution  integrated 
over  the  different  elements. 

For  elements  with  constant  &p,  in  which  no  side  is  parallel 
to  the  x-axis,  the  elements  themselves  can  be  used  as  pressure  and 
upwash  areas.  If  element  sides  are  parallel  to  the  x-axls  and  Ap 
is  constant,  then  the  elements  still  constitute  the  pressure 
areas,  but  the  upwash  areas  must  be  chosen  so  that  one  has  an 
overlap  of  the  singular  lines  within  the  wake;  this  is  necessary 
in  order  to  avoid  infinities  in  the  integrated  upwash. 

Let  N  be  the  number  of  pressure  parameters.  It  is  the 
number  of  corners  in  the  grid  at  which  the  pressure  is  unknown, 
i.e.,  the  number  of  grid  corners  except  for  those  at  the  trailing 
edge  where  the  pressure  difference  is  zero.  It  is  also  the  number 
of  pressure  and  of  upwash  areas.  The  pressure  parameters  (each 

with  a  pertinent  pressure  area)  and  also  the  upwash  areas  are 

2 

numbered  from  1  to  N.  Ultimately,  one  will  generate  an  N  matrix. 
Further  numberings  are  introduced  for  the  elements  and  for  the 
elemental  pressure  distributions.  The  number  of  elemental 
pressure  distributions  is  somewhat  smaller  than  three  times  the 
number  of  elements,  because  the  pressures  at  the  trailing  edge  are 


zero.  The  numberings  of  the  pressure  areas  and  of  the  upwash 

areas  are  the  same.  The  numberings  of  the  elemental  pressure 

distributions  and  of  the  upwash  elements  (with  single  subscripts) 

are  carried  out  independently.  Let  Ng  be  the  number  of  upwash 

elements  and  Np  the  number  of  elemental  pressure  distributions. 

Which  elemental  pressure  distributions  belong  to  certain  pressure 

areas  and  which  elements  belong  to  the  upwash  areas  is  shown  by 

"housekeeping"  matrices  and  of  dimension  Np  by  N  and  N 

by  Np,  respectively.  The  elements  of  these  matrices  are  zero  and 

( 1 ) 

one.  An  element  is  one  If  the  elemental  pressure  distribu- 

bution  with  index  k  pertains  to  the  pressure  area  with  index  l. 

Each  row  of  contains  only  a  single  one,  each  column  a  maximum 

( ) 

of  six.  The  element  is  one  if  the  elements  with  index  j 

belong  to  the  upwash  area  with  index  i. 

(2)  (2) 

An  element  M  .  of  a  third  matrix  M  gives  the  upwash 

J  I  K 

integrated  over  the  element  with  index  j  due  to  the  elemental 
pressure  distribution  with  index  k.  The  matrix  for  the  system  of 
equations  from  which  the  pressure  parameters  are  determined  is 
denoted  by  M.  One  has 


M 


(7) 


(2) 

The  main  effort  is  the  determination  of  the  matrix  M  .  An 

element  in  which  the  pressure  is  prescribed  will  be  called  ^n~ 

element  (because  of  the  independent  variables  for  the  pressure), 

and  an  element  over  which  the  upwash  is  integrated  will  be  called 

xy-element.  In  Section  III  formulae  will  be  developed  for 
(2) 

elements  of  M  for  which  the  xy-element  is  at  some  distance  from 
the  5n“element.  In  subsequent  sections  xy-elements  close  to  ^n- 
elements  will  be  treated. 


SECTION  III 

EVALUATION  OF  THE  ELEMENTS  OF  THE  MATRIX  M^^^  FOR 
xy-ELEMENTS  AT  A  DISTANCE  FROM  A  ^n-ELEMENT 


The  computation  is  based  on  Eq.  (5) 
appears  only  in  the  variable  r.  One  has 


The  z-coordinate 


3  z  3 
3z  “  r  3r 


If  the  xy~element  is  not  adjacent  to  the  5n-element  or  its  wake, 
then  the  function  K  is  free  of  singularities.  One  has 


1^  (zK)  -  K  *  (z'^/r)(3K/3r) 


lira  (|r(zK))  -  lira  K 
z-0  z-0 

One  simply  replaces  |r|  by  ln-y|.  This  is  the  approach  used 
without  restriction  in  Reference  2. 

The  specific  form  of  K  is  found  in  Appendix  B, 


K  -  6  exp(ik(^-x) ) (K^  +  K2) 


No  wake  singularities  occur  in  K2.  The  wake  singularities 
occurring  in  are  displayed  in  Eq.  (B.64a). 

We  have  introduced  in  Eq.  (6) 


Ap  (5,n)  -  exp(ik5) Ap(?;,Ti) 
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Eq.  (5)  will  be  discretized  by  postulating 


B8it  II  exp(  ikx)w(x,y)dxdy 
Aj(x,y) 

(10) 

-  II  z  II  Ap(C,n)(K^  +  K2)d5dn}dxdy 

Aj(x,y)  Aj^(5,n) 

Here  A^(x,y)  and  Aj^C^.n)  refer  respectively  to  the  i'^  upwash  and 
the  %  pressure  area.  For  a  triangular  element  with  corners 
numbered  j,  k,  I  and  pressure  1  at  corner  j,  and  pressure  0  at 
corner  k  and  I,  the  linear  pressure  shape  function  is  given  by 


Ap  (Sih) 


(6-5j) 

(n-hj) 

(n„-nj)  (n,-nj) 


(11) 


This  is  easily  verified  by  simple  operations  on  determinants  if 
one  sets  C  -  5.  ,  n  -  ht,.  and  C  -  Co »  h  -  n.  •  matrix  elements 

is  obtained  by  replacing  in  Eq.  (10)  A. (x,y)  and  Ao(C.n) 

Vh  I'H  ^  ^ 

respectively  by  the  i'^  and  the  surface  elements. 

If  the  xy-element  is  not  adjacent  to  the  Cn-element  or  Its 
wake,  then  the  integrand  of  Eq.  (10)  is  analytic  in  n,  x,  and  y 

and  the  integrations  can  be  carried  out  with  efficient  numerical 
methods.  There  exist  even  formulae  for  Gaussian  integration  over 
a  triangular  area  (developed  for  applications  in  elasticity). 

If  the  xy-element  is  close  to  the  wake  but  at  some  distance 

from  the  Cn^element,  one  proceeds  as  follows.  The  function 
(see  Eq.  (B.64a))  is  written 
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Ku  -  -B^CR(R-(5-x))]‘^-(k^/2)log(2k‘‘' (R-(5-x)))  j" 


-  (n-y)^  ♦ 


The  function  Kg  and  are  analytic  even  in  the  wake.  Their 
contributions  can  again  be  treated  by  Gaussian  integration. 

K3  “  ^31  *  *^32  upon  r  only.  For  the  evaluation  of  this 

part  one  divides  both  the  Cn-eleraent  and  the  xy-element 
respectively  by  a  line  n-const  and  a  line  y-const  into  smaller 
triangles  (see  Fig.  9).  The  contribution  of  these  smaller 
triangles  are  treated  separately.  &p  is  a  linear  function  of  5 
and  n 


Ap(C.n) 


+  5  ♦  Cg  n 


Consider  a  fixed  and  a  fixed  xy  triangle.  Assign  the  index  1 
to  the  corner  opposite  to  the  side  parallel  to  the  C“axis,  and 
subscripts  2  and  3  to  the  other  corners  (proceeding  in 
counterclockwise  direction).  Figs.  10.  We  write 

Ap  -  c«  c,  (e-E,)  +  c^(n-n, )  (17 


here 


^0  “  ^0  h  *  ^2  ’’l 

The  sides  here  1  2  and  1  3  of  the  triangle  are  given  respectively 
by 


and 


^2~^1 


{ n-n^ ) 


Then  one  can  carry  out  the  Integration  over  ^  at  constant  ni  and 
the  integrations  over  x  at  constant  y.  Let,  for  Fig.  (10a) 


f2(h)  - 
♦  (C^ 


f,(y) 

*  ^2 

/2)(n-ni)^ 


-  )(x2-x2)/(y2-y^ ) 

(n-n^  )  ^ 


n,) 


(18) 


(19) 


The  upper  and  lo^'er  signs  hold  for  Figs.  10(a)  and  10(b), 
respectively.  There  are  three  elemental  pressure  distributions 
Ap(5,n)  in  each  element,  and  therefore  three  different  sets  of 
constants  C^,  and  three  different  functions  They  are 

quadratic  in  n« 

Let  be  the  integral  on  the  right  of  Eq.  (10)  with 
+  Kg  replaced  by  .  If  the  triangles  in  the  x,y-plane  and 
in  the  E,n-plane  have  the  orientations  of  Fig.  10(a),  one  has 


and  need  not  be  the  same  as  and  because  the  x,y- 
element  need  not  lie  in  the  wake  of  the  5h-element.  The  inner 
integration  is,  of  course,  carried  out  at  constant  y.  We  write 

f2(n)  -  f2(y)  f2(y)(h-y)  +  f2(n,y)(n-y)^  (2C 


where 


fo(n,y)  -  [fo(h)  -  fo(y)  -  fA(y)(n-y)]/(n-y)' 


The  function  f^Cn.y)  is  analytic  at  y  -  n.  Then  appears 


in  the  form 


^31  “  ^311  *  ^312  *  ^313 


where 


I3,,  -  8^  j  f,(y)r2(y)| 


(n-y)^  ♦ 


B  I  fi(y)f2 
yi 


/^3 

5(y)(  J 

\n, 


(n-y )dn 
(n-y)^  + 


I-,,-.  -  f  f,(y)  f  fo(n,y)dn  -  z  f 


2  "3 


dn  dy  (2U) 


It  i3  Shown  in  Appendix  D  that  11m  4-  (...)  applied  to  the  second 

z-»0 

term  in  the  bracket  of  I313  is  zero.  The  first  term  is 
independent  of  2.  Thus, 


lira  ^  (ZJ313) 


’’'K 


j  f3(n,y>dn)dy 


The  integrands  are  analytic  functions. 

The  inner  integration  In  is  carried  out  analytically. 


(n-y)  +  z‘ 


1  arctg 
z  z 


One  then  has  to  evaluate 


lira  Iz  z  ^farctg  ^  ] 
z-0  ^  2 


(n-y)  + 


nyy  n^-y 

This  result  could  have  been  obtained  by  setting  z  -  0  on  the  left 
side  of  Eq.  (25),  and  substituting  into  the  formal  expression  for 
the  indefinite  integral  (namely  -n  M  the  limits  n.|  and  n^;  this, 
in  spite  of  the  fact  that  the  Integral  obtained  by  immediately 
setting  z  -  0  does  not  exist  for  n  -  0.  Then 


■52  ^311^ 


-  I  f^(y)f2(y)  (  - 


n^-y  n^-y 
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The  integrand  is  singular  for  y  -  and  y 


If  the  xy- 


eleraent  lies  exactly  in  the  wake  of  the  5n-element,  then  » 

^3  ”  ^3’  integrand  is  singular  at  the  two  limits.  If  the  xy- 
element  is  adjacent  to  the  wake,  then  this  happens  at  only  one 
limit.  The  two  singular  terms  are  now  treated  separately. 


^3  fi(y)f2(y) 


n^-y 


-  f^(n3)f2( 


"3>  I 


(26a) 


^3  [’f^(y)f2(y)  ~  f^  (n3)f2(n3)1 

I  L  ho-y 


The  second  term  on  the  right  is  a  smooth  function.  It,  therefore, 
can  be  integrated  numerically  without  difficulty.  It  is  best  to 
write  it  in  the  form 


^3  rfi(y)  -  ^1(13) 


"  1  L' 


)  -  f. (n,)  fp(y)  -  fpCn,) 


The  first  term  on  the  right  of  Eq.  (26a)  gives 


f,(n3)fj(n3)  log  1^^ 


This  term  is  infinite  if  y3  -  n3.  This  happens  because  at  the 
element  boundaries  the  upwash  caused  by  the  tn-element  behaves  as 
(h3-y)  If  Ap  is  continuous  as  one  passes  over  the  element 
boundaries  in  the  n“direction  (as  it  happens  for  the  two 
boundaries  between  the  smaller  triangles  in  Fig.  9  or  for  the 
triangular  elements  with  one  side  parallel  to  the  ^-axis  and 
linear  pressure  distributions),  then  the  singular  term  is  canceled 
by  a  contribution  of  the  adjacent  5n-element. 


We  have 


“S  32  ^31 1> 


ot  h.3~yo  h.  “3^0 

B  (n3)f2(n2)iog  ^  (n^  )f2(n^  )log 


^3  fi(y)f2(y)  -  fi(n3)f2(n3)  ^  ^3  fi(y)f2(y^  ‘  fi  (n,  )f2(nT ) 

J  n^^^^y  ^  J 


n^-y 


(27 


The  inner  integral  of  12^2  Slves 


^  log( ( n-y ♦  z^) 


One  has 


(1/2)lim4v  (z  log((Ti-y)^  +  z^))j 

z*0  I 


-  lim  (1 /2)log( (n-y ♦  z^)  ♦ 


z*0 


z2 

>2  .'-S 


_  (n-y)  +  z  „ 

Hi 


I  J  -  log|n-y|  I 


This  result  would  have  been  obtained,  if  one  sets  in  Eq.  (23) 
z  -  0,  disregards  the  singularity  of  the  integrand,  integrates 
analytically,  and  substitutes  the  limits. 


If  lies  within  the  interval  of  integration,  one  chooses  q  - 
This  ensures  that  (f2(y)  ~  is  an  analytic  function  of 

y,  even  at  y  -  n^*  *13  ii®s  outside  of  this  region,  then  the 

integrand  is  analytic  in  any  case.  It  is,  however,  still 
advisable  to  choose  the  constant  f^(q)  so  that  the  Integrand  is  as 
smooth  as  possible,  but  it  is  not  necessary  that  one  evaluate  Eq. 
(29)  for  values  of  y  outside  of  y^  <  y  <  y^t  or  if  one  does,  it 
need  not  be  done  with  precision.  Therefore, 


I312  ”  “  f4(n3)  )log(n3-y)  ♦  J 


^  ^3  fii(y)  -  f4(n3) 


n,-y 


-  (f4(y)  -  f4(n^ )log(n^ -y)  -  J 


^  ^3  f4(y)  -  f4(ni) 


ni-y 


Let  1^2  be  the  Integral  I  with  K  replaced  by  Eq.  (1M).  We 

omit  a  demonstration  which  would  show  that 


lim(|-(zl,~)  -  lim  I,, 
z*0  z*0 


fc(n-y)  -  I 
^  n-0 


(k(n-y)/2e)‘ 
n! (n+l) ! 


Then 


lim  I. 


^3  /'’3  \ 

I  I  log( n-y )f5(Ti-y)f2(n)dnjdy 


fgCn.y)  -  i  f2(v)f^( v-y)dv 


The  inner  integral  then  becomes 


j  log| n-y  I  Ofg/3n)dn  -  (iog|n-y|fg(n,y)  j  -  f^(y) 
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with 


(y)  -  I 


fg(n,y) 


The  Integrand  in  the  last  expression  is  a  smooth  function  of  n 
which  can  be  treated  numerically. 

For  the  evaluation  of  1^2*  one  then  needs 


{k^/2)  J  f,(y)f7(y)dy 


k^/2  I  f^ (y)Clogln3-ylfg(n3,y)  -  logln^-ylfgCn^ ,y)]dy  (38) 

The  integrand  of  Eq.  (37)  is  smooth.  The  two  terms  in  the  bracket 
of  the  integrand  in  Eq.  (38)  are  treated  separately.  Let 


fg(yfTi3)  -  I  (v)f^(Ti2,v)dv 


Then  one  obtains  instead  of  the  first  of  these  two  terms 


/  Iz  "  ^8^<l'^3^]^o8ln3-yldy 


“  ^og j n^ “ y I C ( y » n3“ t Q » n3 ) I  ~  f^(n3) 
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where 


^3  fQ(n,no)  *  fo(q.no) 

-  I  - -T^ - ^ 


(41  ) 


If  <112^^3*  will  choose  q  -  n^*  Then  the  integrand  of 
Eq.  (41)  is  an  analytic  function,  even  at  y  -  n^*  ^3  I® 

outside  of  this  region  then  the  integrand  of  Eq.  (41 )  is  analytic 
for  any  choice  of  q  in  the  constant  f^(q,n2)»  but  one  will  choose 
q  approximately  equal  to  to  obtain  a  smooth  integrand  of 
Eq.  (41 )  in  the  interval  of  integration.  The  same  procedure  is 
applied  to  the  other  term,  with  replaced  by  • 

Then 


'32 


(k^/2) 


-j)  f,(y)f7 


(y)dy 


♦  log(n2”y)Cf8^y*b2)  "  j 


“  f q ( ^ 


(42) 


-  log(n^ -y)CfQ(y,n^ )  -  fQ(q,n,)] 


♦  fq(n^ ) 


For  triangular  elements  the  functions  f^  and  (Eqs,  (19)) 
are  given  by 


f^  (y.y,  ,c^ )  -  c,  (y-y, ) 


1  ■'I 


f2Chin,  ,c^  ,0^)  -  C2(n-n^)  ♦  C2(n-np' 
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The  constants  ,  0^$  and  are  the  only  parameters  by  which  the 
elemental  pressure  distributions  enter  the  computations.  It  is 
therefore  practical  to  program  their  contributions  to  the 
following  expressions  separately  then 


I'.* 


^3  ■  °3 


and 


^31  3  "■  °1°3^^3~^1  ^ 


One  has 


-1 


(f^(y)f2(y)  "  f^ (v)f2(v) ) (v-y) 

-  [(f,(y)  -  f^ (v))/(v-y)]f2(y)  ♦  [(fgCy)  -  f2Cv))/(v-y)]f^ (v) 

-  -c^{c2(y-np  ♦  c-Cy-n^)^  ♦  (v-y^  )[C2  ♦  c-Cy-n^  ♦  v-n^]) 


Then 


> 

I 

(T, 


|(f^(y)f2(y)  -  f^ (v)f2(v) (v-y)  ^dy 
-  -c^  C2C(v-y)^  (y-n, )  ♦  (y-n^)^/2] 

-c^  CjCCv-y^  (v-n^ )  (y-n^ )  -  ((v-y^  )(y-n^  )^/2)  +  (y-ni)^/3] 

In  the  program  this  is  considered  as  a  function  of  y,  v,  y^ ,  i 
(0^02),  and  (c^c^).  The  expression  occurring  in  Eq.  (26a)  are 
then  obtained  by  substituting  y^  or  y^  and  for  y  and  v 

respectively.  The  lower  limit  in  the  integral  for  f^  does  not 
matter. 

The  integrand  in  fi^  is  given  by 

0,(v-y,)[02  .  203  (y-y,  .  y,  -  i,,)3 


2iJ 


Then 


fijCy)  -  c^CgCCy-y^  )^/2)  ♦  2c^c2[((yi-ni  Xy-y^  )^/2)  +  ((y-yi)^/3)] 

programmed  as  a  function  of  y,  ,  (c^C2),  and  (c^c^)- 

|(fi,(y)  -  fij(v) )  (v-y)“’dy  -  -  (c^  C2/2)  [  (v-y^  (y-y^  ♦  (y-yp^/2)] 

-c^c^CCy^-n^  )((v-yp(y-y^  )  +  (y-y^  )^/2)  +  (2/3)  (v-n^  )^(y-yi ) 

♦  (1 /3) (v-y^ ) (y-y^ ♦  (2/9) (y-y^ )^] 

programmed  as  a  function  of  y,  v,  y^  ,  n.j .  (0^02).  and  (c^c^).  The 
integral  on  the  right  of  Eq.  (30)  is  then  obtained  by  substituting 
y^  and  y^  for  y,  and  for  v.  The  same  procedure  is  carried  out 
for  the  second  term  in  the  integral  on  the  right  of  Eq.  (28);  one 
simply  replaces  by  . 

2 

The  terms  of  1^2  are  0(k  ),  moreover  n-y  is  small  because 
these  computations  will  be  carried  out  only  if  the  xy-elements  are 
in  the  vicinity  of  the  wake  of  the  Cn-eleraent.  Therefore,  an 
approximation  of  f^  by  its  first  term  will  suffice 


Then 

v-n 

fg(n,y)  -  j  f2(v)dv 
v-y 

-  ((c2/2)(n-rii  )^  -  (y-ni)^)  ♦  ((C3/3)  (n-Hi  )^  -  (y-n,)^) 

For  the  evaluation  of  f0  this  expression  is  reordered  in  powers  of 
(y-Yi) 


fg(n,y)  -  (c2/2)[(T)-n, )  -  "  2(y,-ni )  (y-y^ )  -  (y-yi )  3 

♦  (c2/3)C(n-n,  -  (y^-ni)^  -  SCy^-n^  )^(y-yi ) 

-  SCy^TTii )  (y-y^  -  (y-y^  )^] 

For  the  evaluation  of  the  first  form  of  fg  is  used 

'’3  _1 

fyCy)  -  j  fg(n,y)(n-y)  dy 

^1 

fyCy)  -  (c2/2)[((n3-Tii  )^/2)  +  (y-ni  )(n3-ni )] 

>  (c3/3)[((n3-ni)^/3)  ♦  ( (y-n^  )  {n3-n,  )^/2)  ♦  (y-n^  )^(  n3-ni )  ] 

This  is  reorderd,  in  powers  of  (y-y^ ) 

fj(y)  -  (c2/2)[((n3-ni  )^/2)  +  (y^ -hi ) (n3-ni )  ♦  (n3-n, )  (y-y, )] 

>  (C3/3)  {C((n3-ni  )^/3)  ♦  ( (y^ -n, ) (n3-ni  )^/2)  *  (y^ -hi  )^(n3-n^ )] 

♦  C((n3-n^  )^/2)+  2(y^-n,  )(ri3-n^  )](y-yi )  ♦  ( n3-ni ) (y-y^  )^ } 

Then 

^3 

I  f^(y)fy(y)dy  -  (c^ C2/2)  ([ ( ( Ti3-n^  )^/2) 

2  3  , 

+  (y^ -n^ )  (n3-n^ ) ](y3-y^ )  /2  ♦  (n3-n^  )(y3"yi )  /3] 

♦  (c^ C3/3)  IC ( ( n3-n^  )^/3)  ♦  ( (y-n^ )  (n3-ni  )^/2) 

♦  (ni-yi)^(n3-ni)](y3-yf)/2  ♦  L((n3-ni  )^/2) 

♦  2(y.-ni  )(n,-ni  )]((y,-yi  )^/3)  ♦  ( Oo-n, ) ( (y-y,  )^/'«) ) 


fgCYfn)  -  I  (v)fg(n,v)dv 

-  (c,c2/2)i[(n-ni)^  -  (yi-ni)^]((y-yi)^/2)  -  2(yi-ni )  ((y-y^  )^/3) 

-  (y-y,  +  (c^c^/s)  ([(n-n^)^  -  (y^-n^  )^]( (y-yi  )^/2) 

-  3(y,-ni)^((y-yi  )^/3)  -  3(yi-ni  XCy-y,  )V»»)  -  ((y-y^)^/5)} 

This  Is  programmed  as  a  function  with  arguments  y-y^  ,  rrHit  yi'n^i 
(0^02),  and  (c^c^). 

Finally,  one  needs 


f9(niq)  -  j 


3  fg(y,n)  -  fgCq.n) 


y-q 


dy 


but  only  for  n 


Ultimately,  one  must  substitute  q 


'1  • 


f9(qfq)  -  (0^02/2)  {[(q-y^  -  (y^-n^)^][(y^-y^)^/^ 

*  ((q-y^  )(y2-y^  )/2) ]  -  2(y^-n^  )[((y2-y^  )^/9)  +  (q-y^  )((y3-yT  )^/6) 

*  (q-y/)((y3-yi  )/3)]  -  C((y3-yi )  V16)  (q-y^ ) ( (y3-yi  )^/i2) 

+  (q-y^  )^((y-yp^/8)  +  (q-y^  )^( (y3-y^  )/4) ]} 

*  (c^c2/3)  ([(q-n^  -  (y^-n,  )^][((y3-yi  )^/4)  *  (q-y^  ) ( (y3-yi  )/2) ) ] 

-  3(y^ -n^  )^[( (y3-y^  )^/9)  +  (q-y^  )((y3-y^  )^/6)  +  (q-y^  )^(  (y3-y^  )/3) ] 

-  3(y^-n^  )[((y3-yi  ♦  (q-yi  )((y3-yi  2)  ♦  (q-y-,  )^( (y3-yi  )^/8) 

*  (q-yi)^((y3-yi  )/4)]  -  [( (y3-yp^/25)  +  (q-y, )  ( (y3-y,  )  V20) 

+  (q-y,  )^  ( (yo-yi  )^/l  5)  >  (q-y,  )^(  {y^-y,  )^/l  0)  +  (q-y,  )^(  (y^-y,  )/5)  ]} 


SECTION  IV 


THE  LIMITING  PROCESS  z-^-O  FOR  POINTS  (x,y) 
CLOSE  TO  OR  WITHIN  A  Cn-ELEMENT 


A  more  refined  procedure  is  needed  for  points  (x,y)  close  to 
points  (t,n)  because  certain  terms  in  K  have  a  denominator  R, 
which  tends  to  zero  if  (^-x)  and  (n~y)  tend  to  zero  simul¬ 
taneously,  and  because  for  ^-x  negative  there  are  terms  with  a 

2  2  1/2 

denominator  r  -  ((n-y)  +  z  )  .  The  latter  singularity  has  been 

treated  in  a  simpler  situation  in  Section  III.  In  the  present 
section  we  study  the  limit  z+O  for  the  upwash  at  a  fixed  point 
(x.y) . 

According  to  Eqs.  (B.6U)  one  must  evaluate 


w(x,y,z  -  0) 


z-0 


jj  Ap(e.n)(K^ 
A(  c.  n) 


K2)d^dn 


(43) 


where 


w(x,y,z)  -  8TrBexp(ikx)w(x,y ,z) 


The  function  K  -  +  K2  depends  only  upon  ^-x  and  r  with 


r 


( ( n-y)^ 


2,1  /2 

+  z  ) 


The  independent  variable  z  occurs  only  in  r.  Because  of  the 
singularities  in  K  and  because  of  the  differentiation  with  respect 
to  z  the  result  of  the  limiting  process  z*G  is  not  entirely  self- 
evident.  For  part  of  the  discussion  we  introduce  polar 
coordinates 


5-x  -  pcosa,  n-y  -  psina 


One  deals  with  Integrals 


jj  F(p,a)dpda 

Expressing  F(p,o)  in  the  form 

F(p,a)  -  3G(p,o)/3a  (45) 

one  obtains  by  an  argument  familiar  from  the  derivation  of  Green's 
theorem 

jj  F(p,a)dpda  -  jj  ( dG(  p,  a)/3a}dpda  -  -^G(p,a)dp  (46) 

A  A 

Alternatively  one  can  set 

F(p,o)  -  3H(p,a)/3p  (47) 


Then  one  obtains 

jj  F(p,a)dpda  -  J|(3H(  p,a)/3p)dpda  -  ^H(p,a)da  (48) 

A  A^ 

If  the  point  xy  lies  outside  the  ^n-element,  then  p  ^  0 
throughout  the  element,  and  a  returns  after  one  complete  circuit 
around  the  element  to  its  original  value  (Fig.  11).  If  it  lies 
inside,  then  p  varies  between  0  and  p(a),  and  a  between  0  and  2it 
(Fig.  12). 

For  part  of  the  discussion  it  is  convenient  to  write 


P  ^  _ 

H(p,a)  -  I  F(p,o)dp 

0 


O  O 


(where  p  is  a  dummy  variable  of  integration).  Then 


L  p(a)  ^ 

JJ  F(p,a)dpda  “  J  (  |  F( p, a)dp j(da/dt) )di 
A  t-0  '  0  * 

Here  I  is  the  length  measured  along  the  contour  of  the  element  (or 
some  other  parameter  which  varies  monotomically  from  0  to  L).  The 
value  of  a  along  the  contour  is  considered  as  function  of  1. 

We  start  the  discussion  of  the  individual  terms  of  K  with 
Kg,  Eq.  (B.65).  Only  the  first  term,  -R  ^[(exp(ikV)  -  1)/V],  has 
a  denominator  which  tends  to  zero  as  p  and  z  tend  to  zero,  namely 
R  \  The  expression  V  has  the  same  property,  but  (exp(ikV)  -  1)/V 
is  obviously  regular  at  V  -  0.  If  any  term  is  critical  in  the 
limiting  process  z+0,  then  it  is  the  first  one  namely  -(ik)/R.  It 
is  discussed  presently.  Except  for  factors  one  deals  with  the 
expression 


L  p(a) 

|jAp(p,a)R  %dpda  -  J  j  (Ap(p,a)(p^  ♦  z^)  ^  ^^pdp)  (da/dt)  t  (49) 

0  p-0 

For  fixed  a  the  inner  integral  has  the  form 


T  \f  2  2. -1/2  . 

I  -  j  f(p)(p  +  z  )  pdp 
0 


Ultimately  one  has  to  form  lim  (z  ...)  for  this  expression. 

z-^0 

an  integration  by  parts  one  obtains 

p(a) 


By 


I  -  f(p)(p2  +  z2)’/2 


p(  o)  2  2  1/  2 

j  (df/dp)(p‘^  +  z^)^''^dp 


and 


rsr  'r.  r;  ■ 


^  (ZI)  -  f(p)(p2  .  z2)1/2  - 


(df/dp)(p^  +  z^)^^^dp 


p(a)  /V 

♦  z^[f(p)(p^  ♦  -  I  (df/dp)(p^  +  z^)“^^^dp] 


The  limiting  process  z*Q  gives 


a  p(a) 

lim(-|— (zl))  -  f(p(o))p(o))  -  [  (df/dp)pdp 


?  -1  ?  pC“' 

+  lim[z^f(p(a))p(a)  -  zf(0)  -  z^  j 
•»*n  J 


(df/dp)(p2  +  z2)"^^2dp] 


In  the  first  line  one  carries  out  an  integration  by  part  and 
obtains 

p(o) 

I  f(p)dp 
0 

The  first  two  terms  of  the  second  line  vanish  for  z  -  0.  In  the 
last  term  one  has 


p(a) 


I  I  (df/dp)(p^  +  z^)~’^^dp|  <  z"^  max|df/dp|  J  (p'^  +  z‘^)“''‘^dp 
0  0 


2  .  2. -1/2, 


,  z2)-'^=dp  -  log(p  .  (p2  » 


p  (o) 


The  lower  limit  might  be  critical;  it  gives  log  z.  The  entire 

2 

expression  is,  however,  0(z  logz). 

One  thus  obtains 


.  PU 

I 


f(p)(p^  +  z^)" 


p(a) 

I  f(p)dp 


Applying  this  result  to  Eq.  (49)  one  has 


-1  ^  -  \ 
fj  Ap(p,a)R  pdpda  "  f  (  J  (Ap{p,a)dpJ 

0  ^p-0  ' 


(da/dt)dt 


The  results  for  K,  are  less  foreseeable.  As  p+0,  the  first 

’  -2 

term  of  behaves  for  z  -  0  as  p  ,  moreover,  there  is  a  singular 
point  in  the  wake  of  the  point  (^,n),  even  for  p  0 . 

For  a  fixed  point  (x,y)  the  function  Ap(5,n)  is  written  in 
the  form 

Ap(C,Ti)  ~  I  I  c  (x,y){5-x)‘”(n-y)'^  (51) 

m  n  mn 


obviously 

Cqq  -  Ap(x,y) 

In  practice  m  <  1 ,  n  <  1 ,  but  the  analysis  will  be  carried  out  in 
more  general  terms.  The  individual  terms  in  Eq.  (51)  will  be 
treated  separately. 

2  -1 

The  term  3  [R(R  ♦  5-x)]  of  gives  a  contribution  to  the 
upwash  of  the  form 


3  J  (e-x)™(ii-y)  CR(R  +  C-x)]"  d^dn 


•w  *  •.  ^ 


To  thi3  expression  we  apply  the  procedure  implied  by  Eq.  (45), 
(rather  than  Eq.  (47)  because  the  integration  with  respect  to  a  is 
somewhat  more  easily  discussed  than  the  integration  with  respect 
to  p).  Setting 

♦  -  o  -  IT  (53) 

(because  then  an  important  singularity  occurs  at  ♦  -  0),  one 
obtains  for  a  single  term  of  the  expression  in  Eq.  (51) 

-  ,  xm+n.2/f  m+n-1  ^m  ,  ^n 

-  (■)  ®  jj  ^  ®^'''*** 

A 


[( •  (*)'"  (c 


.2  .1/2 


cos(t>  ]|  dpd(t> 


The  function  G  of  Eq.  (45)  is  then  defined  by 


0(p.,)  .  (-)""Vp'"‘"-’a-’g(a.») 


with 


g(a,<>)  -  j  cos‘”((>3in”^Ca  -  coso]'’^d(> 


One  has 


a(2,p)  -  (1  ♦  (2^/p^)^^^ 


a^  -  1  -  2^/p^ 


-1  -2 

3a/32  -  a  p  2 


If  the  point  xy  lies  outside  the  ^irelement,  then  p  >  0  for  all 
points  of  the  5n"element  and  the  denominator  p  in  Eq.  (57)  is  no 
matter  of  concern.  If  the  point  xy  lies  inside  the  Erreleraent, 


then  one  Introduces  a  cut  in  the  ^n~plane  In  the  downstream 
direction  through  the  point  xy.  In  the  plane  one  then  obtains 
the  region  shown  in  Fig.  13.  The  boundaries  AB  and  DE  are 
respectively  lines  ^  -  -ir  and  ^  ♦  ir.  Along  the  contour  of  the 

element,  p  is  considered  as  a  function  of  denoted  by  p(^). 

Along  AB  and  DE,  p  varies  from  zero  to  pC-ir)  -  p(ir). 

The  expression  in  Eq.  (5M)  is  now  evaluated  by  means  of 
Eq.  (46).  One  obtains 

„2,  .m+n+l  ,,  /a  (  m+n-1  -1  \ 

w  -  B  (-)  llm/^  z  j  p  a  g(a,(^)dp^  (59) 

Using  the  definitions  in  Eqs.  (55)  through  (57)  one  has  a  combined 
contribution  of  AB  and  DE  (Fig.  13)- 

p(ir) 

^.jm+nHg2  I  p®**^-1a'^g(a,ir)  -  g(a,--))dp 
0 

For  n  odd  the  integrand  in  Eq.  (56)  is  an  odd  function  of 
therefore. 


g(a,w)  -  g(a,-ir)  -  0  ,  n  odd 

and  one  only  needs  to  consider  the  integral  along  BCD,  which 
amounts  to  the  integral  around  the  contour  of  the  region  in  the 
plane. 


I„n  ■  (60) 

The  function  g(a,^)  is  defined  in  Eq.  (56).  The  variables  z  and  p 
enter  through  the  parameter  "a,"  Eq.  (57).  The  integral  is  formed 
along  the  contour  of  the  region. 

To  evaluate  g(a,^)  (Eq.  (56))  we  set,  for  n  odd 


g(a,((»)  -  g.  (a,(j»)  ♦  gp(a,(|.) 


(61) 


with 


g^(a,^)  -  -  f[cos*”^sin'^  ^-a^Cl  -  ^^^^][cos4>  -  a]  ^sin^d^ 

^  (62) 


ggCa,^)  -  Ja”(1  -  a^) 


2^(n  1)/2^^  _  ^sin^d^ 


In  the  last  equation,  Eq.  (58)  is  substituted  and  the  integration 
is  carried  out. 


gg  -  (-)^""^  ^'^^a®(z/p)”"^log(a  -  cos(j»)  (63) 

In  g^  one  can  carry  out  the  division  by  (cos^  -  a).  One 
obtains  a  polynomial  in  cpsip  with  coefficient  given  by  powers  of 
"a."  Notice  that 

38.J  (1  *  ♦)/ “  coa®^ain'^^(  1  -  cos^)  ^  for  n  3.  n  odd 

(6U) 

3gl(1i^)/3^  -  (cos“(>  -  1)(1  -  cos4»)~^ sin<|>  for  n  -  1 

Let 


I 

mn 


^mn,1  * 


mn,2 


(65) 


with 


.m+n+1 „2 


mn,i 


^pm^n 


-  (-)“ . 6"  r  P‘“  “  'a  ’gidp,  i  -  1  ,2 


To  obtain  the  upwash  w,  one  must  form  lim  Aloog  the 


z+0 


mn 


contour  p  /  0.  In  the  present  case  (n  odd)  one  obtains  the 
result  immediately  by  setting  a  -  1 . 


w 


mn,1  “  (-)'“^"*^  0^  j  p"'^"'^gl(1,4.)(dp/dt)dt 

l-O 


35 


where  dl  is  the  length  element  of  the  contour  and  L  the  total 
length.  Carrying  out  an  integration  by  parts  and  using  Eq.  (64) 
one  obtains 


Wmn  1  "  (- )®*”0^(m+n)“^  |  p®'^'^cos“(j>sin'^*(  1  -  cos<|>)"^  (d<|)/dl)dt 

£-0 

(66) 

n  >  3,  n  odd 


^  )■’  J  p®‘"’(cos“^  -  1)(1  -  cos*)’’sin<fr(d^/d£)d£ 

£-0 

(67) 


n  -  1 


With  Eq.  (63)  one  obtains 


,  •  6^- I  p"a"-'log(a  -  003»)dp 

tull  f  (Cm 


n-1 


Because  of  the  factor  z  ,  o  will  be  zero  for  n  >  3.  The 

mn  ,2  — 

discussion  is  therefore  restricted  to  the  case  n-1.  There  one 


obtains 


w  -  -  -  lim  -57(21-1  5) 

n ,  1  f 2  2'^0  ®1  f 2 


-  lim(-)“  If  p”  a®  ^log(a  - 
z-0  (i 


cos^) (dp/d£)d£ 


(68) 


z‘^(m-1)|  p®  ^  a®  ^log(a  -  cos^)  (dp/d£)d£ 
0 


♦  z 


^  *f  p®‘^  a®‘‘'[a  -  cos^]"^(dp/d£)d£! 


If  ^  ^  0  along  the  entire  contour,  one  obtains  immediately 


'^m,1  2  *  ^  p'^logd  -  cos^)dp 


and  after  an  integration  by  parts 


Wjij  ^  2  ■  ^m+1)"'  p“*'(1  -  cos*)”' 3in^d(|» 


-1  .m+1 


-1 


for  n  -  1 


We  observed  above  that  w  ^  ■  0  for  n  >  3.  One  then  obtains 

m ,  n, ^  — 

from  Eq.  (66)  and  (67) 

-  (-)®*”B^(m+n)  ^  ^ p®*”cos°*3in”*(1  -  cos*)  ^ d* 

We  return  to  a  <■*-•-  v,  and  subsequently  to  the  original 
coordinates  (6-x)  and  (n-y).  One  has 


da  -  C(6-x)dn  -  (n-y)d5D/p‘ 


Therefore, 


mn 


(69) 


3^(m*n)"^  |(C-x)‘“(n-y)"[p  ♦  (5-x)]”%”U  (6-x)dn  -  (n-y)dt] 


-1  -1 


The  result  remains  the  same  if  at  some  point  of  the  contour  *  -  0 
There  n  “  y  ■  0,  and  6-x  <  0.  Then  one  writes 

(p  ♦  (6-x))~^  -  (p  -  (6-x)/(p^  -  (t-x)^) 

.  P  -  (6-x) 

(n-y)^ 


r'vS 


'‘v-_ 
-7fM 


m 

Hi 


-y.i 


1 


One  thus  obtains  an  integrand 


(C-x)‘“(n-y)""^(p  -  (t-x))p“^[(C-x)dn  -  (n-y)de] 


If  n  >  3»  then  there  is  no  singularity.  If  n  -  1 ,  there  is  one 
singular  point  where  the  integrand  behaves  as  (n-y)  The  value 
of  is  then  obtained  by  taking  the  principal  value. 

This  is  seen  from  Eq.  (67).  Let  £q  be  the  value  of  I  for 
which  ♦  -  0.  In  the  vicinity  of  this  point  we  consider  £  as 
variable  of  integration.  Let  £  -  £(♦).  Then  £(-ir)  -  0,  £(ir)  -  L, 
and  £(0)  -  £q.  In  the  first  integral  in  Eq.  (68),  but  with  limits 

£(-e)  L 

1  •  I  • 

0  £(*e) 

one  can  make  the  limiting  process  z+0  and  obtains 
£(-e)  L 

(-)"  I  +1  p“log(1  -  cos^) (dp/d£)d£ 

0  £(+e) 

In  the  region  -e  <  ♦  <  e,  we  take  ♦  as  variable  of  integration 


j  p^a®”^  (dp/d^)log(a'^  -  cos^ -  log(a  ♦ 


cos^)]d^ 


This  is  majorized  by 


I  const  log  {ar  -  cos‘^^)cos^d^  ♦  j  const  d<j) 


^  0^  2 

const  j  log(-^  sin  ♦)d(sin^)  +  const  2e 
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and,  by  an  integration  by  part 

*1  23ln  td(,at,nt)  ]  , 

^  +  sin'^^  J 

^  2 
P 

The  integral  is  0(e),  the  first  term  O(eloge).  As  expected  the 
term  vanishes  in  the  limit  e+0.  The  first  one  can  be  replaced  by 

L  ^ 

logd  -  cos(^} (dp/di)d£ 

0 

where  the  slash  indicates  that  an  e  neighborhood  of  t  -  *,q  is 

omitted.  Actually,  in  this  form  the  integral  is  well  defined  even 

if  one  does  not  omit  this  neighborhood.  But  if  one  carries  out  an 

integration  by  parts,  as  one  does  in  order  to  arrive  at  Eq.  (69), 

the  exclusion  of  such  a  neighborhood  is  necessary  and  leads  to  the 

definition  of  the  principal  part  of  the  integral.  The  second 

integral  in  Eq.  (68)  vanishes  in  the  limit  because  of  the  factor 
2 

z  .  In  the  third  integral  only  an  e  neighborhood  of  the  point  (or 
points)  where  ♦  -  0,  is  critical 

j  [p”  ^  a™  ^(a  ♦  cos<^)  ]^cos~%  +  sin^c^j”^  d(3in(j)) 

-E  ^ 

The  term  in  the  first  bracket  is  bounded.  Carrying  out  the 
integration  for  the  remaining  expression  one  obtains  as  bound 

+  e 

I  const  pz~^arctg(z  %sin(^)  j 


-  const 


’  ^  2 
sinolog  ^^2  *  sin‘^(>^ 


♦  E 


The  arctg  function  is  always  bounded.  For  z  -  0,  the  values  are 
-it/2  and  ♦ir/2.  The  whole  expression  is  then  0(z). 


For  n  even  the  contribution  of  the  sides  AB  and  DE  of  the 
region  of  integration  do  not  cancel.  Therefore,  the  vicinity  of 
p  -  0  must  be  discussed. 

As  before,  g  is  defined  by  Eq.  (56).  We  set  again 

g  -  ♦  gg  ^70) 

but  now  (n  even) 

g^Ca,*)  -  |[  cos^^sin*'^ 

g2(a,4») 

with 

*2  “  ~  co3<fr]  ^d(f  (72) 

Incidentally,  for  n  -  0,  m  -  0, 

g  -  gg  (73) 


-  a"(l  -  a^)”^^][a  -  cos^]  ^d^ 

m,,  2.n/2-  . 

-  a  (1  -  a  )  g2  (71  ) 


One  verifies  that 

gg  -  2(a^  -  1  )^ '^^arctg[(a  +  1)^^^(a  -  1)  ^  ^^tg(  (|)/2)  ] 
Using  the  definition  of  "a"  Eq.  (56)  one  finds 
(a  ♦  1)^^^(a  -  1)"^^^  -  (a  ♦  1)(a^  -  1)’^^  -  z~\p  *  ( 


Therefore, 


zg 


2 


2p  arctgCz  ^(p  +  (p^  ♦  z^)^  ^^)tg((>/2)] 


(74) 


The  expression  Eq.  (74)  will  be  encountered  only  if  p  0.  In  the 

2  2  1  /2 

limit  z+0  one  can  therefore  replace  (p  *  z  )  by  p.  One  then 
obtains 

I^Czig)  -  -  (z  small,  p  0  )  (75) 

z  +  4p  tg  4»/2 

Furthermore, 

lim  ^(zg^)  -  -  cot((J>/2) 
z-0 

Of  course,  the  expression  is  meaningless  for  ^  -  0.  If  p  0  and 

♦  »<  0,  the  last  equation  can  be  obtained  directly  from  Eq.  (72)  in 
the  following  manner. 

limO(zg.)/3z)  -  11m  gj,  -  f(l  -  cos^)~^d4>  -  [  (2sin^  ((|)/2)  ^  d4> 
z-^O  ^  z-»0  •'  ^ 


llm(3(zgp)/3z)  -  -cot(*/2) 
z*0 

For  ♦  -  ir  and  ^  -  -x  one  has  respectively  tg(^/2)  *>  •  and 
tg(^/2)  -  Therefore,  from  Eq.  (74) 


zg^{v)  -  px 

zg2(-x)  -  -px  (76) 

and 

3(zg2/3z  -  3(zg2/a2)  I  -  0 

4i-x  -X 

We  begin  with  the  contribution  of  g2  to  because  there 

the  vicinity  of  p  -  0  requires  some  extra  attention.  To  evaluate 
the  downwash,  one  has  to  form  Eq.  (59).  After  substitution  of 
Eq.  (71)  it  assumes  the  form 


-  2,  3  ,  ,  m+n-1  2.ti/2  -  ,  . 

“mn,2  “  ®  J  f*  ®  ^  g2(a,(f)dp) 


and  with  the  first  of  Eqs.  (58), 


2,  ,  vin+n+l +(n/2) ,  3.f  m-1  m-1  n,  - 


».n,2-  8  <-’> 


llm  ^(J  p“‘'a"‘'z"(zi2(a. ,))<!,)  (77) 


Along  the  portions  AB  and  DE  of  the  path  around  ABODE,  ^  •  -ir  and 
^  -  •*■•»,  respectively;  p(^)  varies,  respectively,  from  zero  to 
p(-v)  -  p(v)  and  from  p(v)  to  zero.  Substituting  Eqs.  (76),  one 
then  obtains  (for  the  two  portions  combined) 


f'*  p-a^-'op  .  (m-Dz"*' 


The  first  term  in  the  bracket  vanishes  for  all  even  values  of  n; 
for  n  -  0,  because  of  the  factor  n  -  0,  and  otherwise  because  of 
the  factor  z””^ . 

For  m  >  2,  the  Integral  in  the  second  term  gives  a  bounded 

quantity.  In  the  limit  z  -  0,  the  expression  vanishes  because  of 
the  power  of  z.  For  m  -  1  ,  the  second  term  vanishes  because  of 
the  factor  (m-1). 

For  m  -  0,  one  goes  back  to  the  first  formulation, 

Fq«  (77) •  Again  substituting  Eq.  (76),  one  obtains  for  the 
integral  along  AB 


n,  2  ^  2,1/2 

vz  ( p  +  z  ) 


Forming  lim  —  (...)  gives  0  for  n  >  2. 
z*Q 


For  m  ■  0,  n  -  0,  one  obtains 


(p2  .  z2)l/2 


pi  w) 


for  this  limit.  Accordingly,  the  only  contribution  of  for  the 
portions  AB  and  DE  of  the  path  is  encountered  for  m  •  0, 
n  -  0,  and  one  obtains 


-2itb‘ 


For  the  portion  BCD  (which  is  the  map  in  the  p4>  plane  of  the 
contour  of  the  5n-element)  the  expression,  Eq.  (77),  is  written  in 
the  form 

jm+n+1  *(0/2)  p  ^^2  a  )dp  +  ^  p  z  a  ■^(zg2)dp( 

11  ( in'*  1 

Here  p  /  0.  The  derivative  3(z"a''  )/3z  vanishes  in  the  lim  z+0 

for  all  m  and  n,  (da/3z  is  found  in  Eq.  (59)).  In  the  second 

integral,  lim  a™  ^ -  1.  If  the  integral  is  bounded  in  the  limit 
z*0 

for  n  -  0,  then  it  will  vanish  for  n  >  2,  (n  even)  because  of  the 
factor  z”.  Already  substituting  Eq.  (75),  one  obtains  for  n  -  0, 

B  (-)®+o'”^^  “2 - —  (79) 

6  (-)  T  P  2,2^  ^tg^^/2 

For  the  vicinity  of  f  -  0  the  integrand  behaves  as  ^  \  It  can 
be  shown  that  one  obtains  the  correct  lim  z+0  by  taking  in  the 
expression 


'(-)”  I 


p‘"'^cot(*/2)dp 


the  principal  part. 

The  use  of  the  principal  value  in  Eq.  (80)  is  justified  in 
the  following  manner.  First,  one  excludes  from  the  region  of 
integration  e-nelghborhoods  of  the  points  where  ^  -  0.  Outside  of 
these  neighborhoods,  one  can  make  the  limiting  process  z  -  0  and 
obtains  Eq.  (80).  It  remains  to  show  that  for  every  value  of  z, 
the  contributions  of  these  neighborhoods  tend  to  zero  as  e->0.  For 
this  purpose  we  set 


the  residue  is  a_^ .  Furthermore, 

log((|)2)  -  log|<j>2l 
log((>^  )  -  log|*^  I  +  iir 
Denoting  by  P  the  principal  value  one  finds 


P  J  a^  +  a^4i)d<(.  «  Clog) 'l>2 1  '  log^i  I  ~  *^1^ 

♦  a,  (<>2  ' 

In  every  situation  of  this  kind  the  integral  will  contain  a 
logarithm.  The  contribution  of  the  residue  is  automatically  taken 
into  account  by  taking  the  absolute  values  of  logarithms. 

This  makes  an  additional  step  possible.  In  the  integral  in 
the  complex  plane 


*  *2  ^2 

n  o  -  P"*"^  cot(^/2)(dp/d(t.)d^ 

m,  u ,  2 

an  integration  by  parts  is  carried  out.  Then 


^  '2  r 

r  * '' 

-  8^(-)"fm' 

p”  cot((^/2) 

L 

'm,0,2  I  - 

^2  m  1 

♦  f  T ^ - r  residuei 

J  1  -  cos  4* 


Taking  the  above  example  and  carrying  an  integration  by 
parts  one  obtains 


••1 

j  (a_^  ♦  ♦  aQ  + 


h  *2 


.  <^(a_^  ^  +  a^4>) 


"1  '*’1 


The  integrand  after  the  integration  by  parts  has  the  same  residue 
as  the  original  integral.  Again,  one  takes  in  the  analytic 
expression  for  the  resulting  integral  the  absolute  value  of  the 
logarithm. 

If  one  carries  out  the  integration  only  along  a  portion  of 
the  contour,  then  one  must  of  course  include  the  expression 
outside  of  the  Integral.  If  one  makes  this  integration  by  parts 
along  all  sections  of  the  contour,  then  the  terms  outside  of  the 
integral  cancel  each  other. 

The  same  integration  by  parts  can,  of  course,  be  made  for  a 
contour  of  the  5n~element  for  which  ^  /  0  everywhere.  Thus  one 
can  always  write 

Wm  Q  2  "  <1*  “  coS((»)  ^d(fr  (n-  0) ,  m  >  1  (81  ) 

if  one  uses  for  the  individual  sections  of  the  contour,  analytic 
expressions  for  the  integrals,  substitutes  the  limits,  and  takes 
the  absolute  values  of  the  logarithms  which  may  occur. 

For  m  -  0,  n  -  0,  the  evaluation  of  the  expression  Eq.  (80), 
now  given  by 

6^  ^  p‘^  cot((>/2)dp 
is  carried  out  as  follows.  One  has 

cot  i<^/2)  -  (1  +  co3^)/sin0  -  -(1  -  co3a)/sincx 


•  J‘  •  _•  »  ^  .*  * 


We  defined  C-x  -  p  cosa,  n-y 


p  sina.  Accordingly, 


a  p’‘'cot(<D/2)dp  -  -  +  L  U-x.)P.dp. 

T  ^  y  1  p'^(Ti-y) 


The  second  integrand  is  transformed  into 


(5-x)[ (5-x)d5  +  (n-y)dn]/(p  (n-y)) 


[(5-x)^d5  +  (?;-x)  (n-y)dn]/(  P^(n-y) ) 


2  2  2 
Here  (5~x)  is  replaced  by  p  -  (n-y)  . 

One  then  obtains 


(dC/(n-y))  +  (C(5-x)dn  -  (n-y)d5]/p‘^)  =  (d^/n-y))  +  dcj) 


For  points  xy  outside  or  inside  the  5n-element,  one  has, 
respectively,  d^  -  0  and  ^  d$  =  2ii.  The  second  expression 
cancels  the  contribution  to  Wq  2  shown  in  Eq.  (78).  Moreover  - 
0  for  m  ■  0,  n  -  0.  Thus  one  obtains,  generally. 


w 


00 


2  r,  (d^/dt)  -  (dp/dO 

T  h-y 


(02) 


For  (n-y)  -  0, 


dp  -  dCsign(C-x) 

If  n-y  -  0  and  (5-x)  <  0  (x  in  the  wake  of  £),  a  singularity  will 
arise  and  the  expression  must  be  interpreted  (according  to  the 
above  discussion)  as  principal  part.  In  the  discussion  following 
Eq.  (78)  we  found  that 

”m, n, 2 


0 


for  even  n  >  2 


(83) 


The  portions  AB  and  DE  of  the  contour  in  the  p,^  plane  start  at 
p  -  0.  There,  ^  -  -ir  and  +ir,  respectively,  and  one  has  to  form 
according  to  Eq.  (59) 


tI  (a,-ir)  -  g^(a,v)]dp| 

.  lim  I  [  p®*”"^  a"^  [g- (a,-ir)  -  g,(a,if)]dp 

9!-»n  W 


p®*n"l  a'l  [g^  (a,-ir)  -  g^(a,if)]dp 


2  j  |_(a“^  [g^  (a,-ir)  -  g^  (a,iT)])dp| 

A  / 


Forming  the  limit  z+0  one  replaces  in  the  first  integral,  a  by  1 

2 

The  second  Integral  vanishes  because  of  the  factor  z  if  m+n  >  3 
We  are  discussing  cases  for  m+n  >  1  and  n  even.  Since  n  is  even 
m+n-1  implies  m  -  1.  Then  by  its  definition,  Eq.  (71) 


g^  -  I  (cos^  -  a)(a  -  cos^)~^d^ 


and  one  has  to  form 


p(-v) 


a-' op) 


lim  ^(P^  ♦  z^)^^^ 
z+0 


p(-ir) 


-  11m  I (p^  +  + 


2 

T77^]  I 


p(-v) 


p(-n) 


This  is  the  same  as  if  one  had  set  a  -  1  in  the  j  a  ^  dp.  For 
m+n  -  2,  one  has  either  n-2,m-0,  orm-0,  n-2.  In  each 
case  one  deals  (except  for  the  sign)  with  an  integral 

2  2 

g^(a,<>)  -  I  ~c|g|~r  'I"  <!♦  *  I  (cos4>  +  a)d(t»  -  ain«()  +  a4> 


Then  g^(a,it)  -  g^(a,-ir)  -  2na  and  one  has  to  form 

lim  -j-  (z  f  p2irdp)  -  iip^(n) 

z-^0  J 

One  thus  obtains  in  all  cases  as  combined  contribution  of  the 
portions  AB  and  DE  due  to  g^ 

B  (-1)  (g^df-ir)  -  (1 ,  ♦»)  (m+n)  p  (if) 

This  expression  cancels  the  term  outside  the  integral  in  Eq.  (85). 
One  thus  obtains  (after  returning  to  a) 

Wjjj  ^  ^  J  p™  *^coa***a  sin'™(o)[1  coso]  ^do,  (86) 

even  n  >  2 

Since,  according  to  Eq.  (83)  2  ”  ®  h  >  2,  this  is 

(rather  than  ,).  Moreover 
mn ,  1 

'^m  0  1  "  ®  ®  1  ^  (cos  o  -  (-1)  )[1  coso]  do  »  m  >  i 

This  is  combined  with  the  expression  *  o*  (31).  One  then 

ID  f  U  f  c 

obtains 

^m,0  ”  ”m,0,1  *  ”m,0,2  “  ^  P^cosVl  ♦  cosoD'^do,  (86a) 

if,  as  was  stated  above,  one  uses  for  the  individual  sections  of 
the  contour  analytic  expressions  for  the  integral,  substitutes  the 
limits  and  takes  the  absolute  value  of  the  logorithms,  if  they 
should  occur.  Actually,  with  this  interpretation  Eq.  (86)  is 
generally  valid. 


These  results  give  the  contribution  to  K  obtained  for  the 
power  zero  of  the  frequency  k,  which  will  be  indicated  by  a 
superscript.  In  the  original  coordinates  one  obtains 

w°n(xy)  “  B^(in*n)“^  (5“x)“*(n-y)''[p(p  +  (£-x) )]  ^  [  (5-x)dn-(n-y)d5] 

(87) 

except  for  m  -  0,  n  -  0 


For  n  odd,  this  is  a  repetition  of  Eq.  (69).  For  n 
follow  the  procedure  described  after  Eq.  (86a). 

Furthermore, 


0 ,  one  must 


WQo(x»y)  "  <|ii(dC  '  dp)(n-y; 


Eq.  (51)  gives  the  result  for  the  first  power  of  the  reduced 
frequency.  Introducing  again  the  form  Eq.  (52)  for  Ap  considering 
a  single  term  and  carrying  out  the  integration  with  respect  to  p, 


one  obtains 

wL  -  (m+n+l)'"'  [(C-x)“(n-y)®"k(C-x)dn-(n-y)d5]  (88a) 

mn  j 

The  results  of  these  discussions  can  be  summarized  in  a 
simple  manner.  We  assume  that  for  x,y  fixed,  Ap  has  the  form  of 
Eq.  (51),  l.e.,  a  development  in  powers  of  (t”X)  and  (n-y). 

Except  for  m  «  0,  n  -  0,  the  correct  expression  for  the  upwash  at 
z  -  0  is  obtained  by  setting  z  •>  0  in  the  expression  for  K.  One 
then  has  to  evaluate  integrals 


I  -  II  F(5-x, n-y)cicdn 


For  z  >  0,  one  has  R  -  p.  Let 


5-x  -  R  cosa 
n-y  -  R  sino 


For  one  of  the  summands  in  Ap  and  one  power  of  the  reduced 
frequency  k,  the  function  F  has  the  form 

F  -  R'^fCo)  (90) 


Then 


I  -  JJ  R‘*f(a)RdRdo 

-  (q+2)  ^  R{  o)^*^f  (o)da 

or  in  terms  of  the  original  coordinates 

I  -  (q*2)‘^  ^  F(c,n)C(e-x)dn-{n-y)d5]  (91) 

If  q  -  -2  (i.e,  for  m  -  0,  n  -  0  and  for  the  power  zero  of  the 

reduced  frequency),  the  procedure  is  not  feasible.  The 

integration  with  respect  to  R  gives  log  R.  The  limits  for  R  are  0 

and  R(o).  The  lower  limit  gives  infinity,  which  is  an  indication 

that  this  simple  minded  procedure  fails.  Here  it  is  necessary  to 

proceed  in  the  manner  described  above;  one  carries  out  limiting 

process  z-»0  only  after  the  integrations  with  respect  to  ^  and  n 

(or  p  and  a)  have  been  carried  out.  This  gives  Eq.  (88).  Another 

limitation  of  this  rule  arises  if  n  •  0,  m  >  1,  the  power  of  the 

reduced  frequency  is  zero  and  if  along  the  contour  n-y  becomes 

2 

zero.  Then  the  Integrand  will  contain  a  denominator  (n~y)  ,  and 
the  integral  becomes  nonsensical.  One  then  has  to  apply  the 
procedure  described  after  Eq.  (86a). 

In  ,  one  finds  terms  which,  after  one  sets  z  «  0,  gives 
rise  to  Integrals 


I 


jj  log  (n-y)F(6-fTi-y)d5dn 


(92) 


where  F  has  the  form  assumed  above.  Then 


I  -  JJ  log(Rsino)R‘^*^  f(a)dRda 

-  jj  (log  R  +  log  sina)R‘**V(a)dRda 

-  (q*2)"^  j  f(o)R‘^*^[logR(o)  -  (q*2)"^  +  [log 
This  is  written  in  terms  of  the  original  coordinates 

I  -  (q-*-2)~^  ^  F(t,n)[log(n-y)  -  (q+2)  ^][(5-x)dn  -  (n- 


sina) ]da 


y)dt]  (93) 


53 


SECTION  V 

THE  UPWASH  AT  A  GIVEN  POINT  (x.y) 


We  evaluate  In  this  section  the  upwash  at  a  fixed  point  x,y 
due  to  a  pressure  distribution  in  a  given  ^n'element.  To  each 
triangular  or  trapezoidal  elements  belong,  respectively,  3  or  4 
elemental  pressure  distributions.  The  results  can  be  expressed  in 
terms  of  elementary  transcendental  functions,  provided  that  the 
elemental  pressure  distributions  have  the  form  of  polynomials  in  5 
and  n.  These  formulae  will  be  needed  if  (x,y)  lies  within  or  in 
the  vicinity  of  a  Cn~element.  For  triangular  element  and  Ap 
linear  in  5  and  n,  one  obtains  pressure  distributions  which  are 
continuous  at  the  common  boundary  of  neighboring  elements.  For 
trapezoidal  elements,  one  can  achieve  only  imperfect  matching  even 
if  one  uses  terms  of  higher  order  in  t  and  n.  Except  for  this  flaw 
in  the  basic  data,  the  analysis  is  Independent  of  the  shape  of  the 
elements. 


For  triangular  elements  the  corner  opposite  the  side  that  is 
parallel  to  the  C  axis  has  the  number  1,  from  thereon  one  proceeds 
in  the  counterclockwise  direction.  The  elemental  pressure 
distributions  are  then  given  by 


where 


(1) 

(i) 

(i) 

(i) 

Ap 

U, 

n)  - 

♦  (5 

)  ♦  ©2  (n-n^ ) 

(1) 

-  1, 

(1) 

(2) 

°0 

-  0, 

■ 

-C  n2“n^ ) 

(2) 

(2) 

-1 

(2)  -1 

-  0, 

°1 

-  D 

(n^-n^ ) 

,  C2  -  -D  ( ) 

(3) 

(3) 

- 

1 

(3)  -1 

-  0, 

.  -D 

(n2“ni 

) .  ^2  “  D  ^ 

D 

-  [( 

,)(n3 

-ni)  - 

)(  ^2“ H'j  )  3 

i  -  1.2,3 


(94) 
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U.11. 


The  upwash  is  evaluated  separately  for  these  three 

.  *-(i) 

expressions  Ap 

We  shall  write  during  the  derivations  of  the  formulae 


Ap(5.n)  -  Zcju^{x,y)(c-x)'“(n-y)” 

For  triangular  elements  one  has  in  particular 

(i)  (i)  (i) 

Ap(5tn)  -  Coo(x.y)  ♦ 

where 

(!)  (i)  (i)  (i) 

Cqo  (x.y)  -  ^0  "  °i 

(i)  (!)  (i)  (!) 

°10  *  **  °01  ‘  °2 

Later,  in  integrations  over  x  and  y,  each  of  the  expressions 

(i)  (i)  (i)  (i) 

Ap  (x,y)  -  Cqq  (x,y)  ♦  c,q  (C-x)  ♦  Cq^  (n-y) 

will  be  used  in  three  different  forms. 

(i)  (i)  (i) 

Ap  (x,y)  -  c^Q  *  °01 

Here  the  subscript  i  refers  to  those  corners  of  the  £n“ 
element  for  which  the  elemental  pressure  is  1  .  (For  the  other 
corners  the  elemental  pressure  is  zero);  j  is  the  subscript  of  any 
of  the  three  corners  when  it  is  used  as  origin  of  a  x,y-3ystem. 

The  discussion  of  the  limiting  process  z+0  in  Section  V  has 
led  to  integrals  around  the  contour  of  the  ^^"element  which  give 

the  upwash  at  a  fixed  point  x,y  due  to  a  pressure  distribution 

If  the  point  xy  lies  in  the  interior  of  the  5n-element 
or  of  its  neighbors  then  the  integrands  of  these  expressions  are 


i 

\4 


•  •  A.*  m  ^  m  i*'*'**#  ’**•  "*•  ^ ^ •’  .  •*-  “*•  *  “  *  •  *'\  •  *  *  *  "  •  \  V* 


analytic  functions  except  for  well  defined  singularities.  The 
author  hoped  originally  that  the  smooth  part  of  these  integrals 
could  be  evaluated  numerically.  But  upon  closer  examination!  he 
realized  that  for  a  point  (xy)  close  to  the  boundary  of  the 
element  the  integration  over  C  and  n  will  encounter  difficulties  in 
spite  of  the  analyticity  of  the  integrands.  Fortunately,  the 
integrals  can  be  evaluated  in  terms  of  elementary  transcendental 
functions,  although  one  obtains  fairly  lengthy  expressions. 

The  difficulties  arise  in  the  following  manner.  The 
integrands  are  analytic  for  real  values  of  the  variable  of 
integration,  but  in  some  vicinity  of  the  origin  they  have 
singularities  (poles  and  branch  points)  at  complex  values.  As  the 
point  (xy)  moves  closer  to  a  boundary  of  the  5n  element  the 
interval  of  Integration  extends  farther  and  farther  along  the  real 
axis  (if  one  keeps  the  singularities  in  fixed  positions).  One 
might  then  divide  the  region  of  integration  into  a  section  close  to 
the  singularities  and  a  remaining  part.  But  even  the  second 
section  is  not  well  suited  for  numerical  integrations.  Let  z  be 
the  variable  of  integration  and  f(z)  the  integrand.  A  suitable 
variable  of  integration  in  the  outer  region  would  then  be  1/z  and 
one  must  form 

-  I  z^  f(z)d(l/z) 

where  (1/z)  is  the  variable  of  integration.  This  works  well  only 
if  z^f(z)  is  a  regular  function  of  (1/z)  for  1/z  -  0.  The  terms 
which  violate  this  requirement  are  fairly  easily  identified,  so 
that  one  could  treat  them  separately.  But  since  an  analytical 
integration  over  the  whole  range  is  possible  it  appears  more 
practical  to  follow  this  course.  Naturally,  once  programmed  one 
will  use  these  formulae  also  if  the  point  (x,y)  is  not  close  to  the 


16 
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element  boundaries.  Later  this  procedure  will  make  It  possible  to 
carry  out  also  the  integration  over  x  and  y  analytically.  The 
integrals  to  be  evaluated  are  given  by 


IJ  Ap(5-n)K(5-x,n-y,k)d5dTi 

The  function  K  is  found  in  Eq.  (B.63) 

K  -  B~^exp(ik(5-x)[K^ (5-x,r,k)  +  K2(5-x.r,k)] 

where,  according  to  Section  IV,  r  can  be  replaced  in  most  cases  by 

|n-y|. 

We  have  introduced 

Ap(5.n)  -  exp(ikt)Ap(5,n)  (97) 

Within  the  elements  Ap(Cih)  is  represented  by  shape  functions  which 
are  linear  in  5  and  r\.  The  infinite  sums  in  Kg  can  be  truncated 
after  only  very  few  terms,  because  k  is  usually  small  and  because 
the  present  discussion  refers  to  points  (x,y)  within  or  close  to 
the  Sn-element.  In  the  following  term  of  Kg  the  exponential 
function  is  expressed  as  a  power  series. 


(RV)"^  (exp(ikV)  -  1)  -  ikR'^  E  (ikV)^“V)ll 

1 


2  2-1/2 

In  Ueda's  formulation  one  finds  a  similar  term  namely  R(V  ♦  r  ) 
exp(ikV).  The  present  simpler  form  arises  from  a  combination  of 
this  term  with  other  expressions  of  K,  (see  the  discussion 
following  Eq.  (B.59).  In  the  present  context  this  brings  about  a 
considerable  simplification  in  the  integrations.  In  the  procedure 
of  Ueda  and  Dowell  the  simplification  is  only  minor.  All  terms  of 


Kg  can  now  be  brought  into  a  form 

P(5-x,Ti-y)  ♦  r"^  P(^-x,n-y) 


where  P(e-x,n-y)  denotes  a  generic  expression  for  a  power  aeries  in 
their  arguments.  In  practice  one  deals  with  polynomials. 

We  shall  derive  integration  formulae  for  the  lowest  powers  of 
the  reduced  frequency  k  of  the  development  of  K,  namely  k°,  k, 
k^logk  and  k^ .  For  a  portion  of  Ap  given  by  ( 6-x)™(  iry .  we 

-  Jt 

shall  denote  by  I  the  integral  over  the  5n-element.  The 
superscript  I  refers  to  the  power  of  the  reduced  frequency  k. 
does  not  contain  the  factor  (8it)~^P  ^exp(-ikx).  After  substituting 
Eqs.  (B.66)  and  (B.67)  one  obtains  the  following  expression 

-  3^  j|(5-x)‘“(rj-y")r"^R"^(R  -  (C-x))dedn;  m+n  >  1  (98) 

(The  case  m  -  0,  n  -  0  has  been  treated  separately.) 

-  “ik  ||(t-x)”(Ti-y)”R  ^d^dn;  m+n  >  0  (99) 

2 

ijn  “  7  Jf(5-x)®(Ti-y)"CC  -  (b"^(5-x)r"‘' ) 

(100) 

♦  (log(k(R  +  (C-x)))]dCdn 
Here, 


-  (n-y)^,  R  -  [(C-x)^  *  (n-y)^]^^^  (101) 

The  function  V  which  occurs  in  Eq.  (B.67)  is  found  in  one  of  the 
Eqs.  (B.61 ) . 

C  -  -(1/2)[4»(1)  *(2)  >  log(l-M)  +  21og2  -  iir  +  3“^M] 


(102) 


Then,  according  to  Eqs.  (89),  (91),  (92),  and  (33) 


i°n  -  8^(ra+n)‘‘'j^  (C-x)”(n-y)"'2R'''(R-(5-x))[(5-x)dn-(n-y)d5]  (103) 
-  -ik(m+n+1)  ^  ^  (5-x)'”(n-y)'’R  ^  [  (  5-x)dn- ( n-y)d5]  (104) 

“  (K^/S)  (m+n+2)“''  ^  (5-x)‘"(n-y)'’lc-(6’^(5-x)R~‘' ) 

(105) 

+  (log(k(R  +  (C-x))  -  (m+n+2)  M[(5-x)dn  -  (Ti-y)dC] 


The  path  of  integration  is  the  contour  of  the  element.  Practically 
this  means  that  we  have  separate  expressions  for  the  sides  of  the 
triangular  or  trapezoidal  element.  Accordingly  we  write 


12jJl  ^  23jt  ^  31jt 
mn  mn  mn 


(106) 


where  the  left  superscripts  refer  to  the  end  points  of  the  sides 
(as  one  travels  around  the  ^n-element  in  the  counterclock-wise 
direction).  We  shall  derive  separate  formulae  for  the  sides 
parallel  to  the  ^  axis  and  sides  inclined  to  the  ^  axis.  For  a 
aide  of  the  ^p-element  between  the  point  i  and  i+1,  not  parallel  to 
the  ^-axia,  (where  for  triangles  i+1  -  4  is  understood  to  refer  to 
point  1 ) ,  we  define 


■  ‘'l.l  ■  -  '>!> 

The  subscripts  of  a  will  frequently  be  omitted.  The  expressions  I 
appear  first  as  indefinite  Integrals.  After  substitution  of  the 
limits  of  integration  (expressed  by  and  one 

obtains  two  expressions  (x-^^^,  y-yj^)  and  (x-^^  +  i  tyn^  +  i )  •  Those 
expressions  will  be  denoted  by 


so  that 


mn  "  "^inn 


’“l.l*! 


<>1,1.,)  008) 


Information  about  the  side  for  which  the  expression  is  being 
evaluated  enters  through  the  argument  . 

For  the  side  parallel  to  the  5  axis,  we  insert  n/2  as 

argument  of  a.  The  expressions  I(..o)(to  be  found  later)  tend  to 

infinity  as  a  +  ir/2,  but  I(..ir/2)  remains  finite.  This  happens 

because  of  the  occurrence  of  a-dependent  integration  constant  which 

tends  to  infinity  as  a-»iT/2.  They  cancel  as  one  forms  the  integral 

2  2  1/2 

between  points  i  and  i+1 .  The  expressions  R  -  ((^-x)  +  (n-y)  ) 

will  introduce  branch  points  of  the  second  order  in  the  integrands. 
We  write 


I 


I 

m.n 


♦ 


(109) 


The  additional  superscripts  1  and  2  refer  respectively  to  the 
absence  and  presence  of  such  branch  points.  For  a  -  ir/2,  one  has 
n  -  hi  ■  ■  const,  dn  ■  0.  Then  we  shall  introduce, 


-  y  -  V 


(110) 


(^-x)  -  w.v 


(111) 


and  consider  w  as  variable  of  integration.  For  the  power  zero  of 
the  reduced  frequency,  one  has  according  to  Eq.  (90) 


^00  ■  ‘8^1  ^  /  %] 

R  is  defined  in  Eq.  (101).  For  -  ir/2  one  has  to  evaluate 
(because  dn  *  0) . 

i,i  +  1  ^0  i,i  +  1 .01  .  0,i  +  1  .02 
■^00  “  ■‘■00  ■‘■00 
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^00  ■  "  ®  J  (n,  -  y) 

r  * 


Hence , 


i,i+1  02  2  [  (£-x)d£ 

00  ■  J  RTn/  -  yl 


1,1  +  2,02  .  A2  £  \ 

^00  -  ni  -  y) 


and  with  Eq.  (Ill), 


■  “0^  slgn(nj  -  y)  j  w(w^  +  1)“’^^dw 


2  2  1/2 
-6  slgnCrij^  -  y)  (w  +  1  )  ' 


If  one  goee  back  to  the  original  variables,  then  one  obtains 

(w^  +  -  sign(nj  -  y)R/(nj  -  y) 

because  R  Is  always  positive.  The  sign  functions  cancel  in  the 
original  coordinates.  This  happens  in  all  formulae.  For 


simplicity  we  shall  omit  in  future  computations  such  intermediate 
sign  expressions.  For  future  work  it  is  convenient  to  introduce 


X  -  q,  y  -  y  - 


U  -  -u,  V  -  -V 


(113) 


Then, 


-  e^X/Y 


(114) 


lQQ(X,Y.Tr/2)  -  e^R/Y 


(115) 


According  to  Eq.  (103)  one  has  to  evaluate  for  a 
zero  of  the  reduced  frequency  k,  and  m+n  >  1 


i,i+1j0  _  i.i+liOl  ^  i,i+1=02 

■*’mn  "  ■'■mn  "^mn 


with 


“  -B^(®^n)'^(ni  '  y)""''  }^\^-x)‘"d5 


■6^(m+n)  ^  (m+1 )  ^(n^  ^  y)"^  ^(C“x)”*^ 


ir/2,  the  power 


i,i>l70,2  „2,  ,-1,  ,n-1  f  ,  ,m+1„-l 

-  6  (m+n)  (n.  -  y)  I  (e-x)  R  d£ 


with 


w  -  (i-x)/(x\^  -  y) 


one  obtains 


1,1+1 jO, 2 
■‘‘mn 


P^(m+n)~^ (n. 


{  \  2  ,  v-l  /2  m+K 

I  (w  +  1 )  w  dw 


Formulae  for  the  integrals  occurring  here  are  found  in  Appendix 
C.  Specializing  to  [m.n]  -  [1,0]  and  [0,1]  one  obtains 

I°Q  (X,Y,ir/2)  -  (g^/2)X^/Y 
0?  ? 

IJ'q  (X,Y,ir/2)  -  (B  /2)[(RX/Y)  +  Ylog(k(R  -  X))]  (116) 

I^]  (X,Y,ir/2)  -  B^X 
Iq^  (X,Y,Tr/2)  -  B^R 


The  factor  k  in  the  logarithm  amounts  to  a  change  of  the  constant 
of  integration.  If  the  expression  is  written  in  this  form,  it  will 
not  change  if  one  changes  the  reference  length  L. 

For  o  -  ir/2  and  the  power  1  of  the  reduced  frequency  k  one 
obtains  from  Eq.  (104) 


-  *  ik(m+n+1)‘^  J  r"'' ( e-x)”"( -  y)'™^^dC 
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l.i*lT23  . 
mn 


■>(k^/2){m+n+2) 


(5-x)'”(nj  -  [log(k(R  ♦  5-x))-(in+n+2) 


Hence 

IqJ  (X,Y,ir/2)  -  -(k^/4)CXy 

(X.Y,ir/2)  -  (k^/12)CX^y  (118) 

Iq]  (X,y,ir/2)  -  +  (k^/6)cxy^ 

22 

In  one  Introduces  w  -  {5“x)/(ti^  -  y)  In  an  Intermediate  step 
and  thus  obtains  integrals  listed  in  Appendix  C.  Hence 


Iqq  (X,y,»/2)  -  -(k^/i»)6‘^yR 

I^Q  (X,y,x/2)  -  (k^/12)B”^[yXR  ♦  y^log(k(R  -  X)]  (119) 

Iq^  (X,y,ir/2)  -  (k^/6)B'^y^R 

Furthermore, 


i.i+lT2,3 

mn 


-2(k/2)^(m+n+2)“'' (hj  -  y)"*^  I 


with 


(5-x)‘“[log(k(R  ♦  5-x)) 


(m+n+2)~^  ]d5 


1 


I  -  (m+1)  {(e-x)  Clog(k(R  +  5-x))  -  (m+n+2)“  ] 


-  /*'  <5-x)"*'R''dU 

With  w  -  (5~x)/(nj  ~  y)  one  obtains  for  the  last  integral 
expressions  listed  in  Appendix  C;  in  particular 

J  (e-x)R"^dC  -  R 

J  (5-x)^R'’d5  -  (l/2){(c-x)R-(n^  -  y)^log(i<(R  ♦  ^-x))} 
Therefore 

^00  “  -(k/2)^(Tij  -  y){(C-x)[log{k(R  ♦  e-x))-l/2]  -  R} 


Iq^  (X,Y,ir/2)  -  -(k/2)^{yXClog(k(R  -  X))  -  (1/2)]  +  YR) 
{X,Y,it/2)  -  (1/3)(k/2)^{YX^ClOg(k(R  -  X))-1/3] 

+  (1/2)  XYR  +  (1/2)Y^log(k(R  -  X)) 

(X,Y,if/2)  -  (2/3)(k/2)^{Y^X[log(k(R  -  X))-1/3]  ♦  Y^R} 


This  terminates  the  case  a  -  it/2. 


For  a  ^  ir/2  auxiliary  systems  of  coordinates  are  used.  The 
"sheared*'  coordinate  system  is  defined  by 

u_  -  (5-x)  -  (n-y)tga,  i-x  -  u_  ♦  tga  (121) 

3  3  9 

Vg  -  n-y.  n-y  -  Vg 

Along  the  side  iti^l  of  an  element  one  has 

Ug  -  Ug^  -  -  const.  -  -  x)  -  (n^  -  y)tga  (122) 

Note  that 

(e-x)dn  -  (n“y)d5  -  u„dv„  -  v_dua  -  u_dv_  (123) 

S  3  9  9  9 

(The  last  simplification  occurs  because  u_-  const  because  along  the 

3 

side  i,i*1,  and  therefore  du_  -  0.) 

9 

The  "rotated"  system  of  coordinates  is  defined  by 
u  -  (C-x)cosa  -  (n-y)sina,  5-x  -  u  cosa  ♦  v  sina 

(124) 

V  -  (5-x)sino  ♦  (n*y)cosa,  n-y  ■  -u  sino  +  v  cosa 

Along  the  side  i,i+1  of  an  element  one  has 

u  -  Ui  -  Uj^^  -  const  (125) 

u„ .  -  cos  ^  a  u.  (126) 

SI  i 


Note  that 


Therefore, 


(C-x)<16  ♦  (Ti*y)dn  -  u  du  ♦  V  dv  -  V  dv  i 

The  last  simplification  occurs  because  u  -  const  along  a  side 
i,i+1.  Moreover 

(5-x)dn  -  (n~y)d5  ■  u  dv  -  v  du  -  u  dv  i 

We  shall  furthermore  introduce 

w  -  v/u,  I 


(128) 


(129) 


0  -  arctgw  -  arctg(v/u^) 


(130) 


(131 ) 


Then 


sine  -  v/R 


(132) 


cose  -  Uj/R 


Hence , 


sign(cos0)  -  sign  Uj 


(133) 


Moreover, 


sin(0-a)  ■  (v  cosa  -  UjSino)/R  -  (Ti~y)/R 


cos(0-a)  -  (u^coso  +  V  sina)/R  -  (5-'x)/R 


(13i») 


The  basic  formula  for  o  ^  ir/2  and  the  power  zero  of  the 
reduced  frequency  is  quoted  in  Eq.  (112).  We  introduce  the  rotated 


.V •  • 


system  of  coordinates  Eqs.  (124)  through  (129)  and  Eq.  (130). 
Writing 


i,i  +  1=0  i,i  +  1::01  i,i  +  1=02 

00  “  -^00  00 


one  obtains 


w 


i  +  1 


iji+I^OI  o2  f  4  ,  .  v-1 

Iqq  -  6  I  sino(w  COSO  -  sino) 


dw 


1,1+1=02 


1  +  1 


00 


sign  Uj  I  (wcoso  -  sina)~Vw^  +  1)”^'^^wdw 


w. 


Eq.  (135)  gives 


w 


log(w-tga) 


1  +  1 


w. 


Introducing  9,  Eq.  (I31),  one  obtains 


"  S^tgo  log  (3in( 9-a)/cos0 


1  +  1 


In  Eq.  (136),  0  is  introduced  immediately.  One  observes  that 

(w^  +  l)^'^^  -  sign(co30)co3  % 


(135) 


(136) 


(137) 


ti; 


lyy  V  "jj  ■■'-*  v*jf"ji  ■- 1-'  *’^\’^i’*v»i’f  i’^  I  ^'.’^  '■’^v^  tf*  'jv  i™ 


I 

\'k 


s 


Then  with  Eq.  (133) 

1,1*1x02 

^00 


1*1 

oosLfnCe-a)  ■  I  ‘^'"8“  3in(e-a)''  *  C03''a  cos“'0]cie 


’1*1 


-1 


-1 


“0^{tgo  log[3in(0-o) (1 *co3(e-a))  ^]  ♦  cos  ^ olog[( 1 +3in0)co3  ^0]} 


Combining  the  expreselone  Eqe.  (137)  and  (138)  one  obtalne 

“  B^{tgalog((1 +co3(0-o))co3  ^0) 

6i 


C' 

(138) 

•  *  < 

0]d0 

a 

C] 

(139) 


1  +  1 


-  C03  ^ alog(  1 +sin0)co3  %)) 


We  return  to  the  original  coordinates,  but  partially  retain  u^ . 
Again  the  definition  Eq.  (113)  are  Introduced.  In  addition  we  set 


V  -  -  V  -  Xslna  +  Ycoso  (140) 

U  -  -  u  -  Xcosa  -  Ysino 

Then  (with  a  different  constant  of  integration) 

Iqq(X,Y,o)  -  0^{tga  log(k(R-X))  -  cos"’a  log((R  -  V)/|Ul)} 


R  -  (X^  + 


(141  ) 


The  basic  formulae  for  the  power  of  the  reduced  frequency  0 
and  m+n  >  1  l3  Eq.  (103).  We  write 


I®  . 

mn  mn  mn 


70 


V, 

s.' 


m 


m 


ill  -  e^(m+n)‘’  J(5-x)®(n-y)""^[{C-x)dTi  -  (n-y)de]  (l‘»2) 

i°n  -  -B^(m*n)‘’  J(5-x)“*^n-y)"‘^R‘‘'[(C-x)dTi  -  (n-y)dc]  (143) 

For  n  -  1,  the  integrand  contains  a  factor  (ri“y)  The 
integrals  are  then  interpreted  as  principal  values. 

-2 

For  n  -  0,  one  has  a  factor  (n-y)  .  In  Eq.  (103),  (from 
which  Eqs.  (142)  and  (143)  arise)  no  singularity  is  encountered  if 
5-x  >  0.  The  singularities  in  Eqs.  (142)  and  (143)  appear  because 
of  the  separation  into  terms  with  and  without  R.  The  separation  i 
necessary  for  the  derivation  of  analytical  formulae.  In  the 
evaluation  of  the  integrals  one  follows  the  procedure  described 
after  Eq.  (86a) . 

Introducing  sheared  coordinates,  Eqs.  (121)  through  (123), 
one  obtains 


i,i+1 jOI 
■‘'mn 


B^(m+n)“^ 


i+1 


(u 


si 


n-2 


Developing  (u^^^  +  v^tgo)”  one  obtains  elementary  integrals. 
Specializing  immediately  one  obtains 


i,i+1 ^01 

1,0 


I't-(u32/v^) 


tga  Ug^lOg|kVg|] 


3i  +  1 


si 


(144) 


i,i+lT01 

01 


S  UgjlogjkVgl 


si  +  1 


si 


These  expressions  will  be  rewritten  in  terms  of  Y  and  U  in 
Eqs.  (147). 


The  evaluations  of  ^  ^mn  rotated  system  of 

coordinates  Eq.  (124)  through  (129). 

Y 

^'^*^^mn  "  ~B^(ra'*’n)  ^  |  (u^cosa  ♦  vsina)*”*^ 

(-Uj^sina  +  V  coso)"^  ^^*^1^  *  ^^^Uj^dv 

and  with 


w  -  v/u^ 

‘'i  +  1 

B^(m+n)  'u°*  I  sina  +  coso)*** 

(w  cos  -  slna)”  ^(w^  ■*•  1)  ^  ^^dw 

For  n  >  2  one  writes  the  product  of  the  first  two  terms  in  the 

Integrand  as  a  polynomial  in  w.  One  then  obtains  the  expressions 

treated  in  Appendix  C.  For  n  <  2,  one  must  carry  out  the  division 

which  gives  a  polynomial  (of  rather  low  order  in  w)  and  remainder 

2  1/2 

terms  with  denominators  (wcosa  -  slna)(w  +  1)  and 
2  2  1  /2 

(wcosa  -sina)  (w  +  1 )  .  Also  these  integrals  are  treated  in 

Appendix  C.  For  the  case  m  -  1,  n  -  0,  one  has 

2  "2 
(wslna  ♦  cosa)  (wcosa  -  sina) 

2  —2  —I  — 2  — ? 

-  tg  a  +  2slna  cps  a(wcosa  -  sina)  +  cos  a(wcosa  -  sina) 


i,l+1 j02 
■‘‘mn 


Then  with  Eqs.  (C.16)  and  (C.17) 


sign(Uj)Uj^{tg^a  log(<w^  +  *  w) 

—2  ~  1 
♦  3ino  C03  a3ign(co30)log[3in( e-a) ( 1  ♦co3( 0-o) )  ] 

-C03  ^ a3ign(co30)3ln  ^(0-a)} 

The  3ign  functlone  cancel  because  of  Eq.  (133).  as  can  be  seen  from 
the  following  discussion 

sign  (Uj^)log((w^  +  1)^'^^  w)  - 

(sign(Uj^)log[(u^^  ♦  (signuj^  )v)/ |Uj  I  ] 

The  sign  function  can  obviously  be  omitted  if  sign  (u^^)  -  +1. 
Consider  now  sign  (u^^)  -  -1.  Then  one  obtains  for  the  right-hand 
side 

-log[(u^^  ♦  -  V)/|Uj^|]  -  -10g[Uj^/[  |Uj^l  ( (Uj^  +  v)]] 

-  log[((Uj^^  +  v)/iuj^| 


Therefore , 

(signUj)log(  (w^  +  1)^^^  +  w)  -  logCCu^^  ♦  +  v)/|Uj^l] 

One  thus  obtains  (with  a  change  of  the  constant  of  integration) 


l,i+1 ^02 


j  ^Uj{-tg^a  log[k((u^^ 


+  +  V)] 


(145) 


-2  - 1 
-  sino  cos  alog[( n-n. )/(R  ♦  5-x)]  +  cos  aR/(n-y)} 


One  has  for  the  case  m  ■  Ot  n  ■  1 

(wsina  +  cosa)(wcosa  -  sina)  ^  ■  tga  +  cos  ^o(wcosa  -  sina) 


Therefore 


-  +  0^u^{-tgalog((R+v)  |Uj^  I  b 


-  cos  log((n-y)[R  +  (t-x)]  ') 


i+1 


-1 


(146) 


■ 


Collecting  these  results  (Eqs.  (144),  (145),  and  (146))  and 
replacing  in  by  cos  U  one  obtains 


I°J(X,Y.a)  -  6^{cos'‘'a  (U^/Y)  -  sina  cos  ‘o  U  log|kY|} 


.2,  -2  ...2  .  .  ...-2 


(147) 


lQ](X,Y,a)  -  -B^cos'^a  U  log|kY| 


(In  these  equations  a  term  log  (cos  a)  has  been  omitted  because  it 
can  be  incorporated  into  a  constant  of  integration.) 


I°Q(X,Y,a)  -  6^u{tg^a  log(  (R-V)/ |U  | )  -  sina  cos  “"a  log(  (R-X)/ |  Y  |  ) 

(148) 


-2 


+  cos  ^a  R/Y) 


lQ^(X,Y,a)  -  B^U{tga  log( (R-V)/ (U | )  -  cos"’a  log( (R-X)/ | Y | ) 

See  Eqs.  (113)  and  (140)  for  the  definitions  of  the  variables. 

For  the  power  1  of  the  reduced  frequency  k,  one  starts  from 
Eqs.  (104).  Introducing  the  rotated  system  of  coordinates 
Eqs.  (124)  through  (129)  one  obtains  for  a  /  ii/2 


74 


I 

k  •  • 

.nV 

►  •  W 
.  ^ 

.  *-  « 


M 


yxox 


1 

mn 


-ik(m+n+1 ) 


-1 


2v-1/2 


I  (Uj  +  v^)  '^^(Uj,cosa  ♦  V  sino) 


m 


(-UjSina  +  cosa)”uj^ 


dv 


and  with 


w  -  v/u. 


an 


. ,  /  ,  V -1  m+n+1 

-ik(m+n+1 ) 


i+1 


(w^  ♦  l)^^^(wsina  +  cosa)”* 


w. 


(wcoaa  -  3ino)''dw 


Hence 


iJqCX.Y.o)  -  ♦ik  U  log(R-V)/|Ul 

iJqCX.Y.o)  -  +i(k/2){3ino  UR  -  co3a  U^  log( (R-V)/ |U | )  (149 

lJ^(X,Y,a)  -  +i(k/2){co3a  UR  ♦  3ina  U^  log[(R-V)/ |U | ] 

2  2 

Eq.  (105)  i3  the  baeic  formula  for  terma  of  order  k  and  k  logk. 
We  write  aa  before 


-21 

i3  treated  in  the  rotated  coordinate  ayatera.  The  reaulta  can 
mn  ' 

be  obtained  in  the  "aheared"  ayatem  of  coordinatea,  but  thia 
complicatea  the  limiting  proceaa  a-»'iT/2. 


i,i+1j21  _  (k^/2)C(m+n+2)~^  f  (sinav  +  cosau,)® 
mn  J  ^ 

'^i 

(coaav  -  sinaUj  )%j^dv 

Hence, 

IqJcx.y.o)  -  (k^/l»)C  UV 

I^J(X,Y,a)  -  (k^/6)C{-coso  U^V  -  (1/2)3ina  UV^}  (150) 

Iq](X,Y,o)  -  (k^/6)C  +  (sino  U^V  -  (1/2)co3a  UV^) 

I  i3  treated  in  the  rotated  ayatem  of  coordinatea.  With  w  -  v/u^ 
one  obtaina 

''i  +  l 

-  -0^(k^/2)  (m>n+2)“^uj"’”''^  J  (w^  +  1)“^^^ 

^'i 

(waino  +  coao)”*^^  (wcoao  -  3ina)’^dw 

One  haa 

for  m  -  0,  n  -  0 

p  _ 1 /p  p  1 /2 

j(w  +  1)  '^(waina  +  cosa)dw  -  3lna(w  +  1) 

+  coaa  logC(w^  ♦  1)^^^  +  w], 

for  m  -  1 ,  n  -  0 

p  _  1  /p  p  2  2  1/2 

|(w  +  1)  '^(waina  +  coaa)  dw  -  (l/2)3in  aw(w  +  1) 

2  1  /2  2  2  2  1  /2 
23inaco3a( +  1)  +  (coa  a  -  (1/2)3in  a)log[(w  +1)  +  w]it 


and  for  m  -  0,  n  -  1 

p  —I  /p 

|(w  +  1)  (wsina  +  coaaXwcosa  -  sina)dw 

p  -1  /2  2  2  2 

-  (1/2)sino  C03aw(w  +  1)  +  (cos  o  -  sin  a)(w  +  1 

1  /2 

-  (3/2)sina  cosa  logC(w  +1)  +  w] 

Therefore 

Iqo  “  (sina  UR-  cosa  log(  (R-V)/ |U  | ) } 

ifo  “  -6^(k^/6) {(1/2)sin^a  UVR  +  2sina  cosa  U^R 

-(1/2)(3  cos^a  -  1)U^  log(R-V)/|U| ) 

^01  "  -6^(k^/6) {1 /2  sina  cosa  UVR  +  (cos^a  -  3in^a)U^R 

*  (3/2)sina  cosa  U^  log( (R-V)/ |U | ) 

Using  again  the  rotated  coordinate  system  one  obtains  from 
Eq.  (104) 

i,i+1j23  .  (k^/2) (m+n+2)~.  f  (u.cosa  +  vsina)®(-u.sina  + 
mn  J  1  1 

[log(k(R  +  u.cosa  ♦  vsina))  -  (m+n+2)  ^u,]dv 


1  /2 


(151 


vcosa) 


(152 


and  with 

w 

w 

where 

Let 


Then 


1,1+1=23 

00 

i,l+1r23 

10 


i.i+1f23 

01 


One  has 


Here 


w  -  V/Uj^ 

( 


1  +  1 

I  (wslna  +  coso)“(wcoaa  -  3ino)'^Cf(w)  -  (m+n+2)~^]dw 
1 


f(w)  -  log(kUj)  ♦  logC(w^  ♦  ♦  wslna  +  coso] 

Qj  -  I  w'^f(w)dw  ( 

-  (k^/i»)[u^^(QQ  -  (l/2)vuj]j^’ 

-  (k^/6)[cosa  Uj^Qq  ♦  slna 
-  (1/3)coso  Uj^v  -  (l/6)slno 

-  (k^/6)[-slnaUj^  +  coso 

♦  (1/3)slno  -  (l/6)coso 

Qj  -  ( J  +  1  )’kw'^^^f(w)  -  j  '(w)dw] 

+  wslna  ♦  cosa]~^[w(w^  + 


f'(w)  -  [(w^  + 


+  slno] 


Multiplying  numerator  and  denominator  by  (w  -^1) 
cosa)  one  obtains 


f  C. 


(wsina  + 


f  ..p  2  1  /2 

f  (w)  -  [wcosa  -  slna]  [(w  •‘■1).  -  (wsino  +  cosa)] 


[w  ♦  (w^  ♦  1  ^^sina](w^  1)  ^ 


I  2  -1  /2 

f  (w)  -  [wcoso  -  sina]  [cosa  -  (w  *1)  ] 


Then,  for  j  -  0,  after  dividing  w  by  [wcosa  -  sina] 


I  2—1  /2 

jwf  (w)dw  -  |[co3a  .  ♦  tgaCwcosa  -  sina)  ][cosa  -  [w  +  1 ) 


and  with  Eq.  (C.16) 


t  *^12  1  /2 

Jwf  (w)dw  -  w  +  tga  log(wcoaa  -  sina)  -  cos  a  log((w  +1) 


“1 


-tga  log[3in(0-a)  (1  cos(0-a))  .] 


here 


0  -  arctgw 


Moreover,  for  j  -  1,  again,  after  carrying  out  the  division 

tr _ —  1  .  _ _ —2  . 2  / _  \— 1 


jwf  (w)dw  -  j[wcos  a  ♦  sina  cos  ^ <x  ♦  tg‘^a(wco3a  -  sina)  ] 


[cosa  -  (w^  ♦  1)  ^^^]dw 


,  2  ’  2  2 

jwf  wdw  -  (w  /2)  +  tga  w  ♦  tg  a  log(wcosa  -  sina) 


-1  ?  1  /?  ?  1  /? 

-  cos  'a(w  +1)  -  sina  cos  log((w^  +1)  ♦  w) 


-  tg^a  log[sin( 0-a)  ( 1  cos(0-a))  ^] 


By  substitution  of  w  -  v/Uj^  one  obtains 

p  PI  /P 

f  -  log(k|Uj^|)  +  log[((Uj  ♦  V  ) .  +  Ujcosa  +  V3ina)/|uj|] 

f  -  log(k(R  +  £-x)) 

j  wf  (w)dw  -  v/Uj  ♦  tgo  log[k(n-y)]  “  cos  log[  (R+v )/ |u |  ] 

-  tga  log[(n-y)(R  +  t-x)  ^3 

I  wf  (w)dw  -  v/Uj^  -  cos~lo  log(  (R+v)/ |u  I )  +  tga  log(k(R  ♦  (£-x)) 

J  w^f*(w)dw  -  <1/2)(v/Uj)^  +  tga(v/Uj^)  -  cos  ^o(R/Uj^) 

-  sinacos  log[(R+v)/ ju^^  1 3  ♦  tg^o  log[k{R  +  (5~x))3 

Qq  -  (V/U  -  tgo)log(k(R-X))  -  V/U  +  cos’^o  log((R-V)/|U| ) 

-  (1/2){(V/U)^  -  tg^a)log{k(R-X))  ♦  sina  cos'^o  log( (R-V)/ |U |  ) 
-  (1/2) (V/U)^  -  tga(V/U)  -  cos’^aCR/U)} 

and  rewriting  Eqs.  (155) 

lQ^(X,Y,a)  -  (k^/4){U^QQ  -  (1/2)VU)} 

I^^(X,Y,a)  -  -  (k^/6){cosa  U^Qq  ♦  sina 

-  (1/3)cosa  U^V  -  (1/6)sina  UV^} 

lQ^(X,Y,a)  -  (k^/6){*sina  U^Qq  -  cosa 

-  (1/3)3ina  U^V  +  (1/6)cosa  UV^} 
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SECTION  VI 

INTEGRATIONS  WITH  RESPECT  TO  X  AND  Y 

With  the  formulae  derived  so  far,  the  upwash  at  a  given  point 
x,y  can  be  evaluated  in  terms  of  the  parameters  that  describe  the 
pressure  distribution.  One  still  has  to  carry  out  the  integration 
over  the  xy-elements.  This  can  be  done  numerically  or 
analytically.  The  upwash  has  singularities  as  the  point  x,y 
approaches  the  boundaries  of  the  5n~element  or  of  its  wake,  and 
special  provisions  must  be  made  for  the  terms  which  give  infinities 
in  the  upwash.  Aside  from  this  one  will  obtain  good  results  if  one 
uses  a  sufficiently  large  number  of  points  xy. 

The  analytic  procedure  will  be  shown  in  this  section.  In  a 
numerical  approach  one  will  add  all  (or  nearly  all)  contributions 
to  the  upwash  before  the  integrations  are  carried  out;  in  an 
analytic  procedure  one  must  keep  them  separate.  This  gives  quite 
lengthy  lists  of  formulae. 

The  upwash  formulae  have  arisen  from  contour  integrals  around 
the  Cn-element.  The  individual  expressions  are  functions  of 
X  -  x-5^,  Y  -  y-hj  and  variables  U  and  V  are  considered  as 

functions  of  X  and  Y.  The  and  are  the  corners  of  the  in¬ 
clement,  gives  the  slope  of  the  side  of  the  in-element  for 

which  the  integration  has  been  carried  out.  The  general 
expressions  (for  a  ^  ir/2)  fail  if  one  tries  to  substitute  a  -  ir/S. 
Separate  formulae  for  a  -  11/2  have  therefore  been  derived.  The 
limiting  process  a-»ir/2  is  shown  in  Appendix  F. 

The  contribution  to  the  upwash  for  the  form  of  the  pressure 
distribution  assumed  here  is  given  by 

p(x,y)lQQ(X,Y,a)  +  c^Ql^Q(X,Y,a)  +  Cq, Iq^ ( X , Y , o) , 


but 


p(x,y)  -  p(ij,nj)+(x-i^)c^Q+(n-yi)cQ^ 
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Thus,  one  has  to  evaluate  expressions 


pCSihi)  jj  Ioo<3xdy  +  c^Q  jj  ^’56) 


with 


^10  "  ^10  ^  ^^00 


(157) 


^01  ■  ^01  ■‘^  ^^00 


The  integration  is  carried  out  over  an  xy-element.  The 
expressions  I  consists  of  a  number  of  summands  all  of  which,  except 
those  that  contain  a  factor  log  (k(R-X))  appear  in  the  form. 


G(x,y)  -  R’^fO) 


(158) 


if  written  in  polar  coordinates  R,  0.  Inspecting  the  expressions 


one  finds,  that 
mn 


t  -  (I  ♦  m  +  n) 


The  Integral  over  an  xy-element  then  assumes  the  form 


||G(x,y)dxdy  -  JJ  R^"’^ f.( 0 )dRde 

Here  the  Integration  with  respect  to  R  can  be  carried  out;  the 

contribution  of  the  lower  limit  R  -  0  vanishes  because  t  >  -2; 

_2 

subsequently,  dO  is  expressed  by  d0  -  (x  dy  -  ydx)R  .  Therefore, 

J  G(x,y)dxdy  -  (1*2)"^  G(x,y)(xdy-ydx)  (159) 

where  the  integral  on  the  right-  is  taken  along  the  contour  of  the 
xy-element. 

The  expression  containing  the  factor  log(k(R-X))  has  the  form 


G(x,y)  -  g(x,y)log(k(R-X) ) 


where 


g(x,y)  -  R^f(8) 


Therefore, 

G(x,y)  -  R^f(0)[log(kR)  ♦  log((R-X)/R)] 

The  second  term  in  the  bracket  depends  only  upon  e.  Therefore, 

J|  G(x,y)dxdy  -  ||  R*'*^  f  (0)[log(  (R-X)/R)  +  log(kR)]d0 

-  ^  f(0){[R*"‘^[log(R-X)/(R)  +  log(kR)])  -  dR) } d0 

0 

Hence, 

j|  G(x,y)dxdy  -  (1+2)”^^  [G(x,y)  -  (t+2)"’ g(x,y)] (xdy-ydx)  (160) 

We  Introduce  an  angle  Y  which  gives  the  direction  of  one  side  of 
the  xy-element. 


tgY 


i,  i*l 


Xj^-X 


Yizy 


ill 

i^-1 


(161  ) 


For  integrations  which  are  to  be  carried  out  along  the  sides  of  the 
xy-eleraent  we  introduce  coordinates  p,  q; 


X  -  p  cos  Y  ♦  q  sin  Y  ;  p  -  X  cos  Y  -  Y  sin  Y 
Y  -  -p  sin  Y  ♦  q  cos  Y;  q  -  X  sin  Y  ♦  Y  cos  Y 


(162) 


Along  such  a  side  p  is  constant;  the  variable  of  integration  is  q. 
The  expressions  U  and  V  had  been  defined  by 


X  -  U  coa  a  V  sin  a  ;  U  -  X  cos  o  -  Y  sin  a 

Y  -  -U  sin  a  +  V  cos  a  ;  V  ■  X  sin  a  +  Y  coa  a 
For  the  integrations  U  and  V  are  expressed  by  p  and  q 

U  -  p  coa  6  +  q  sin  6  ;  p  -  U  cos  6  -  V  sin  6 

V  -  -p  sin  6  +  q  coa  6  ;  q  -  U  sin  6  +  V  cos  6 


(163) 


(164) 


with 


6  -  Y  -  o 


(164a) 


In  the  integration  along  the  sides  of  the  xy-elements  expressions, 
xdy-ydx  will  occur.  One  finds  in  general 

XdY-YdX  -  pdq-qdp  (165) 


and  since  p  -  const 

XdY-Ydx  -  pdq 

Now  we  summarize  the  results  of  Section  V  (i.e.,  the  formulae 
for  the  functions  Iqq  and  I^g  and  Iq^ )  in  a  form  suitable  for  the 
integration.  For  this  purpose  we  introduce 

♦  ^  -  log  (1<(R-X))  -  log  (R-X)/|Y|)  +  log  (k|Yl) 


log  (R-V)/|U|) 


k  w-j 


•  ' 


These  expressions  will  be  regarded  as  functions  of  p  and  q.  The 

expressions  I’’  are  linear  combinations  of  these  functions  with 

•  • 

cofficients  which  are  homogeneous  in  U,  V,  X  and  Y.  (The  explicit 
dependence  upon  X  and  Y  occurs  because  we  have  to  form 

^10  “  ^10  *  *^00  ^01  “  ^01  *  ^^00*^  These  coefficients  are 

considered  as  scalar  products,  written  as  products  of  row  matrices 
and  column  matrices.  The  row  matrix  depends  only  upon  a,  the 
column  matrices  are  homogeneous  functions  in  (U,V,X,Y).  In  the 
following  expressions,  only  those  terms  actually  encountered  in  the 
expressions  I  are  included.  We  introduce 


[w,]  -  [U,X,Y]* 

[W2]  -  [U^,UV,UX,UY]* 


(167) 


Cw^]  -  [U^,U^V,UV^,U^X,UVX,U^Y,UVY]'" 

The  row  matrices  belonging  to  a  function  lf„(or  I^^)  are  denoted  by 
c  with  the  same  indices  as  the  function  I  and  a  third  subscript 
referring  to  the  function  to  which  they  belong.  Notice  that  the 
coefficients  occurring  in  a  specific  function  I  may  be  of  different 
degree,  because  the  functions  ♦j  have  different  degrees  in  U  and  V. 
Only  some  of  the  function  ^  will  occur  in  a  given  expression  I. 

For  Iqq  the  matrix  [Wq]  and  the  vectors  Cqq,..  are  scalars. 
These  expressions  are  written  down  directly.  One  then  obtains  the 
following  list. 


2  -1 
0  tga^^  -  cos  0(>2 


(168) 


m+n  -  1 


(cont'd) 


^°00,1  ^*-“2^**’l  *  ^°00,2^^“2^*^2  *  ^°00,5^^”2^%  ■*■  *- °00 , 6 


The  vectors  [c],  ordered  according  to  their  dimensions  (which,  of 
course,  matches  the  dimension  of  the  vectors  [w]),  are  given  by 


0  2  ~2  2 
[c^Q  sinacos  a,  $  tga,  0] 

r  0  T  r«2^  2  „2  -1 

[Cio,2J  “  LB  tg  o,  -B  cos  a,  0] 

Cc°o,n3  “  CB^cos'^a,  0,  0] 

[Cqi^i]  -  [-B^cos’^o,  0,  B^tga] 

0  2  2  -1 

^°01,2^  -  CB  tga,  0,  -B  cos  a] 

CcJo,2^  “ 

[c|o,6^  -  [(il</2)sin  a,  0,  0] 
*-°01,6^  -  [(ik/2)cos  o,  0,  0] 

[Cqo,6^  -  [ (k^/M) B"^sin  o,  0,  0] 

0  2  -2 

[c^0,3^  -  [B  cos  a,  0,  0,  0] 

[c^0,2^  -  [-(ik/2)cos  a,  0,  ik,  0] 

Ccl,  ,]  -  [(ik/2)sin  o,  0,  0,  ik] 


]  -  [-(l<"/4)tga,  (k^/4),  0,  0] 


^°00,1 

[Cqo  2^  "  (I<^/4)(co3  ^a-e  ^cos  a),  0,  0,  0] 

[c^0,5^  -  [0,  (k^/4)(C-3/2),  0,  0] 

"  [(k^/6)((l/2)tgo-28‘^slna  cosa) ,  (k^/12) T^sin  o, 

’  (k2/4/B-23ina.O] 

"  C  (k^/6)  ( (l/2)-8  ^(cos^a-3in^o)  ,-(k^/2)  6  ^3inoco3a,  0, 
(k'^/4)0  ‘^sina] 

[(k^/6)3ln  a  (l+(l/2)tg^a) ,-(k^/6)co3  a,  -(k^/12)3ina, 

-  (k^/4)tgoi.  k^/4,  0,  0]  (169) 

(cont*d) 

[(k^/6)(-(l*(l/2)tg^o)+{l/2)r^(3co3^o-l),0,0, 

(k^/4) (cos"^o-6~^co3  a),  0,  0,  0] 

[CiQ  gD  -  [0,{k^/6)((4/3)-C)cos  a^(l/2)3lnatgo) , 

(k^/6)((5/12)-(C/2))sin  a,  0,  (k^/4) (C-3/2) ,  0,  0] 

[c^.  i]  -  C-(k^/12)3inotga,  (k^/6)sin  a,  -(k^/12)co3  a,  0,  0, 
-(k2/4)tg  a,  (k2/4)] 

[Cq^  2^  “  t~('<^/12)  (tgo-3B  ^3inaco3a,0,  0,  0,  0, 

-(k^/4)tg  a,  <k^/4)] 

tCn?  r]  “  [0,(k^/6)(C-5/6)3ln  a,  (k^/6) ( (5/12)-(C/2)cos  a,  0,  0, 

U 1  f  D  ^ 

0,  (k‘^/4)C-3/2)  ] 


'■“10.2^  ■ 


The  factor3  tga  and  cos  which  occur  in  some  of  the  coefficients 
are  indications  that  some  of  these  formulae  are  not  applicable  for 
a  -  ir/2.  The  results  necessary  for  a  -  ir/2  have  been  obtained  by 
direct  computation  in  Section  V,  They  are  brought  into  a  form 


analogous  to  that  for  a  ir/2.  The  functions  ^  remains  the  same. 
The  function  does  not  appear.  The  vectors  [Wj]  assume  different 
forms 


Cw,]  -  [X,Y]\  [Wg]*  -  [X^.XY.Y^]",  [W3]  -  [X^Y.XY^.Y^]^  (170) 

Then 


1^0 


Iqi  (ir/2) 


looCif/S) 


ljQ(ir/2) 


Iqi  (ir/2) 


lJo<’f/2) 


[e^,o][w^  34.3  + 

[0,a^/2][w^]^^ +[36^/2, 0,0][w2](fr3  +  [3e^/2,0][w3^](l>4 

[2B^,0][w^  +  2B^*g 

[0,  -ik][w^](|»^ 

(171) 

[0,  -ik,  0][w2]«|»i  ♦  C0,-ik/2][w^  ](j»g 

[0,0,  -ik/2][w2]«i 

[0,-k^/l»,0][w2]<fri  ♦  C0,(k^/4)(1-C)][w2]^5 

+  [0,(k2/4)(1-B"2)][w^]4>5 


I^Q(it/2)  -  [-k^/6,0,(k^/12)(B"^ 

♦  (1/2))][w3]*^  +  [7k^/72,-C/6,0,0][w^](fr5 

+  [0,k^(-(B'2/6)-(5/24),0][w2]^6 

Iq^(v/2)  -  [0,-k^/12,0][w3]<fr^  ♦  C0,0,k^(-(C/12)  +  (5/1 8) )  ][W3]<>5 

♦  [0,0,-(b"^/12)  -  (1/12)][w2]^6 


So  far  this  amounts  only  to  repetitions  in  a  different  form 
of  the  formulae  of  Section  V.  To  carry  out  the  integration  along 
the  side  of  the  xy-elements  for  the  case  a  -  ir/2  one  must  express 


the  functions  U,V,X,  and  Y,  occurring  in  the  vectors  Cwj],  by  p, 
and  q.  Along  the  path  of  integration  p  is  constant  and  q  is  the 
variable  of  integration. 

Therefore,  one  writes  for  a  ^  ir/2 

-  Mj^Cp.q]^ 

CW2]  -  M2Cp^,pq,q^]* 

[w^]  -  M^Cp^.p^q.q^]"" 

where  ,  and  M2  are  respectively  3  by  2,  4  by  3  and  7  by  4 
matrices.  The  rows  of  the  matrices  are  the  coefficients  obtained 

by  expressing,  one  element  of  the  column  matrix  on  the  left  in 

2 

terms  of  p  and  q.  (One  has  for  instance  in  [w*]  U 

22  22  ^ 

cos  6p  +  2coa6sin6pq  ♦  sin  6q  )  and  one  will  indeed  find  in  Mp  as 

2  2  ^ 
first  row;  cos  6,  2cos6sin5,  sin  6.  One  obtains  the  following 

expressions. 


M 


1 


cos  6  sin  fi 

cos  Y  sin  Y 

sin  Y  cos  Y 


(172) 


M  - 

2 


2 

COS  6  2sin6cos6 

2  2 

-sinScosfi  cos  6-sin  6 

cos6cosY  sin6cosY-*-oos6sinY 


3in^6 

sin6cos6 

sin6sinY 


(173) 


-3in6cosY  cos6cosY-3in6sinY 


3in6co3Y 
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is:  >  s: 


For  a  -  ir/2  one  sets 


■*  [M^][p,q]* 

Wg  -  [MgDLp^tPqtQ^]* 

-  [MjlCp^.pq.q^]^ 

The  rows  of  ,  M2,  and  are  the  coefficients  of  the  development 
of  X,  Y,  X^,  XY  etc.  in  terms  of  p  and  q. 
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Take  for 


Now  consider  some  term  of  one  function  . 
instance  a  term  containing  the  column  matrix  [W2].  Then  one  has 

[c  JCw-]*.  -  {[c  .]CM_]}{[p^,pq,q^]%.} 

The  expression  within  the  first  brace  on  the  right  is  a  row  vector 
which  depends  solely  on  a  and  Y.  The  integration  along  a  side  of 
the  xy-element  operates  only  on  the  elements  of  the  vector 


Applying  Eq.  (129)  and  observing  that  xdy-ydx  -  pdq  one  obtains  as 
a  contribution  to 

mn 


For  -  log(k(R-X)  one  has  to  apply  Eq.  (130);  the  function  g(xy) 
'  2 

is  given  by  1,  q,  and  q  .  One  obtains 


( l)  2  ■ 

We  denote  by  the  vector  C^^dq,  q^^dq,  q  ^.dq...] 

j  th  J  J  J 


truncated  after  the  1 


term. 


The  contribution  to  the  Integral  over  the  xy-element  of  one 
endpoint  of  one  aide  of  the  xy-element  with  angle  Y  is  denoted  by 
J**(atY),  with  the  same  indices  as  the  functions  1‘*  from  which  it 


originates.  Ultimately,  J*‘  depends  on  X 

•  • 

(besides  a  and  Y) .  Then  one  has 


and  Y 


Yj-tli 


*  P‘ 


Ci|»j],j-2  (178) 


i  i 

^  is  then  the  sum  of  all  pertinent  expressions  (See 

Eqs.  (168)  and  (171).)  The  dimensions  of  the  row  matrix  [w],  the 

diagonal  matrix  containing  the  powers  of  p  and  of  the  vector 

have  not  been  shown.  The  subscript  of  w  is  found  in  Eq.  (168). 

The  number  of  rows  in  w  gives  the  dimension  of  the  following  square 

matrix  and  of  the  vector  C^].  For  J  -  1  one  has 


ffln,1 


m+n+i 


m 

p  ♦n+l-l 


r  ,  m+n+tn 

J 


(179) 


As  preparation  for  the  evaluation  of  the  vectors  we  first 

list  some  recurring  auxiliary  expressions.  In  the  computation 
they,  and  also  the  vector  [i|i],  can  be  evaluated  immediately  as 
definite  integrals  between  two  limits  for  q. 


j  R"’dq  -  -log  (k(R-q)) 

I  R  ^  qdq  ■  R 

I  R'^^dq  -  ^  ((Rq*p^  log(k(R-q)) 

J  R'^^dq  -  R[(1/3)q^  -  (2/3)p^] 

(180) 

I  R'^q^dq  -  R[(lA)q^  -  (3/8)p^q]  -  (3/8)p**  log(k(R-q) 

J  Rdq  -  1  (Rq  -  p^  log(k(R-V)) 

J  Rqdq  -  1 /3  R^ 

I  Rq^dq  -  RC(1/H)q^  ♦  (1/8)p^q]  ♦  1/8  p^  log{k(R-V)) 

/  bI  ■  /  »(qco3V%3lnYT  ’  ?  -  .  i  lcgCR-X)/Y) 

The  last  relation  follows  from  Eq.  (E.2). 

The  general  formula  for  j  q'°^^dq  is  Eq.  (E.3)>  Here  we 
specialize 

j  ♦^dq  -  l(q  -  ptgY)log(k(R-X))  -  pcoa~S  |  R~^dq  -  q} 

j  ♦^qdq  -  ^  l(q^  -  (ptgY^ )log(k(R-X) )  -  pcos'N  j  R~^  (q  +  (ptgY)  )dq 
-(q^/2)  -  (ptgY)q} 

J  ♦^q^dq  -  (1/3)  l(q^  -  (ptgY) ^log(k(R-X) ) 

-  pcos'^Y  I  R~^(q^  ♦  (ptqY)q  ♦  (ptgY)^)dq  (181) 

-  (q^/3)  -  (ptgY)(q^/2)  -  (ptgY)^} 

J  ♦^q^dq  -  (1/4)  iCq**  -  (ptgY)  )^log(k(R-X) ) 

-  pcos  S  J  R  ^(q^  ♦  (ptgY)q^  ♦  (ptgY)^q  ♦  (ptgY)^)dq 

-  (q^/4)  -  ptgY(q^/3)  -  (ptgY)^(q^/2)  -  (ptgY)^q} 
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From  the  factor  of  log  (k(R-X)),  that  is  from  q*'-(ptgY)'^,  one  can 
spilt  off  a  factor  (q-ptgY)  -  Ycos  Furthermore  from  Eqs.  (E.4) 

I  “  (-pcot6)log{  (R'V)/U)  +  psin  ^6  J  R  ^dq} 

I  -  (1/2)  {(q^  -  (-pcot6)^)log(  (R-V)/U) 

♦  pain”^  6  R~^(q  ♦  (-pcot6)dq}  (182) 

J  ■  (1/3)  l(q^  -  (-pcot6)^)log((R-V)/U) 

♦  psin  6  j  R~  (q  ♦  (-pcot6)q  ♦  (-pcot6)  dq} 

J  ■  <1''^H(q^  -  (-pcot«)^)log(R-V)/U) 

■“1  *^13  2  2 

+  psin  6  J  R  (q"^  ♦  (-pcot6)q  ♦  (-pcotfi)  q 

+  (-pcot6)^)dq} 

fi  n  "*  1 

From  q  -  (-pcot6)”  a  factor  q  ♦  pcos5  -  Usin  6  can  be  split  off. 

I  ^Ylog(kY) 

I  ♦^qdq  -  co3~S[q  +  (ptgY)log(kY)]  (183) 

j  N[(q^/2)+(ptgY)q+(ptgY)^log(kY)] 

j  N[(q/3)+(ptgY)(q^/2)+(ptgY)^q  +  (ptgY)^log(k(Y))] 

j  ♦4<lq  “  J  RY'^dq  -  p^cos'^Y  |  R~’ Y~’ dq*cos~S[  |  ^  dq+(ptgY)j 

j  ♦ijqdq  -  p^tgYcos'^j  R~^Y  ^dq+cos  N  j(R  \^  +  (ptgY)R  \)dq 

♦  p^cos'^Y  I  R~^dq 

j  ♦sq'^dq  -  ^  q"*’  (n  -  0,  1,2,  3) 

I  ♦eq'^dq  -  j  q^Rdq,  n  -  0,  1,  2 
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One  obviously  needs  special  formulae  in  those  integrals  where 
cos  and  tga  or  sin  ^6  and  cos  5  occur,  if  respectively  or 

6^0.  They  are  found  by  direct  integration.  For  Y  -  ir/2  one  has 
q  -  X,  p  -  -Y.  Then 

J  ♦^q"dq  -  I  x"  log{k(R-X))dx  -  -jjJj  [x"*’  log(k(R-X))  ♦  j  X^^^R'^dx] 

The  integrals  are  Immediately  found  from  Eqs.  (180)  by  replacing  q 
by  X,  and  p  by  -y.  Therefore,  for  Y  -  ir/2 

J  ♦^dq  -  X  log(k(R-X)  ♦  R 

I  ♦^qdq  -  1 /2{[log(k(R-X))  ♦  j  (RX  -*■  Y^log(k  (R-X) ) } 

I  ♦^q^dq  -  (1/3) {[X^log(k(R-X))  ♦  (1 /3)R(X^-2Y^) }  (184) 

I  ♦^q^dq  -  l/4{X**log(k(R-X))+R(l/4)X^-(3/8)XY^)-(3/8)Y^log(k(R-X))} 

I  ♦3q”dq  -  y'^  j  x”dx  -  Y‘\n*l)'^ x”^^  ,  n  -  0,  1,  2,  3- 

I  ♦4q"dy  -  y'^  j  X^Rdx 
I  ♦jjdy  -  Y'U(RX-Y^log(k(R-X))] 
j  ♦4qdy  -  y"‘'(1/3)R 

The  limit  6-0  requires  a  revision  of  the  integrals  Involving  ♦2* 
Then  p  -  U,  q  -  V,  and  one  obtains 

J  q"log(R-V/|U| )dq  -  |  v"log( (R-V)/U)dv 

-  (n+1  )"^[v"'"^log(R-V)/U)  +  j  v”*^R"’dV 

j  ♦2dq  -  {Vlog(R-V/U)  +  R  (185) 

I  q^gdq  -  (1/2) {V^log((R-V)/U)  ♦  1/2((RV  ♦  U^log( (R-V) /Y) } 
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SECTION  VII 

SINGULARITIES  OF  THE  UPWASH  FIELD 


In  carrying  out  the  integrations  over  X  and  Y  one  must  be 
aware  of  the  singularities  in  the  upwash  field  and  how  different 
terms  contribute  to  them.  Consider  a  single  tn-element. 
Singularities  will  probably  occur  at  its  boundaries  and  at  the 
boundaries  of  the  wake  and  at  the  wake  of  points  of  the  contour  of 
the  ^n-element  for  which  the  contour  is  discontinuous,  (see 
Figs.  3).  The  wake  has  already  been  studied  in  Section  III,  but 
one  must  ask  how  these  wake  singularities  express  themselves  by  the 
formulae  found  in  Section  V. 

Because  of  the  denominator  Y  the  strongest  singularities  are 
caused  by  the  functions  and  Assume  first  that  Ap  is 

constant  throughout  the  element,  then  only  Iqq  will  contain 
<>ij.  One  has  as  the  only  contribution 

iQQ(ir/2)  -  0^y‘\x+r)  -  B^(Y-ni)‘^C(x-tj)  +  ((x-Ci)^*(y-ni)^)’^^] 

For  the  triangular  element,  shown  in  Fig.  3,  this  must  be 
evaluated  for  the  side  2,3.  To  recognize  the  analytic  behavior,  we 
rewrite  this  expression.  One  has  for  x  <  5^  ^  ^2 

y"^ (X+R)  -  Y/(R-X)  -  Y/2R 

Here  the  flow  field  does  not  have  a  singularity.  One  has  for 

X  >  ^2  >  ^3 

Y"hx+R)  -  y‘’(2X  +  Y(R-X))  -  2XY~^  for  Y  small 

This  displays  the  singularity  caused  by  one  of  the  limits 
We  know  from  Section  III,  that  except  for  a  factor  exp(ikx), 
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the  singularity  is  independent  of  x.  This  is  seen  if  one 
substitutes  the  limits  ^3)*  replaces  X  by  X-C2 

x-^^*  One  then  obtains 

2y"’[(x-53)  “  ■  2y~^(^2"53^ 

The  situation  is  more  complicated  if  Ap  is  linear  in  ^  and  n 
The  contribution  of  ®nd  then  appear 

in  Iqq  (it/2),  namely,  6^  Y  ^X+R) 

in  I°Q  (x/2),  namely,  (  B^/2)XY"\x+R) 

0  2-1  -2  2  -1 
and  in  I^q  (x),  namely,  B  Y  (cos  oU  ♦  cos  oUR) 

The  last  term  is  rewritten  for  the  vicinity  of  Y  -  0.  There 
U  -  X  COSO.  Therefore,  except  for  terms  of  higher  order  in  Y 

6^Y“^  (cos"^oU^  ♦  cos'^oUR)  -  B^y"^X(X-*-R) 

As  before,  no  singularities  will  be  encountered  for  x  <  ^2  ^  ^3* 

In  the  following  discussion  for  x  >  ^2  ^  ^3  w®  retain  only  the 
singular  parts. 

B^Y"’ (X+R)+  2B^y"’x 

2-1  2-12 
B  /2Y  X(X+R)-  BY  X 

B^Y"^  (cos'^oU^  +  cos'^oUR)-*’  2B^y"’x 

Now  choose  a  function  Ap(C.n),  whose  integral,  taken  for 
n  -  n2  •  ^3  from  to  is  zero.  Such  a  function  should  not 
generate  a  wake  singularity; 


Ap( n,) 


c,nU-(5p^5,)/2) 


In  order  to  evaluate  the  upwash,  one  must  first  evaluate 
Ap(x,y)  -  Ap(x,n2)  -  o^q(.x-U2*^^)/2) 

Then  one  forms 

Ap(x,y)lQQ  +  °10^10 

The  integrals  must  be  formed  around  the  entire  ^n-element.  Only 
the  limits  due  to  the  points  2  and  3  will  contribute  to  the 
singularity.  Point  2  is  the  upper  limit  for  the  integration  along 
the  side  1,2  and  the  lower  limit  for  the  side  2,3.  Point  3  is  the 
upper  limit  for  the  integration  along  2,3  and  the  lower  limit  for 
the  integration  along  3i1.  One  obtains  the  following  contributions 


from 

iJot'/S) 

from 

b^y‘ 

from 

^00*“l2* 

B^Y- 

from 

-B^Y 

B  y'  G^Q2[x-(C2+C3)/2)]((x-e3)-(x-52)^ 


The  sum  of  these  expressions  cancels,  indeed. 

In  the  analytic  integration  over  x  and  y,  and  (j>jj  give 
separate  formulae.  The  Y  ^  singulatity  (which  appears  in  the 
individual  expression)  becomes  a  log  Y  singularity  after  the 
Integration.  If  Y  -  0  is  one  of  the  limits  of  integration,  then 
one  obtains  infinity.  It  is  true  that  these  infinities  will  always 
cancel.  This  happens  for  points  upstream  of  the  5n-element  because 
no  singularities  are  present,  for  wake  points  if  Ap  is  continuous 
because  of  the  contribution  of  an  adjacent  Cn-element;  for  wake 
points  if  Ap  is  discontinuous  because  the  region  of  integration 
will  then  extend  across  the  wake  and  the  Y~^  singularity  changes 


V*  ‘  •  V.  •*,  , 
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sign  as  one  crosses  the  wake.  In  any  case  the  formulae  are  needed 
for  the  numerical  procedure  which  show  the  singularity  separately. 
They  will  be  derived  later. 

For  elements  where  a  ^  ir/2  for  all  sides,  no  Y  ^ 
singularities  will  appear.  As  mentioned  above,  this  is  of  interest 
if  one  works  with  elements  in  which  Ap  is  constant  for  then  one  can 
choose  the  individual  elements  as  pressure  and  upwash  areas. 

Weaker  singularities  occur  in  and  (l)^.  In  =  log(k(R-X)) 
the  argument  of  logarithm  becomes  zero  for  Y  =  0  and  X  >  0.  To 
display  the  singularity  one  writes 

-  log(k(R-X)  -  logkY^(R  +  X)"^  -  log(kY^)  -  log(R+X) 

In  some  of  the  expressions  41^  is  multiplied  by  powers  of  Y,  which 
makes  the  singularity  less  pronounced.  Upon  integration  with 
respect  to  x  and  y  one  obtains  finite  quantities.  One  does  not 
depend  upon  cancellation  of  different  terms,  but  one  must  make 
sure,  that  the  formulae  for  the  integrals  do  not,  inadvertantly, 
give  the  difference  beween  two  infinite  quantities. 

For  I?,  (•it/2),  the  expression  41^^  does  not  introduce 
singularities  except  at  the  boundary  of  the  element.  This  is 
obvious  for  points  Y  -  0  if  x  <  ^  ^3*  x  >  ^2  >^3 

obtains  after  substitution  of  the  limits 

[Y"log(kY^)  -  Y"log(2(x-C3))]  -  [Y"log(kY^)  -  log( 2 ( x- ^2 ) )  ] 

The  logarithmic  terms  cancel. 

However,  the  same  argument  cannot  be  made  for  <  x  <  ^2 • 

The  expression  due  to  ^2  does  not  generate  a  singularity,  the 
singularity  of  the  expression  due  to  remains  uncompensated. 

2  3  2 

In  Iqq  (11/2),  the  function  41^  is  multiplied  by  X,  X  etc. 

Then  the  singularity  will  not  cancel  for  x  >  >  ^3 !  ® 

logarithmic  contribution  to  the  wake  singularity.  This  is  to  be 


expected  if  one  examines  the  portion  of  the  kernel  K;  I.,  arises 
from  the  lowest  order  term  in  k  that  contains  the  logarithm. 

In  the  expressions  for  a  /  n/2,  the  function  (|)^  is 
encountered  in  the  following  forms 

2 

S  tga<J)^ 

2  -2 
-&  sinacos  aU<i>^  , 

.2  -1  „ 

-S  cos  aU(j»^ 


in  Iqq  (a) 


•  t02  ,  , 

in  Q  (a) 


^  •  t02  ,  , 

and  in  Iq.  (a) 


Only  an  integration  limit  for  which  Y  -  0  will  contribute  to  the 
singular  part.  In  general,  there  is  no  second  compensatory  term. 

These  are  the  terms  responsible  for  the  wake  singularities  for 
elements  sketched  in  Fig.  3b. 

The  terms  (^>2  =  log((R-V)/U)  occur  in  integrations  along  lines 
U  -  const.  They  occur  multiplied  by  different  powers  of  U.  Terms 
log  U  can  be  considered  as  constants  of  integration  which  cancel 
(even  for  n  -  0  and  in  the  limit  U+O).  Therefore,  it  suffices  if 
we  discuss  expression  log  (R-V).  One  best  returns  to  the 
original  coordinates  u  =  -U,  v  =  -V.  For  a  point  to  the  line 
U  =  0,  one  must  evaluate  (-U)'^[log(R2  +  V2  )-log(R^ +v^  )  ].  If  the 
point  (x,y)  lies  close  to  the  line  u  =  const  and  v  >  v^  >  v^,  then 
R2  +  V2  =*  2R2;  *  2R^  the  terras  are  not  singular.  For  v  <  v^ 

one  writes  log(R^+v^)  -  log(  — ),  if  v  <  V2  the  singular  term 
2  ^ 

log  u  is  canceled  by  the  contribution  of  point  2.  For  v^ <  v  <  V2, 
this  cancellation  does  not  occur.  The  terms  (^>2  express 
singularities  in  the  upwash  along  element  sides  for  which  a  *  Tr/2. 

According  to  these  considerations  one  must  examine  how  the 
formulae  for  the  integrals  derived  from  4)^  (combined  into  a  vector 
)  appear  in  the  limit  Y+O,  and  how  the  formulae  for  4)2  appear  in 
the  limit  U  »  0.  One  should  always  obtain  finite  expressions. 

Next,  consider  4)^  in  the  case  a  ^  Tr/2,  6  /  n/2.  The  origin  W[^ 

lies  at  a  point  t  ■  n  -  :  one  integrates  along  a  line  p  --  const  from  q. 


to  q  Singularities  in  the  vector  arise  if  one  of  the 

endpoints  of  the  interval  of  integration  lies  at  Y  -  0.  The 
components  of  the  vector  \i)^  are  the  integrals  shown  in  Eq.  (I8l). 
The  term  which  becomes  singular  for  Y  -  0  is  log(l<(R-X),  it  behaves 
as  log  Y).  The  factor  q"  -  (ptgY)”  of  this  term  vanishes  as 
(Ycos'^Y).  The  term  is  simply  disregarded  if  Y  -  0  at  one  of  the 
limits.  For  Y  -  ir/2,  one  has  p  -  -Y.  The  expressions  are  then 
found  in  Eqs.  (18M).  They  do  not  vanish.  But  the  element  of  the 
diagonal  matrix  with  powers  of  p  »  -Y  now  vanish.  Again  one  simply 
disregards  this  limit. 

Similar  considerations  apply  to  therefore  4<2)* 

6  y  0  the  singular  element  log  (R~V)/U  in  Eqs.  (182)  vanishes  if  V 
-  0  at  one  of  the  limits.  For  6  0,  p  “  u,  and  the  diagonal 

matrix  with  elements  given  by  powers  of  p  vanishes. 

Now  we  identify  the  singularities  due  to  =  Y  and 
.  RY~^  which  appear  after  the  integration  with  respect  to  x  and 
y  has  been  carried  out. 

We  begin  with  cases  a  -  •tr/2. 

Then  one  has  and  6^R<J>4  as  contributions  to  lQQ(ir/2); 

and  (3B^/2)X^(P2  and  (3B^ /2 )XR(|.i,  as  contributions  to  I^q(tt/2). 

Discussion  for  <^2* 

Case  Y  -  Tr/2. 

One  remembers  that  the  areas  of  integration  over  x  and  y  are 
triangles.  Therefore  one  forms 


■jj  X(p^(lx(iy  -  XY~^dxdy 


-1  j  Xdx  -  -(1/^)X‘ 
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The  negative  sign  arises  because  the  convention  that  one  integrates 
around  the  xy-eleraent  in  the  counterclockwise  sense.  If  the 
integration  extends  from  point  j  to  point  j+1  in  Fig.  14  then  the 
element  lies  above  the  line  j,  j+1  and  the  contribution  of  the 
triangle  must  be  subtracted.  In  the  notation  of  Eq.  (178)  one  then 
has 

(ir/2,  ir/2)  -  -(e^/4)X^ 

Similarly 

ir/2)  -  -(e^/6)X^ 

In  spite  of  the  fact  that  the  area  of  the  triangle  tends  to  zero, 
the  integral  over  the  triangle  is  finite.  The  same  behavior  will 
be  found  for  J  |,  ( ir/2 ,  ■ir/2 ) . 

For  and  Y  A  ir/2  one  has 

I  XY"’dq 

-  0  (p/2) J(qsin6  +  pcos6)(qcos6  -  psin6)  dq 

Hence , 

Jqq  ^(ir/2,Y)  -  (0^/2)[tgYpq  +  p^cos  ^Ylog(k|Yl)] 

If  one  of  the  limits  of  integration  should  be  Y  -  0,  then  the  last 

2  -2 

term  gives  a  logarithmic  singularity.  For  y  -  0,  one  has  p  cos  Y 
2 

-  X  +  0(y).  Thus  one  can  write 


Jqq^2(ii/2,Y)  -  (6^/2)[tgYpq  +  (X^ 


0(Y)log(k|Yi )) 


This  displays  the  singularity  in  terras  of  X  and  Y.  Of  course,  in 
all  applications  these  infinities  raust  ultiraately  cancel.  The  only 
purpose  to  retain  them  in  the  coraputations  as  separate  terras  is  to 
provide  a  check. 

For  the  contribution  of  to  J^qCtt/Z.Y)  one  raust  evaluate 

-  (36^/2)(p/3)|  X^y'^dq 

-  (36^/2)  (p/3)  J  (^^(sin^Yq^  +  2sinYcosYpq  +  cos^Yp^)dq 

The  individual  integrals  are  found  in  Eq.  (183).  Singularities 
arise  because  of  the  factor  log  (kY)  if  one  of  the  limits  is  Y  -  0. 
Collecting  the  terms  contributing  to  the  singularity  one  obtains 

(36^/2)  (p^/3)G0s~S[3in^Ytg^Y  2sinYcosYtgY  +  cos^Y) log(ky ) 

-  (36^/2)  (p^/3)cosY(tg^Y  +  2tg^Y-*-1)  -  (3B^/2)  (1 /3)  (p/cosY)^log(kY) 

Thus,  substituting  the  remaining  terras  of  Bqs.  (183) 

J?o_3(./2.’r)  . 

( e^/2) {(sinYtg( Y/2)pq^  +  p^qsinY(2  +  tg^Y)  +  (X^  +  0(Y) )log(k/Y) } 

This  terminates  the  evaluation  of  I  ,  for  a  =  ir/2.  (The  results 
can  be  summarized  by  stating  that  it  suffices  that  one  disregards 
the  singular  terras,  provided  of  course,  that  one  uses  analytical 
expressions  in  which  they  cancel  on  theoretical  grounds.) 

Now  we  discuss  the  corresponding  terms  for  for  a  »  tt/2 


Jqq  ^^(^r/2,^r/2)  -  -6^/2)  |  Rdx  =  -  (  B^/4 )  [  RX- Y^l  og  ( k  ( R-X )  ] 

For  Y  -»  0,  X  >  0,  we  write  logk(R-X)  »  logkY^-log(R+X) .  For  Y  -  0, 
the  expression  remains  finite  (although  the  area  of  the  triangle  is 
zero)  one  obtains  for  Y  -  tt/2 


'00,4 

,0 


(it/2, ir/2)  -  -(B  /'DRX  for  Y  -  0, 


j'j'Q^i,(ir/2)  ,it/2)  -  -(3B  /2)(l/3)  J 

J°Q^^(Tr/2.ir/2)  -  -  (B^/6)R^ 


XRdX 


For  Y  /  11/2,  one  has  to  evaluate 


•'OO.K  ■  IV” 


This  Is  one  of  the  expression  Eq.  (183).  The  individual  terms  on 
the  right  of  Eq.  (183)  are  found  in  Eq.  (180).  After  substitution 
one  obtains 


jJqj^(it/2.Y)  - 

-  (6^/2)  {p^cos"^Ylog((R-X)/Y)  +  cos“S[pR  -  p^tgYlogk(R-q)  ]} 

-  (B^/2)  {((X^+0(Y))(log(R-X)-logkY)+cos"S[pR-p^tgYlog(k(R-q)]} 


( 6^/2) {(X^+0(Y) )(log(kY)-log(k(R+X)+COs“’ Y[pR-p^tgYlog(k(R-q) ]} 


The  first  version  is  the  general  formula,  the  second  version 

dislays  the  singularity  for  X  <  0,  and  the  third  version  the 

singularity  of  X  >  0.  The  last  term  is  singular  for  p  -  0,  but 

2 

then  it  vanishes  because  of  the  factor  p  .  Next  we  evaluate 

Next  we  evaluate 

J°q^j^(tt/2,Y)  -  (36^/2)(p/3)  I  X(1)^dq 

•  ( B^/2)  [cosYp^  f  (j)i,dq  +  sinYp  |  q^^^dq] 


Here  expressions  from  Eq.  (I83)  are  substituted: 

-  (g^/2){p^cos  j  R  N  ^dy  ♦  ptgY  j  R  \^dq  +  p^cos  ^Y  |  R  \dq 
+  p^tgYd  +  cos  ^Y)J  R  ^dq} 

Only  the  first  term  within  the  braces  has  a  denominator  Y,  which 
will  give  a  singularity  for  Y  -  0.  Substituting  expressions  from 
Eq.  (180)  one  obtains 

-  R^/2{p^cos~^Ylog(R-X)/Y)  +  (ptgY/2)  (Rq+p^log(k(R-q) ) 
+  p^cos  ^YR  -  p^tgYd  +  cos  ^Y)log(kR-q) ) ) 

The  singularity  for  Y  -  0  is  displayed  by  writing 

p^cos'^Y  -  ♦  0(Y) 

Then 

p^cos"^Ylog(R-X/Y)  -  (X^  +  0(Y)(log(k(R-X)  -  log(kY)) 

-  (X^  0(Y)(-log(k(R+X)  +  log(k/Y) 

The  first  form  is  suitable  if  X  <  0,  the  second  if  X  >  0. 

For  a  ^  7r/2  one  finds  terms  with  and  in  I^q.  namely 

2  -2  2  2  -2  -1  2  -1  2  -1  -1 
B  cos  a  U  (l)^  -  B  cos  a  U  Y  and  B  cos  a  U(j>j^  -  B  cos  a  URY 


For  Y  -  n/2,  one  has 


u 


-Y  3ino  +  X  cosa 


Y^sin^a  -  2XY  sina  coso  ♦  X^cos^a 


Hence 

J°q  -  B^coa  ^a(p/3)  | 

-  (B^/3)co3  ^o[-3in^oY^  j  dx  *  23inoco3oYj  Xdx  -  coa^oj  X^dx 

J°0^3(«.if/2)  -  (B^/3)C-tg^aY^X  ♦  tgaYX^  -  (co3^o/3)X^] 

Here  no  aingularity  for  Y  -  0  occura;  the  expreaaion  doea  not 
vaniah,  even  though  the  area  of  the  triangle  from  which  it  ariaea 
ia  zero. 

Next  we  diacuaa 

J^Q  ij(a,Tr/2)  -  (6^/3)p  coa  U<>^dq 

-  (B^/3)co3  ^aCYaina  jRdX-coaa  jRXdx] 

-  (6^/3)  {((tgo)/2)Y[RX  -  Ylog(l<(R-X)  ]-  j  R^} 

The  aingularity  which  ia  preaent  for  X  >  0,  ia  of  the  character  Y 
log  Y.  For  Y  -  0  one  obtaina 

J°  i.(o,Tr/2)  -  -(S^/9)R^ 


For  Y  ^  ir/2  one  has 


J°q(o.Y)  -  6^(p/3)cos~^a  |U^(|»3dq 

•  6‘^(p/3)cos  a  {cos  6  p  J  ♦^dq  •*■  2sin6cos6  p  J  ♦^qdq 

2  2 
+  sin  6  I  ^^q  dq} 

Here  expressions  from  Eq.  (183)  are  substituted.  We  collect  terms 
of  the  same  character;  the  terms  logk(Y)  is  of  particular  interest 

(e^/3)cos  ^acos  ^ Y{(2sin6co36  ♦  3in^6tg6)p^q  +  ( (sin^6)/2)pq^ 

+  p^(co3^6  +  2sin6co36tgY  ♦  3in^6tg^Y)log(kY) } 

The  factor  of  log  kY  simplifies  to 

p^cos^fid  +  2tg6tgY  +  tgfi^tg^Y)  -  C0S^6(1  ♦  tg6tgY)^ 

-  p^cos  ^Ycos(Y-6)^ 

But 

Y  -  6  -  a 

Thus  one  finds 

J^Q  3(0, Y)  -  (0^/3) {cos  ^ocos  N[sin6C036(2  +  tg6tg6)p^q 
+  ( 1 /2)sin^6pq^]  +  p^cos  ^Ylog(kY)) 

One  has  p^cos  ^Y  -  +  0(Y). 

This  shows  the  character  of  the  singular  terra;  for  Y  ♦  0  it  behaves 
as  (e^/3)X^log(kY). 


a 


Consider  finally 

J^O  ij(Qit'Y)  -  0^003  ^a(p/3)  I  U^jjdq 

O  —  1 

-  0  cos  a(p/3)[cos6  p  I  ♦ijdq  +  sin6  j  q<j»^dq] 

Here  the  expressions  from  Eq.  (183)  are  substituted.  One  then 
obtains 

J°0  ij(oi,Y)  -  (0^/3)coa  ^  a{p^co3~^Y(cos(  6-Y)  j  R  V  ^dq 

2  ”“2  ~  1 
+  p  cos  Ycos(6-a)  j  qR  dq 

♦  p^cos  ^Ycos6(3inYcosY  ♦  1 )  J  R  dq 
-K. _ _  f  2^., 


cos"  Y  3in6  p  I  r"  q  dq} 


Here  expressione  from  Eq.  (180)  and  6-Y  -  o  are  substituted 
J?0,1»(“’Y)  -  (0^/3)  {p^coa"^Ylog{R-X)/Y)] 

+  p^cos  ^YR-p^cos  ^ocos  ^Ycos6(sinYcosY-*-1  )log(l<(R-q) ) 
-1-1  3 

+  (1/2)co3  ocos  Ysin6[Rpq  +  p  log(k(R-q) ) ]} 

3-3  3 

Here  again  p  cos  Y  -  X  +  0(Y),  so  that  the  singular  term  due  to 
log(R-X/Y)  is  readily  recognized. 


SECTION  VIII 
CONCLUDING  SURVEY 


The  sheer  volume  of  details  and  formulae  may  make  it 
difficult  to  extract  from  this  report  the  information  needed  for 
numerical  work.  We  describe  here  the  main  concepts,  the 
definitions  and  list  the  important  equations. 

One  observes  that  the  formulae  become  longer  (although  not 
really  more  complex)  as  one  takes  higher  terms  in  the  development 
with  respect  to  the  reduced  frequency  into  account.  But  this 
applies  only  if  the  xy-element  coincides  with,  or  lies  in,  the 
vicinity  of  the  ^n-element.  For  elements  at  a  distance  the 
different  terms  of  the  development  with  respect  to  frequency  are 
lumped  together,  one  obtains  rather  smooth  functions,  the 
Integrations  are  carried  out  numerically  and  this  difficulty  does 
not  arise.  Ir  any  case  one  must  choose  a  subdivision  of  the  wing 
into  elements  so  that  the  linear  elemental  pressure  distribution 
give  a  sufficient  accuracy.  This  implies  that  for  an  xy-element  in 
the  neighborhood  of  an  ^n-element  the  terms  of  low  order  in  the 
development  with  respect  to  the  reduced  frequency  will  suffice. 

Most  likely  the  powers  zero  and  one  will  be  enough.  (For  elements 
at  a  greater  distance  one  should  go  further  in  the  development. 

But  we  mentioned  already  that  then  the  integrations  can  be  carried 
out  without  difficulty.) 

The  overall  arrangement  of  the  computation  is  symbolized  by 

Eq.  (7).  The  definitions  of  the  "housekeeping"  matrices  and 

( 3 ) 

M  are  given  in  the  paragraph  preceding  Eq.  (7).  An  element  of 
the  matrix  gives  the  average  upwash  in  some  xy-elsment  due  Vo 

one  elemental  pressure  distribution.  (In  each  triangular  Cn- 
element  there  are  three  elemental  pressure  distributions.) 

The  elemental  pressure  distributions  are  expressed  in  terms 
of  coordinates  5  and  n;  the  coordinates  within  an  upwash  element 
are  x  and  y.  The  pressure  difference  between  the  upper  and  lower 


side  of  the  wing  is  denoted  by  ApC^.n).  to  make  analytical 
integrations  possible  one  must  consider  Ap(^,n)  •  Ap( C. n)exp( ik^) 
as  unknown  functions  (Eq.  (6)). 

To  each  triangular  element  there  belong  three  elemental  shape 
functions,  shown  in  Eq.  (11).  They  are  linear  functions  which 
assume  the  value  one  at  the  corner  of  the  triangle  with  subscript  j 
and  zero  at  the  other  corners. 

For  an  xy-element  at  a  sufficient  distance  from  the 
element  one  substitutes  the  expressions  Eq.  (11)  into  Eq.  (10)  and 
carries  out  the  integrations  numerically.  To  make  analytical 
integrations  feasible  Eq.  (10)  contains  a  weight  factor  exp(ikx). 

For  xy-elements  at  a  distance  from  the  ^n-element  but  within 
its  wake  or  close  to  it,  one  applies  the  procedure  given  by 
Eqs.  (20)  through  (31)  (if  one  restricts  oneself  to  the  powers  zero 
and  one  of  the  reduced  frequency  and  by  Eqs.  (20)  through  (42)  for 
the  general  case).  In  the  equations  (without  number)  starting 
after  Eq.  (42)  and  proceeding  to  the  end  of  the  section  the 
procedure  is  specialized  to  triangular  elements  and  carried  out 
analytically.  One  stops  with  the  evaluation  of  f^  if  the  powers  of 
the  reduced  frequency  are  only  zero  and  one. 

Before  the  integrals  of  the  upwash  over  an  xy-element  close 
to  a  En-element  are  evaluated,  we  have  derived  formulae  for  the 
upwash  at  a  given  fixed  point  x,y.  These  results  are  of  an 
intermediate  nature.  There  is  no  need  to  evaluate  these  formulae 
unless  one  wants  to  show  the  details  of  the  upwash  distribution  due 
to  one  elemental  pressure  distribution.  We  describe  these  results 
because  they  form  an  important  part  of  the  overall  procedure.  The 
pressure  difference  Ap  are  written  in  the  form  Eq.  (5i).  The 
functions  C[njj(x,y)  are  expressed  in  terms  of  the  elemental  pressure 
distribution  given  by  Eq .  (11).  (In  practice  only  the  combineti ens 
(0,0),  (1,0),  and  (0,1)  for  (m,n)  are  used.) 

If  one  chooses  to  work  only  with  (m,n)  -  (0,0),  then  one  ,/u;, 
elements  with  constant  Ap ,  and  if  it  is  not  possible  to  maintain 
continuity  as  one  moves  from  one  element  to  its  neighbor.  Then  if 


the  5n-element  has  a  side  parallel  to  the  ^  axis,  one  has  at  the 
wake  line  originating  at  this  side  (n-y)  ^  singularity  in  the 
upwash.  To  obtain  finite  averages  for  the  upwash,  the  upwash  areas 
(over  which  the  averages  are  formed)  must  overlap  these  wake  lines. 
If  none  of  the  sides  of  the  5n-element  is  parallel  to  the  5  axis, 
then  the  singularities  are  of  the  character  log(n-y)  and  the 
average  remains  finite  without  an  overlap.  One  can  then  develop  a 
method  in  which  the  pressure  and  upwash  areas  coincide  with  the 
pressure  elements. 

The  general  procedure  obtained  after  a  limiting  process  z-»0 
has  been  carried  out,  is  summarized  in  the  remarks  following 
Eq.  (88a).  The  case  m  -  0,  n  -  0  is  exceptional,  the  procedure  is 
given  by  Eq.  (88).  To  obtain  the  upwash  at  a  fixed  point  xy,  one 
has  to  carry  out  integrations  around  the  contour  of  the  ^n-element, 
separately  for  each  individual  term  in  the  representation  for  the 
pressure . 

Carrying  out  the  procedure  one  obtains  definite  integral  over 

the  aides  of  the  triangle,  denoted  by  The  subscripts 

ra,n 

(here  m  and  n)  refer  to  the  term  in  the  development  of  Ap,  the 
superscripts  on  the  left  (here  i  and  i+1 )  refer  to  the  numbering  of 
the  corners  of  the  5n~eleraent.  The  indefinite  integral  is  denoted 

0 

by  The  slope  of  the  side  enters  by  Eq.  (107) 


The  indefinite  integral  depends  upon 


Furthermore,  we  have  used 


R  -  (X^ 

U  »  X  cosa  -  Y  sina 
V  -  X  sina  -  Y  cosa 


(see  Eq.  (140)). 

(  o') 

The  indefinite  integrals  are  denoted  by  ^  (X,Y,a).  Here  the 
limits  of  5  and  n  for  the  individual  sides  of  the  triangle  must  be 
substituted.  One  obtains  a  number  of  functions,  that  depend  upon 

I 

t  h~yj  t  and  a.  The  expression  I  „  ^  are  divided  because 

1  1  1,1+1  m,n 

of  the  different  analytical  characters  of  the  expressions,  for 
instance 


I 


02 

00 


The  results  are  found  in  Eqs.  (114)  through  (120)  for  a  -  Tr/2,  and 
Eqs.  (141),  (147)  through  (155). 

Formulae  for  the  integration  with  respect  to  x  and  y  are 
derived  in  Section  VI.  The  integrals  appear  now  as  contour 
integrals  around  the  xy-element.  A  summary  of  the  results  before 
this  integration  is  found  in  Eqs.  (167),  (16«),  and  (169),  they  are 
a  repetition  in  a  more  condensed  form  of  the  equations  listed 
above,  except  for  the  introduction  of  I  ^  and  ,  defined  by 
Eq.  (157). 

In  the  integration  the  slope  of  the  sides  of  the  xy  element 
appears  given  by  Eq.  (161) 


tga 


i ,  i  +  1 


1 1  4 


V 


kb  ' 


The  angle  a  defined  for  the  aides  of  the  5n-element  appears  again. 
We  have  furthermore  (Eq.  (I64a)) 

6  »  Y  -  a 

Regarding  this  last  phase  of  the  computation  we  refer  to  the 
description  given  in  Section  VI. 

In  evaluating  these  expressions  one  will  find  infinities. 
These  and  other  singularities  are  discussed  in  Section  VII. 
Critical,  of  course,  are  the  lines  of  the  wake  that  start  at  the 
corners  of  the  ^n-elementa ,  but  also,  although  the  singularities 
are  less  pronounced,  the  aides  of  the  ^n-triangle  and  their 
straight  extensions  through  the  5n-plane.  For  methods  in  which  the 
effect  of  such  infinities  cancel  on  theoretical  grounds,  they  can 
simply  be  disregarded,  when  they  appear  during  the  computations. 
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APPENDIX  A 
THE  BASIC  EQUATION 


The  classical  equation  for  oscillatory  flows  is  rederived 
here  as  a  convenience  to  the  nonspecialist  who  would  like  to 
know  how  the  equation  comes  about,  but  does  not  want  to  go  back 
to  the  original  literature. 

The  formulation  uses  the  accleration  potential.  Usually, 
and  also  in  the  present  case,  the  acceleration  potential  is  used 
in  connection  with  the  linearized  flow  equations.  But  according 
to  an  observation  (which  was  made  to  the  author  by  Ernst  Hoelder) 
it  has  also  a  meaning  in  a  nonllnearized  isentropic  flow;  it  is 
given  by  the  negative  enthalpy  of  the  particle. 

The  following  derivation  is  made  for  the  linearized  flow 
equation.  Denoting  by  p,  p,  and  4»,  respectively,  the  deviation 
of  the  pressure  from  the  free  stream  pressure,  the  free  stream 
density,  and  the  acceleration  potential,  one  has 

p  =  -pif'  (A.l) 

In  a  linearized  flow  with  free  stream  velocity  U,  and  velocity 
components  U  +  u,  v,  and  w  in  the  x,  y,  and  z  directions*, 
respectively,  of  a  Cartesian  system  of  coordinates,  one  has  the 
components  of  the  acceleration 

Ou/3t)  +  U(3u/9x), 

(3v/3t)  +  U(3v/3x), 

and 

(3w/3t)  +  U(3w/3x). 

Expressed  in  terms  of  a  velocity  potential,  4),  the  velocity 
components  are 


for  simplicity  we  use  in  this  Appendix  x,  y,  z,  etc.,  before  a 
Prandtl-Glauert  coordinate  transformation  has  been  carried  out 
although  in  the  main  body  of  the  report  the  notation  x,  y,  z  is  used 


u  =  0(j>/3x),  V  =  3(j)/3y,  and  w  =  3({i/3x 


Therefore , 


*  =  If  * 


The  velocity  potential  is  then  expressed  in  terms  of  an  acceleration 
potential  \^(x,y,z,t) 

X 

(l)(x,y,z,t)  =  U“^  I  \|»(5,y,z,t  + 


(The  velocity  potential  will  be  used  in  an  intermediate  step  to 
evaluate  the  upwash  due  to  a  given  acceleration  potential.)  One 
has  Indeed 


A 

U(3(j)/3x)  =  )i)(x,y,z,t)  -  U“^  |  i(>^(5,y,z,t 


=  ij<(x,y,z,t)  -  3(ji/3t 


Choosing 


<t)(x,y,z,t)  *  (t(x,y,z)exp(ivt) 

«v 

i|j(x,y,z,t)  =  \j^(x,y,z)exp(ivt) 


One  obtains 


(|)(x,y,z)  =  U 


I  V>(C,y, 


z)exp(iv^^)dC 


(j)  and  ip  are  made  dimensionless  with  UL  and  U  ,  respectively,  where 
L  is  some  characteristic  length. 

0  «  UL| 
a  u^ijJ 


Furthermore 


X  =  x/L,  y  =  y/L,  z  =  z/L,  and  ~ 


•*   •*  a'".  «*  .•* • 


Then  one  has 


<J)(xL,yL,zL)  =  UL|(x,y,z) 
t|»(xL,yL,zL)  *  U^ii5(x,y,z) 

and  one  obtains 

X 

♦  (x,y,z)  =  I  ii/(^,y,z)exp(lk(C-x)  )d^ 
One  has,  from  Eq.  (A.l) 

p  =  -pU^ijj 


(A. 2: 


(A. 3 


Setting 


p  =  Ap/(pU^/2) 


(A. 4 


one  obtains 


p  «  -2ij; 


(A. 5 


From  now  on  the  bars  will  be  omitted. 

4)  and  ip  satisfy  the  same  linearized  differential  equation. 
Originally, 


(1  -  +i|)  +1^  =  0 

“^^xx  ^yy  ^zz  2  ^xt  „2  ’^tt 

a  a 


(A. 6 


(where  a  Is  the  free  stream  velocity  of  sound  and  M  =  U/a  Is  the 
Mach  number).  Introducing  nondimensional  quantities  one  obtains 
for  the  oscillatory  case  (after  omission  of  the  bars) 

(1  -  M^)ip  +41  +41  -  21k  4j  +  k^M^4j  =  0  (A. 7 


Let 


2  2 

1  _  M  =  B 

2  2^2 

r  =  y  +  z 

d2  2  -2  2 

R  =  X  +  $  r 


A  solution  of  Eq.  (A. 7)  for  outgoing  waves  which  has  an  R" 
singularity  at  the  origin  is  given  by 


(A. 8) 


source 


=  -  R~^  exp^^^^  (Mx  -  R)j 


(A. 9) 


Eq.  (A. 9)  can  be  verified  by  substitution  into  Eq.  (A. 7).  If 
the  expression  would  be  used  for  the  potential  it  would  represent 
an  oscillating  source.  The  expression  for  to  be  used  here  is 
given  by 


U)  =  /3z 

^  ^source 


One  has 


The  potential  pertaining  to  it  is  obtained  from  Eq .  (A. 3) 


(A. 10 


4>(x,y,z)  -  -z  I  I  2  2  2  ‘^/2  ^  2  2  2  1/2 

expjp^  (MC  -  +  lk(5-x)jdC 

or  using  the  first  of  Eqs.  (A. 8) 


(()(x,y,z)  =  -z  exp(-lkx) 


(E^  +  (E^ 


(A.il 


exp[^[5  -  M(c^  +  B^r^)^^^]jd^ 


To  simplify  the  argument  of  the  exponential  function,  one  Introduces 


and  (for  the  upper  limit  of  the  integral) 


V(x,r)  =  6"^[x  -  M(x^  + 


(A. 12) 


(A. 13) 


(In  Ref.  1 ,  X  is  used  instead  of  V.  The  author  has  changed  the 
notation  because  the  variable  V  does  not  solely  depend  upon  x.) 

Then,  from  Eq.  (A. 12) 


and  from  Eq .  (A. 13) 


C  =  V  +  M(r^  + 


X  =  V  +  M(r^  + 


(A.li4) 


Hence 


d?  =  dv 


(v^  +  +  Mv 

(r^  +  v2)1/2 


R(C,r)  =  (C^  +  B^r^)^/^  =  (v^  +  r^)^^^  +  Mv  (A. 15) 


analogously 


t(x,r)  =  (x^  +  B^r^)^^^  =  (V^  +  r^)^^^  + 


(A. 15a) 


Then  by  substitution  into  Eq .  (A. 11) 


<l>(x,y,z)  = 


r 

-z  exp(-ikx)  exp(lkv)  - p - r  ^ - 5 — 5 - — ^ 

.i  +  Mv]^(r^  +  v^)^' 

.  _ IkM  _  1  , 

C(v^  +  +  Mv](r^  +  j 


In  the  second  term  of  the  bracket  in  the  integrand  an  integration 
by  parts  is  carried  out; 


())(x,y,z)  =  -z  exp(-ikx)  - p - g  ?/2 — ^ 

[(v'^  +  +  Mv](v^  + 

V 

f  1 

+  exp(ikv)dv  - 2 - 2  1/2  ? — ? - TTTJ 


[(v^  +  +  Mv](v^  + 


The  contribution  of  the  lower  limit  in  the  term  outside  of  the 
integral  vanishes.  With  Eq .  (A. 15)  one  obtains  for  the  upper 
limit 

M  exp(i.kV) 

R(x,r)(V^  + 

The  integral  simplifies  to 


[  exp(ikv) 

J  (v^  +  r^) 


Thus,  one  finally  finds  for  the  velocity  potential  that  pertains 
to  the  acceleration  potential  given  by  Eq.  (A. 10) 


i(x,y,z)  =  -z  exp(-lkx)  — 


R(x,r)(V^  + 


I  _|2EU|v^dv 

i  (v^  +  r^)3/2  J 


(A. 16) 


Where  V  is  defined  in  Eq .  (A. 13),  r  and  R  in  Eq.  (A. 8).  The  upwash 
is  then  given  by  3if>/3z.  Eq.  (A. 10)  is  a  fundamental  solution  with 
the  singular  point  at  the  origin.  The  flow  field  is  represented 
by  a  superposition  of  such  fundamental  solutions  but  with  singular 
points  (?,n)  lying  in  the  plane  of  the  wing.  One  then  must  replace 
x,  by  x-C,  and  y  by  y-n .  Consequently,  one  has  now 


(A. 17 


R 

V 


[(y-n)  +  z  ] 

[(x-C)^  + 
e"^[(x-C)  -  MR] 


We  note  that 

(V^  +  =  e“^(R  -  M(x-C)) 


) 


(A. 18 


Let  f(C,n)  be  the  strength  of  the  doublets  assigned  to  the 
particular  solution  Eq.  (A.  10).  Then 

iKx.y.z)  =  -z  I  dCdnf(C,n)J^  +  exp[p55(M(x-C)  -  R)] 


Here  A  Is  the  wing  area.  We  determine  lln  (i|;(x,y,z).  The  limit 
is  obviously  zero,  as  long  as  R  ^  0,  for'*^?nstance  for  all  points 
outside  of  the  wing.  Assuming  that  the  origin  lies  within  the 
wing  area  we  evaluate  (for  simplicity)  the  expression  for  x  =  0, 
y  =  0,  z  =  e  >  0),  l.e.,  we  approach  the  plane  of  the  wing  from 
above.  Accordingly  we  consider 


,0,e)=  -e  I  dCdnf(C,n) 


A 
Mik 

cc"  -H  t  c2)1F^ 


exp! 


rikM  (  I  f.2  ,  „2 ,  2  , 

(-MC  -  (C  +  6  (n  +  G  ) ) 


(A. 20 

2._2  ,  2, a/2 


We  cut  out  from  the  wing  area  a  small  ellipse,  given  by 

C  =  Y  cos  a,  n  =  Y  8”^  sin  a 


where  y  is  constant  and  Independent  of  g.  For  points  outside  of 
this  ellipse,  the  integrand  is  bounded;  because  of  the  factor 
e  in  front  the  expression  vanishes  for  e  =  0.  Retaining  within 
the  ellipse  only  the  dominant  terms,  one  finds 


v' 
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The  acceleration  potential  Eq .  (A. 20),  therefore,  gives  at  point 
(C,n)  the  value  -2ir  f(5,n)  and  for  the  nondimensional  pressure 
defined  in  Eq .  (A. 5)  4irf(C»n).  The  contribution  to  the  lift  is 
the  pressure  difference  between  the  lower  and  the  upper  side.  It 
is  denoted  by  Ap .  Thus  one  finds 

Ap  =  -  8Trf(£;,n) 


The  upwash  due  to  this  lift  is  then  found  from  a  corresponding 
superposition  of  expressions  Eq.  (A.l6) 


w(x,y,z) 


“  ^  ^  j 


Ap(c,n)K(x-c,  y-n,  z,k)dCdn 


(A. 21) 


where 


K  =  exp(-ik(x-C) 


V 


M  exp(lkV) 

R(x-?,r)(V(x-?,r)^  +  r^)'^'''^ 


f  exp(lkv) 

I  /„2  .  „2,3/2 


(v  +  r  )■ 


(A. 22) 


r,  R  and  V,  and  (v^  +  r^)  are  given  inEqs.  (A. 17)  and  (A.l8). 
This  equation  is  due  to  Kuessner  it  holds  throughout  the  flow 
field.  Identifying  w(x,y,0)  with  the  upwash  found  from  the 
boundary  conditions  at  the  wing  surface,  z  =  0,  one  obtains  an 
Integral  equation  for  Ap(C,Ti)-  Frequently  this  is  written  as 

w(x,y,0)  =  ^  I  Ap(5,n)  K(E-x) (n-y ,k)dCdn 


l.e.,  one  makes  the  limit  z-*-0  immediately.  Because  of  the 
singularities  of  K,  this  expression  must  then  be  Inteipreted  by 
some  special  technique  (e.g.,  the  one  due  to  Mangier),  but  then 
one  has  to  verify  that  the  assumptions  made  by  Mangier  are 
applicable.  The  present  analysis  is  based  on  Eqs .  (A. 21)  and 


(A. 22).  The  evaluation  of  the  kernel  K  as  a  development  with 
respect  to  k  Is  due  to  Ueda  (Ref.  2).  (See  Appendix  B.) 


APPENDIX  B 

REDERIVATION  OF  SOME  FORMULAE  DUE  TO  UEDA 


The  derivation  of  the  formulae  on  which  the  present  report  is 
based  seem  to  be  available  only  in  Japanese  (Ref.  2).  At  least  to 
some  readers,  a  rederivation  may,  therefore,  be  desirable.  The 
results  shown  here  are  somewhat  more  detailed  than  those  of  Ueda. 
The  following  expression  is  to  be  evaluated: 

V 

B(k.f,v)  -  av  (B.i) 


A  slight  simplification  is  obtained  by  setting 


k  -  kr 

V  -  V/r 

V  -  ru 


Then  one  obtains 


with 


B(k,r,V) 


r"^  B(k,V) 


B(k,V) 


(u 


exp( iku) 

2  .  1)3/2 


du 


(B.2) 


(B.3) 


(B.ii) 


In  essence  the  formulae  are  obtained  by  a  development  of  the 
exponential  function  in  the  integrand.  But,  as  the  development 
generates  higher  and  higher  powers  of  u,  the  integrals  will  not 
converge  as  the  lower  limit  tends  to  negative  infinity.  Therefore, 
we  write 


B(k,V)  -  B, (k)  +  B„(k,V) 


(B.5) 


with 


B,(k)  .  I  du 


B  (k.V)  -  f 

^  J  (u^  ^  1) 


exp(lku) _  ^ 

''  3/2 


(B.6) 


(B.7) 


A  development  of  the  exponential  function  in  B2  leads  to  converging 
integrals.  Real  and  imaginary  parts  must  be  treated  separately. 


B^sCk.V)  -  (-)"  k^"  U2„(V) 


(B.8) 


Bp.Ck.V)  -  I  (-)"  + 

n-0  ' 


with 


1  V  ,.n 
“  hr  \  — 5“^ - TFp 


du 


(B.10) 


These  integral  can  be  expressed  by  elementary  functions.  Let 
be  the  pertinent  indefinite  integrals  including  the  factor 
(n! )  ^ .  Then 


Wq(u)  -  u(u^  ♦  1  )"’ 


W, (u)  -  -(u^  +  1  )"^ 


(B.11  ) 
(B.12) 


W2(u)  -  (1/2)[log((u^  +  +  u)  -  u(u^  +  1)“^^^]  (B.13) 
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sr 
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Further  expressions  can  be  obtained  recursively.  The  relation 


(u'^  +  1) 


-  (m-2)- 


(u^  -  D- 


tn-2 

+  - , 

(u^  D- 


(in-2)m!  (dW„/du)  +  ( m-1  )  (m-2 ) !  ( dW  „/du) 
m  in-2 


yields  the  recurrence  relation 


w„(u) 

m 


! m-2 )m!  ^^2  ^  1)^^^  m(m- 


m-2 


(B.14) 


Ueda's  formulae  are  based  on  this  relation. 

A  more  direct  formulation,  which  will  be  used  in  this  report 
is  based  on  formulae  in  Appendix  C.  there  the  following 
definitions  (for  indefinite  integrals)  are  Introduced 


iw)  -  J  w®(w^  ♦ 


--3/2, 


w)  -  j  w^'Cw'^  ♦ 


2  .  ,^-3/2 


(B.15) 


They  are  related  by 


l“^^^(w)  -  -w^'^Cw^  +  1)“’’''^  +  (m-1  )I~''^^(w) 
m  m-2 


First  we  evaluate  ^21’  Setting 


(B.16) 


m  -  2n  +  1 


one  obtains 


r-3/2,  ,  2n,  2  ,^-1/2  „  T-l/S,  , 

.  (w)  -  -w  (w  1)  +  2n  I-  .(w) 


(B.17) 


Furthermore,  from  the  equation  top  of  page  149 


1  /2 


n-1 

I  ” 
*.-0 


28, 


For  the  821'®’  recurrence  relation  Eq.  (C.8). 

obtain  U«  one  must  substitute  the  limits  w  -  V  and  w  - 
2n+l 

Eqs.  (B.17)  and  (B.I8),  and  multiply  by  ((2n+1)!)~  .  One 


-(V^  +  1)  -  (-1) 


U 


2n+1 


■[(2n  +  1)l]  ^ 


+  c 


2n+1 


(V^  * 


'2n  +  1  ;  n  > 


with 


2n 


'2n  +  1  "  (2n'H  )  l(2n-1 


2n-2 


One  has  the  following  recursion  relation 


°2n+1  2n  2n-3  (2n-1 ) !  ^2n-4 


C2^_^  2n^  STTT  (2n  +  1  )  !  020-2 


One  has  from  Eq.  (C.8)  (setting  k 


‘2n-4 


n-2) 

2n-2 


0 


2n-2 


2n-3 


Therefore 


1  /(4n  -1  ) 


(B.18) 

To 

0  into 
obtains 

(B.19) 

(B.20) 

1 

(B.21  ) 


(B.22) 


V.  WL^  L’^  JWVTW_”WV7W^  V-W  IT-fc  _ih  ki~m  *.-»  , 


Because  of  the  factor  2n  in  the  numerator  of  Eq.  (B.21)  the 
recursion  starts  with  n  -  1 ;  then 

^  2  1 
°3  “  31  “3 

This  suffices  to  evaluate  Eq.  (B.20).  But  on  the  basis  of  Eq. 

(B.22),  one  can  also  start  with  n  -  0  and  c^  -  -1 .  Then  Eq.  (B.20) 

encompasses  also  Eq.  (B.19),  if  the  sum 
k-0-1 

I  which  then  arises  is  understood  to  be  zero. 
k-0 

The  expression  consists  of  three  terms:  the  sum  in 

Eq.  (B.9)  formed  with  the  first  of  these  terms  gives 


V  /  xn  1  r2n+1  ^2n,r,2  .  . 


-1/2 


-  V"’ (V^  +  1 )“’^^sin(kV) 


One  thus  obtains 


§2^  -  -  1)^^^  sin  (kV) 


(V^  .  ,)’'2  (-)■> 


(B.22a) 


,-2n+1 


-  y  C-)'^  c  k 
n-0  2n>1 


The  C2n  +  1 ®2J,*®  evaluated,  respectively,  from  Eq.  (B.22) 
with  c^  -  -1  and  from  Eq.  (C.8)  with  Bq  -  1 . 

To  evaluate  B^j^  one  needs  expressions  Ugj.  They,  too,  are 
baaed  on  the  formulae  of  Appendix  C.  Eq.  (C.13)  for  m  -  2m^ 
assumes  the  form 


•"1  ^  ^ 
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m 


m 


.V, 


#■  J 


Examining  this  equation  and  Eqs.  (C.9),  (C.10),  and  (C.11)  one 
finds  that  there  is  no  contribution  of  the  lower  limit  in  Eq. 
(B.10).  One  then  finds 

Uq(V)  -  -  V“^ (V^  +  (V^  +  1)^^^  (B.23) 


U2(V) 


^2n(V) 


v"^  (V^  +  1)“"'^^  ^  ^  log  [(V^  +  1)^'^  +  v]  (B.24) 


2! 


'2n' 


J.-1 


-  C-  log[(V^  ♦  1 ♦  V]  n  >  2 


1/2 


'2n 


(B.25) 


with 


2n-1  1 

°2n  “  (2n) I  (2n“2)B 


(B.26) 


2n-3 


The  Cgj's  are  expressed  explicitly  in  Eq.  (B.27).  This  formula  is 
obtained  as  follows: 


cjj  -  1/16 


(B.26a) 


From  Eq.  (C.5)  which,  by  definition,  is  also  satisfied  by  the  B's 
one  finds 


B 


2n-3 


B 


2n-1 


2n-1 

2n-2 


Then,  one  obtains  from  Eq.  (B.26)  the  recurrence  relation 


°2n*2/°2n  ‘  '  t2n(2n*2)] 


-1 


with  the  initial  condition  (B.26a)  this  is  solved  by 


°2n  ■  2(-1/'l)"  [nKn-D!]"'' 


(B.27) 


•i--' 


'> , 


'•j.v 


Pi 


■■  •-  .j: 
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For  n-1 ,  this  gives 


°2  "  "  ? 

The  last  results  show  that  Eq.  (B.25)  covers  also  Eq.  (B.2^1)  if  the 
i-0 

sum  I  is  understood  to  be  zero. 

*,-  +  1 

Thus  one  obtains 

§2^  -  -  cos(i<V)  +  +  1)^^^ 


+  2(V^  + 


^2  nTTn^TTr  ®2t-i  ^ 

n-2  l-i 


2  log((v"  .  1)^^"  .  V)  XHTtW 


(B.28) 


with 


®2t>1^®2l-1  "  “  *  (1/2))  ,  B,  -  1 

The  independent  variable  in  B^  is  k.  To  suggest  this  we  write 
temporarily 


k  -  y 


B^  -  Q 


(B.29) 


Then 


Q(y)  .  r  -.exp(.izu) 

(u^  ♦  yy'^ 


(B.30) 


Q(y)  satisfies  a  differential  equation  closely  related  to  the 
Bessel  equation.  One  has 


u  exp(lyu) 


aw  ,  f  u  expv 


(B.31 ) 


g  w  _  _  f  u  exp 
dy^  '  (1  ♦  u 


9  exp(tyu) 


2j3/2 


(B.32) 


2  2 

d  Q/dy  is  not  absolutely  convergent;  but  the  integral  is  well 
defined  if  one  reformulates  Eq.  (B.30)  as 


Q(y) 


exp( iyu) 


a-»«  i  (1  +  u 


2^3/2 


i:- 


One  obtains  from  Eqs.  (B.30)  and  (B.32) 


exD( iyu) 


f  expi ly 
Cu^  *  1 


and  from  Eq.  (B.31)  by  an  integration  by  parts 


(B.33) 


I  -  y  1  1“  (B-S*) 

♦  1)^'^  (u^  ♦  1)’'2 


Combining  Eqs.  (B.33)  and  (B.3^)  one  obtains  the  following 
differential  equation  for  Q 


d"Q  1  dQ  p  i 
(ly2  y  dy  “  y 

We  separate  Q  into  Ints  real  and  imaginary  parts  Q 
Then 


(B.35) 


Qr  ^  i  Qj. 


1 

-  -  Qr.  -  0 

y  dy 


1  1 
J.  _ i.  -  0  -  _ 

y  dy  y 


(B.36) 


(B.37) 


Obviously  Q(*)  -  0,  furthermore 


Q(0)  -  f  — g --— 


(u*"  1) 


(B.38) 


These  are  the  boundary  conditions  for  Eq.  (B.35).  The  solution 
will  arise  in  the  form  of  power  series  in  y,  which  makes  it 
difficult  to  recognize  whether  the  boundary  condition  QC-)  -  0  is 
satisfied.  It  will  be  shown  below,  that  the  definition 


Qi(y)  -  I 


(B.39) 


Implies 


(B.40) 


The  derivation  is  shown  below.  This  condition  will  be  used  instead 
of  Qj(«)  -  0.  This  procedure  cannot  be  applied  to  Q^;  (there  one 
would  find  d^Qj^|dy^  |y_Q  -  From  here  on,  the  real  and  imaginary 
parts  are  treated  separately.  The  following  derivation  of  Eq. 
(B.40)  ends  with  Eq.  (B.47).  One  has 


-  / 


^  u^  sin(yu) 


du  -  -Q,.  (y)  -  Qi2(y) 

(y^  >  1) 


(B.41 ) 


with 


,  7  sin(yu)  j.. 

*Ti  ■  — 5 *'  (  ?5  au 


r  sin(yu) 


(B.42) 


/  D  II  O  ^ 


It  Is  not  permissible  to  form  the  limit  y  -  0  under  the  integral 
sign,  for  within  the  region  of  integration  there  are  always  values 
of  u  for  which  sin(yu)  is  not  small.  We  write 


-a/y 


‘12 


/  •  ,  a  >  0 

i  A,.,  •  1)3'^ 


-a/y 


Now 


I T  ^  - 1  <  T  ^ 


-a/y 


(u"  1) 


372 


du 


u 


(u2  .1)1/2 


-a 


T~T7T7? 


+  1 


(a‘  +  y") 

This  expression  vanishes  for  y-*0.  Furthermore, 


sin(yu) 


-a/y 


(u^  1 


du  < 


'  {  (u2  ^,13/2 


-a/y 


(u^  ♦  1) 


du  -  y(u^  +  1 )“^ 


u*- 


-  yCl  -  ((a/y)^  +  1)“''^^] 


Therefore 


lim  Q  (y)  -  0  (B.44) 

y-»0  ’ 


is  transformed  as  follows 

Qjl  -  J  ^  2^ ^^^^1/2  “  ^)3in(yu)du 


-  1 

I  u  sin(yu)du  (B.45) 
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•  *  V  •_*  •-*  -u"  -  '  •  ■  •  •  ^ 


The  second  integral  can  be  evaluated  by  the  calculus  of  residue, 

0  +« 

I  u~'  sin(yu)du  -  (1/2)  |  v” '  (sinv)dv  -  ir/2  (B.J16) 


The  integrand  in  the  first  integral  on  the  right  of  Eq.  (B.45)  is 
rewritten 


(u2  .  l)-1/2  ^  -1  .  (u^  >  u _ 1 _ 

^  ^  (u2.  (u^  .  1)’/V(u2  .  1)1/2 


In  the  integrand  of  the  following  expression,  one  has 


u  <  -  a/y  <  0 


Therefore 


I 


ain(yu)du 
1)1/2uC(u^  ♦  1)' 


(y/a)  I 


.  1 


The  limit  y  -  0  of  this  expression  is  zero.  Furthermore, 


I  77T1 

•a/y 


sin(yu)du 

,1/2  r  /  ^  ,  ,1  /2 

)  u[ (u  ♦  1  ) 


—  h »  I 


(B.1I7) 


Also  this  integral  vanishes  in  the  limit  y+0.  Substituting 
Eqs.  (B.44),  (B.45),  and  (B.46)  into  Eq.  (B.41)  one  obtains, 
indeed,  Eq.  (B.40). 

One  notices  that  the  homogeneous  part  of  Eq.  (B.37)  is  the 

same  as  in  Eq.  (B.36).  The  indicial  equation  for  the  singular 

point  y  -  0  in  Eq.  (B.37)  gives  the  exponents  0  and  2.  Therefore, 

there  exists  a  homogeneous  solution  in  the  form  of  a  power  series 

2 

which  starts  with  y  .  We  set  for  this  particular  solution 


•'>y  y-  • 
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"  I  ^2j.2  y 


(B.48) 


This  gives  the  recurrence  relation 

(2>2)2J  -  ajj  -  0 

which  is  satisfied  by 

®‘2j+2  “  (j  +  1  )!]“^  (B.l»9) 

if  one  chooses  a2  -  1 .  One  solution  of  the  inhomogeneous  equation 
(B.37)  can  be  obtained  as  a  power  series  in  y. 

Q4(y)  -  I  (B.50) 

•J  ^ 

Hence,  by  substitution  into  Eq.  (B*37) 

l  (-)'^C(2jM)(2j-1)b2j^1  ♦  ^2j-1^y^'^~^  “  y"^ 

j  *  0 

The  equation  gives  for  j  -  1 

-b^y-'  -  y"’ 

Hence 

b^  -  -  1  (B.51 ) 

In  addition  one  has  the  recurrence  relation 

^2j  +  /^2j-1  "  "  J  +  W2)(  j-1/2)]"''  -  -[(2j)^  -  1]“''  (B.52) 

Then  because  of  Eq.  (B.i<0) 


QjCy)  -  QhCy)  *  If/**  Qo(y) 


(B.53) 


Returning  to  the  original  notation  one  has 


B 


1 1 


(k) 


I  ^  ♦  If  I  (k/2)^J  +  ^ 

j-0  '  J-0 

(B.54) 


Now 

Bj(V,k)  -  B^j(k)  +  B2j(V,k) 

Bgj  is  found  in  Eq.  (B.22a).  The  last  term  in  Eq.  (B.22a)  cancels 
the  first  term  in  Eq.  (B.54),  for  the  and  the  C2j^^  satisfy 

the  same  recurrence  relations  (Eqs.  (B.52)  and  (B.22)),  the  same 
starting  values  -  c^  *  “1 .  and  the  sums  have  opposite  signs. 

One  thus  obtains 


Bj  -  -  V~\v^  ♦  sin(kV) 


(B.55) 


+  IT  I  (k/2)^'^*^[n!  (n+1  )!] 
J-0 


-1 


where  the  coefficients  C2^^^  and  the  are  obtained  respectively 
from  Eq.  (B.22),  with  c^  -  -1  (or  c^  -  1/3).  and  Eq.  (C.7)  with 
Bq  -  1. 

To  evaluate  Qp(y)  we  use  a  result  about  Bessel  functions. 
Setting 


QR(y)  -  y  z(y) 


one  obtains  from  Eq.  (B.36) 


d^Z  1_  ^ 

77  y  ay 


(1  *  ■  0 

y 


(B.56) 
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which  is  solved  by 


Z  -  Z^(iy) 

where  is  some  linear  combination  of  Bessel  functions  of  order  1 
the  conditino  Z+0,  as  y+0  is  satisfied  by  the  Hankel  function  of 
order  1,  (iz).  One  has  (see  for  instance  Ref.  3) 

H^(2)  -  J^(z)  ♦  i  Y^(z) 


where 


z  -  X  ♦  iy 


,  /_^  "  (-(z/2)^)" 

J^(z)  (z/2)  -n!(n+iyr' 


(z) 


Here 


Then 


-  I  -(z/2)“^  ♦  2(log(z/2))J^ (z) 

-  (z/2)  (*(„M)  .  ♦(„.2))i-m^} 


,2sn 


~  Y 


n-1 


K  •  -  y 


- 1 

♦  I  m  ,  n  >  2 


m-l 


Y  -  .5772156649.. 


n  ,  - 

-  -  2Y  ♦  2  I  m" '  +  (n+1 ' 


(B.57) 


Observing  that 


log(iy/2)  -  log(y/2)  ♦  lir/2 


one  obtains 


H^(ly)  -  J^(iy) 
-  {-(y/2)‘‘' 
According  to  Eq 


+  i  y^(iy) 

+  I  (♦(n+l) 
n-0 

(B.56) 


+  ij)(n+2) 


21og(y/2) ) 


n! (n^l ) ! 


Qp  -  const  y  ( iy) 


and  this  expression  satisfies  the  condition  for  Satisfying 


the  condition  Qply^g  “  ^  returning  to  the  original  notation  one 


then  obtains 


(k/2) 


2n+2 


B.j,(k)  -  1  -  I  (iKn+l)  +  i|;(n+2)  -  21og(k/2)) 
n-0 


With  Eq.  (B.28)  one  then  obtains 


Sr  -  B,r  (S)  '  Brr 


B.  -  -  5"'(V^  ♦  003(kV)  ♦  v"’(v^  ♦  1)'^^  •  1 

R 


(B.58) 


*  *  0^/2  y  (k/2)  y  g 

*  2(V  ♦  1)  I  n!  (n>1  )!  ®2£-1 

n-i  *"-• 


2n+2  n 


^21-1 


n-0 


(k/2) 


2n+2 


(k/2) 


2n+2 


-  I  *  K.(n+2)  -  21og(k/2) 


(Here  the  summation  subscript  In  the  second  row  of  Bq.  (B.28)  has 
been  changed  by  setting  j  -  (n+1).  The  next  step  would  be  to 
return  to  the  original  variables  k  and  r,  and  to  form  B  -  Bp  i  B 

The  first  terms  in  these  expressions  then  combine  to  form  as 
contributions  to  Bp  ♦  i  Bj. 

-V"l  (V^  +  1  expdicV) 


or  because  of  Eq.  (B.2) 


(B.59) 


-  V~’(V^  +  exp(ikV) 


Within  the  kernel  K,  B  appears  in  combination  with 


M  R*^(V^  +  exp(lkV) 


This  and  the  expression  Eq.  (B.59)  can  be  combined.  We  recall  the 
following  relations  (Eq.  (2)) 


r  -  [(Ji-y)^  ♦ 

R  -  [(^-x)^  ♦ 

V  -  -b‘^((^-x)  ♦  MR) 


(B.60) 


Then 


(V^  +  r^)^''^  -  b"^(R  ♦  M(t“x) 


log(V^  -  V)  -  log[B'^(R  +  M(C-x)  +  (^-x)  ♦  MR)] 

-  log[(1-M)“^  (R  5-x)]  -  log[(1+M)r^(R  -  (^-x))"''] 


(B.eOa) 


Using  these  relations  one  obtains 


exp(ikV)(V^  +  r^)"^^^  f-  r  ♦  rl-  “  rr  exp(ikV) 


1 


•  • 


-2  -2  -1  /2 

The  vanishing  of  the  factor  (V  +  r  )  is  a  significant 
simplification  for  subsequent  integrations  with  respect  to  K  and  n. 
The  factor  V  ^  occurs  only  in  the  first  term  of  a  development  with 
respect  to  k.  To  display  this  we  write 


exp(ikV)(V^  +  r^) 


1/2  fi  .  Ml.  .  1_  c^,p(i,v)  -  1]  - 

Lv  rJ  rv  rv 


The  term  in  the  bracket  la,  of  course,  of  0(k).  Combining  the  term 

-  (RV)~^  with  the  other  terms  of  order  k®  (in  Eq.  (B.58)),  one 
obtains  after  introduction  of  the  variables  k  and  r,  using 
Eqs.  (B.60)  and  (B.6l) 

-  1_  .  i-  .  1  (V^  ♦ 

RV  r^  ?  V 

1  1  r  -2  5,r,2  ~2.^/2^ 

-  -=■  +  “5 —  I-  r  +  R(V  +  r  )  J 

r^RV 

-  ^  ♦  -^-4 —  [-  ♦  R^  +  MR(E-x)] 

r  6  r  RV 

■  ^  *  -s-?",-  [-  +  B^r^  +  MR(E-x)] 

e'^r'^RV 

Hence 


RV  r 


1  (V2  . 


(l-x) 


R(R  +  ^-x) 


The  last  formulation  shows  that  for  ^  -  x  >  0,  and  r  -  0  (i.e.,  if 

the  point  (x,y)  lies  upstream  of  the  point  (C»h)»  there  is  no 

singularity  in  the  flow  field,  and  for  (^-x)  <  0  (i.e.,  if  the 

point  (x,y)  lies  in  the  wake  of  the  point  (Cfh)).  one  has  a 
~-2 

singularity  -  r  .  Assembling  all  this  information  one  obtains. 


-•W..VJV^V 


exp(ikV)  +  Bj^(k,V)  +  i  Bj{k,V) 


-  -  [RV]""' (exp(ikV)  -  1  )  +  R  ~— 


^  2(v2  ^  y  (l</2)^^^^  y  y2a-1  -2(n--0 

^  n!(n>1)!  ®2£-1  '' 

n-i  £*i 


+  2(k/2)2[l0g((k/2)((V^  + 


(k/2)2  [i|>(n  +  1)  +  4-(n^2)  -  iir] 


.  i(v2  .  ?2j1/2 


°2n.1 


£-0 


Here  the  expressions  containing  V  can  be  expressed  in  terms  of  ^-x 
and  R,  by  Eqs.  (B.60)  R  in  turn  is  expressed  by  ^-x  and  r.  The 
expression  therefore  appears  in  terms  of  the  basic  independent 
variables. 

In  the  last  minor  step,  we  introduce  the  Prandtl-Glauert 
coordinate  distortion.  The  purpose  is  to  simplify  the  expression 
R,  which  introduces  in  certain  integrations  branch  points  of  the 
integrand.  We  set  accordingly: 


X  -  X ,  5  ■  C 

y  -  B  V,  n  -  B  z  -  3  (B.6l) 

r  -  3~^r,  r  -  [(n-y)^  ♦ 

R(e-x,r)  -  R(C-x,r)  -  ((C-x)^  + 

V(5-x,r)  -  V(5-x,r)  -  -8'^(MR  +  ^-x) 

(V^  >  -  3"^CR+M(5-x)]  -  b“^[R>M(C-x)] 


mk 


log(V^  ♦  -  V)  -  log[(1-M)~^  (R-^e-x))  -  log[(1-M)"^  (R+^-x)] 

-  log[(1-M)‘''  r^CR  -  (5-x))“'']  (B.61) 

(cont'd) 


The  basic  expression  for  the  upwash  now  appears  in  the  form 


w(x,y,z)  -  (8x)  ^  jf  Ap(5,n)K(k,C“X,r)d5dn) 


(B.62) 


K  -  8  exp(ik(C-x)(K^  +  K^) 


(B.63) 


K,  ■  ^  .  (k=/2,[log(k/2)r=(R  -  (e-x))-’]  I 


(B.6M) 


-  28^r“^  -  8^CR(R  -  (C-x))]“^  +  (k^/2)[log( (kr)^/2) 


-  log(k(R  -  (f-x))]  I 


(B.6Ha) 


-  8‘^R‘^(R  ♦  (t-x))"’  ♦  (k^/2)log[(k/2)(R  +  (^-x))]I 

(B.6Mb) 


-  (RV)"'  exp(ikV)  -1) 


(B.65) 


♦  8"^(R+M(  -x))  2  X 

n-1 


(k/2) 


2n+2  n 


,2t-1  x2(n-l) 


im-  J,  »21-1 


i  J  (-)"  0,„.,  k‘ 
n-1 


i  6,. 

l-O 


-  (k/2)^  [ip(n+1)  +  4»(n+2)  ♦  log(l-M)  - 


where 


c _ ,  -  -1/(4n  ♦  1),  c,  -  1 


During  the  derivations  also  constants  C2^  have  been  encountered, 
but  in  this  formula  they  are  directly  expressed  by  factorials  (see 
Eq.  (B.27)).  Moreover 


®2n*2'^®2n 
®2n+1 ^®2n-1 
4i(n*1 )  ♦  i(»(n+2) 


(2n*1 )/(2n+2) ,  Bq  -  1 
2n/(2n+1 ) ,  B,  -  1 

2Y  ♦  I  +  (n+1 )“^ 

m-1 


Y  -  .5772156649 


r,  R,  and  V  are  defined  in  Eq.  (B.61).  Eq.  (B.64a)  displays  the 
singularities  which  occur  for  (5“x)  <  0  and  r+0.  Eq.  (64b)  shows 
that  there  are  no  singularities  for  (C“x)  >  0.  For  points  (x,y) 
close  to  points  (C,n),  for  instance  for  points  (x,y)  within  or 

close  to  a  (5,n)  element  only  very  few  terms  of  the  infinite  series 

2 

are  needed.  Including  powers  up  to  k  one  has 

-  B^r“^R‘\R  -  (C-x))  ♦  (k^/2)log  kCr^(R  -  (C-x))“^]  (B.66) 

-  B^R‘^(R  (C-x))“^  +  (k^/2)  log  k(R  (C-x)) 

Kg  -  -  IkR*’  ♦  (k^/2)R~W 

-  (k/2)^[«(1)  +  i|»(2)  +  2  log  2  ♦  log(l-M)  -  in] 


(B.67) 


APPENDIX  C 

SOME  INDEFINITE  INTEGRALS 


Let 


-1/2 


m 


=  dw 


(C.l) 


Setting 

I  =  (w^  +  I  a™  +  e®  f(w^  +  dw  (C.2) 

m  g  K  j 

one  has 

-1/2/.  /  2  ^  -v-l/2r'"v^  k+1  ^  ,k-l  ^  m, 

dl  /dw  *  (w  +1)  [  2,  a.  (k+l)w  +  I  a.  k  w  +  c  ] 

™  k=0  ^  1  ^ 

(C.3) 

We  replace  k  in  the  first  sum  by  j-1  and  in  the  second  sum  by 
J+1.  Then 

dljjj“^^^ldw  =  (w^  +  aj^.^(j+l)w^  +  o'"] 

(C.4) 

For  1  <  J  <  m  one  obtains  by  comparison  of  the  coefficients  of 
the  power  w*^  in  Eq.  (C.3)  and  in  the  derivative  of  Eq.  (C.l) 


J  -  0 


(C.5) 


Moreover,  from  the  power  w 


0 


„m  .  m  - 
+  c  *  0 


(C.6) 


and  from  the  power  w 


m 


m  aj"  =  1 
m— i 


(C.7) 


1-47 


Eq.  (C.4)  amounts  to  two  recurrence  relations,  one  for  even  and 
one  for  odd  values  of  m.  Eq.  (C.7)  shows,  that  the  subscripts 
of  the  will  be  odd  for  even  m,  and  even  for  odd  m.  Consider 

odd  values  of  m  first  and  set 

m  =  2m^  +  1 


then  from  Eq.  (C.7) 

2m^ +1  ^ 

Let  a  sequence  6^^  with  6^  =  1  satisfy  the  recurrence  relation 
Eq.  (C.5)  namely 

®2k+2/®2k  “  “  (2k+l)/(2k+2)  (C.8) 

The  B2k*®  expressed  in  terms  of  factorial  for  half-integral 

argument,  but  this  has  no  practical  significance,  as  one  will  even 
then  evaluate  the  82j^*s  recursively.  The  first  few  of  the  32jj’s 
are  found  in  the  following  table 

0  1 

2  -1/2 

'I  3/8 

6  -5/16 

8  35/128 

Because  of  Eq.  (C.7a)  one  has 

2m,  +1 


=  [(2m,+l)R  „  ]-1r  „ 


Therefore,  from  Eq.  (C.2) 


m-i 


-1/2  _ 

2m^+l 

We  write  down  some  of  these  expressions 


-1/2  _ 


(w2  +  1)1/2 


-1/2 


(w^  +  l)^/^[(l/3)w^  -  (2/3)] 


-1/2 


(w^  +  l)l/^[(l/5)w^  -  (i4/15)w^  +  (8/15)] 


-1/2 


(w^  +  l)l/^[(l/7)w^  -  (6/35)w^  +  (8/35)w^  -  (16/35)] 


For  m  even,  we  set  m  =  2m^.  Then 

2m 


2m 


1.1  =  (2mi) 


-1 


Setting  J  =  2k  in  Eq.  (C,5)  one  obtains  the  recurrence  relation 


^2k+l/®2k-l  "  "  (2k)/(2k+l)  ,  k  >  1 


we  choose  6^^  =  1.  Then  one  obtains  the  following  values 


i 

1 

3 

5 

7 


e 


jt 


1 

-2/3 

8/15 

-^<8/105 


=  -16/35 


2"^1  -1 
^2k-l  "  '-^"’l^2mi-l^  ^2k-l 


Then 


Finally,  from  Eq.  (C.6) 


2m,  2m, 


-1 


Thus , 

(C.9) 

=  |(w^  +  l)~^^^dv  =  log[(w^  +  +  w]  (C.IO) 

Some  such  expressions  are 

^  („2  ^  1)  1/2 3^/2  )w] 

=  (w^  +  l)^^^[(l/4)w^ 

=  (w^  +  l)^''^[(l/6)w^ 

Expressions 

1-3/2  ^  f  „in(^2  ^  i)-3/2  (C.12) 

m  J 


-  (1/2) 

-  (3/8)w]  +  (3/8)  (C.ll) 

-  (5/2H)w^  +  (5/16)w]  -  (5/16) 


could  be  treated  independently  by  a  similar  procedure.  In  the 

present  context,  one  is  led  to  more  useful  formulae  for  the  upwash 

-3/2 

if  one  carries  out  an  integration  by  parts  to  express  I  in 
-1/2 

terms  of  I  .  One  has 
m 


r-3/2  _ 


m 


m-1,  2 
-  w  (w 


+  1) 


-1/2 


+  (m-1)  I 


,-1/2 


m-2 


m  >  1 


(C.13) 


For  m  =  0,  one  obtains  directly 


I 


-3/2 


w(w  +  1) 


-1/2 


(C.14) 


Also  needed  Is  the  following  Integral 


I  =  |(wcosa  -  slnot)“^(w^  +  l)~^^^dw 


Setting 


w  =  tg0 


(C.15) 


one  obtains 


I  =  slgn(cos0)|[(tgecosa  -  sina)cos0]~^d0  =  sign(cos0 )|sin(0-a)“^d0 


Then 


|(wcosa  -  sina)“^(w^  +  l)”^'^^dw 


(C.16) 


-1- 


sign(cos0)log[sln(0-o) (1  +  cos(0-a))”  ] 


Furthermore 


|(wcosa  -  sina)”^(w^  +  l)”^'^^dw 


=  sign(cos6) jsin(0-a)”^[cos(0-a)cosa  -  sin(0-ot)sina]d0 


Hence 


|(wcosa  -  sina)~^(w^  +  l)~^^^dw 


(C.17) 


*  -  slgn(cos0 ) {[cosa  sin(0-a)I^  +  sina  log[sin(0-a) (1  +  cos(0-a) )~^]} 


r 


f^(T))dT) 


In  both  integrals  we  replace  f^Cn)  by  the  maximum  of  If^Cn)!  - 
|f^|max.  One  has 


dn  , 

I  (n-y)^  ♦  z^ 

^1 


-1 

z 


arctg 


As  the  arctg  remains  bounded  even  when  its  argument  tends  to 

infinity,  (as  it  does  for  z+0)  the  first  term  is  0(2)  because  of 
2 

the  factor  3z  in  front.  Moreover 


I 


dn 

(n-y)^  ♦  z^ 


z(n-y) 
(n-y)  ♦  z^ 


♦  arctg 


Also  the  second  term  is  0(z). 


APPENDIX  E 


THE  EVALUATION  OF  INTEGRALS  Jq"4»^dq,  |q'’()>2dq,  AND  Jq"(|»j^dq 

First  we  derive  some  auxiliary  relations.  Considering  U 
V  as  functions  of  p  and  q  (Eqs.  (I6i»))  one  has 

3U/3q  -  sin  6 

3V/3q  -  cos  6 

Furthermore,  since  R  -  (U^  + 


-R 


-1 


U  -  const 


Itt  10g(R-V)|  -  (R-V)  ’uR  ^  -  (R+V)U  ^R  ^ 

V  -  const  .  .  , 

-  U"  +  VU  R 


Therefore, 


|-  log(R-V)  -  -r‘^cos  6  *  R“’vu‘^sin  6  +  if^sin  6 
3q 

-  r“^u'^(-Uco3  6  +  V  sin  6)  ♦  u“^sin  6 


Then  with  Eqs.  (16U) 


1^  log(R-V)  -  -pR"’u“’4.  u'^sin  6 
1^  log((R-V)  jU)"^  )  -  -pR'''u“‘' 


(E. 


Similarly 


1^  log((R-X)|Y|"’  )  -  pR‘"'y*'' 


(E. 


and 


1) 


2) 


(Notice  the  difference  in  sign.) 


Consider  now 


|q"4»idq  -  jq'^log(k(R-X))dq  -  |q"[log(  (R-X)/Y)  +  log(kY)]dq 
With  Eqs.  (E.2)  and  Eqs.  (162),  viz. 

Y  -  q  cos  Y  -  p  sin  Y 


one  obtains 

Jq”^^dq  -  (n+1)  ^  log(k(R-X) )  -  pjq'^^^R  ^dq  -  cosYjq’^^^Y  ^  dq 
Carrying  out  the  division  one  finds 

-  cos“S[q"  +  (ptgY)q"“^+{ptgY)^q''‘^  . . .  (ptgY)”]  +  (ptgY)"*^  Y“^ 
There  are  two  integrals  with  the  factor  (ptgY)’^*\  namely 
-  (ptgY)”*Hp  jR”^  Y“^ dq  +  cos  Y  jY^^dq] 

Hence  with  Eq.  (E.2) 

-  -(ptgY)"'"klog((R-X)|Y|"b  +  log(kY) 

-  -(ptgY)"^^  log(k(R-X)) 

Therefore , 

|q%^dq  -  (n+1  -  (ptgY)"*’ )log(k(R-X) )  (E.3) 

-  p  cos'Nc  Jr“^  (q"  +  (ptgY)q"'‘^  . . .  (ptgY)'^]dq 

-  [(n  +  1)“’q"^^  +  (ptgY)n"\"  +  (ptgY)^(n-l  )“''q""\ . .  (ptgY)"q] 

The  treatment  of  jq^(^2dq  ■  log(R-V)|U|  is  nearly  the  same. 
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Beside  Eqs.  (E.1)  one  has  Eq.  (164),  viz. 


U  -  q  sin  6  +  p  cos  6 

jq^^gdq  -  (n+1  )"Mq'^*’log(R-V)  |U|“M  ♦  pjR'"' (qsin6  +  pcos6)‘''dq} 
q”*^  (qsin6  +  pcos6)~^  -  sin~^ d[q”-*-(-pcot6 )q*'  ^  +  (-pcot6)^q’™  ^ 

+  (-cot6p)*^  +  (-pcot6)"*^U~^ 

Hence , 

|q”*2dq  -  (n+1  1  (q"*^  -  (-pcotY)"''hog(  (R-V)  |U  |“’ ) 

(E.4) 

+  sin'^6jR“kq'^+(-pcot6)q"‘’+(-pcot6)^q"‘^  +  (-pcot6)"]dq| 
Finally  we  consider 

Jq”*4dq  -  |q"Ry“"'dq  -  ICp^q"^  ♦  q""^  )Y"'' R""*  dq 

Now 

p^q"Y“^  -  Pacos’* WCq"“’  +  (ptgY)q""^  ♦  (ptgY)"]"’  +  p^ (ptgY) "y“^ 

_  cos’SCq"*^  ♦  (ptgY)q"  ♦  .  .  .  (ptgY) ]  +  (ptgY)’^'"^Y“’ 

One  notices  that  terms  with  the  same  power  of  p  can  be  combined. 

For  instance 

n+2,.  ^vn  -2^ 

-  P  (tgY)  cos  Y 


2/  ..  ,  .  ^^n+2 

p  (ptgY)  ♦  (ptgY) 


(E.5) 


but  in  general  the  result  appears  rather  more  complicated. 


Jq"<|)ndq  -  cos'NjR'kq”'^’  +  (ptgYq”  +  (ptgY) ^q""\ .  .  (ptgY)""’’ 

+  P^  q""^ - p^(ptgY)'’~^  ]dq  ♦  p"""^  (tgY)"cos~^Ylog((R-X)/|Y|“b 


In  the  last  term  Eq.  (E.2)  has  been  applied. 


APPENDIX  F 
LIMIT  o-»Tr/2 


Only  some  of  the  expressions  for  the  upwash  written  down  for 
a  f  ir/2  allow  one  to  form  the  limit  a+ii/2  directly.  The  others 
give  either  infinity  or  the  difference  between  two  quantities  that 
tend  to  infinity.  The  results  for  o  »  it/2  have  been  obtained  by 
direct  computation.  This  suffices  for  practical  work.  The 
limiting  process  a+ir/2  has  some  value  as  an  exercise,  besides  it 
provides  a  cross  check  of  the  formulae. 


The  type  of  difficulty  can  be  seen  in  the  following  examples 
Consider 


dx 

2x  +  e 


in  the 


limit  e+0. 


/ 


dx 

ex+c 


One  has 

£  ^  log(ex+c) 
x/c 


for  e  /  0 
for  e  -  0 


Now 


e  ^  log(ex+c)  •  e  ^  logc+e  ^  log(1  +  (ex/c)) 

The  first  terra  can  be  regarded  as  a  constant  of  integration, 
tends  infinity  as  e  tends  to  zero.  Developing  the  logarithm 
respect  to  e  one  obtains  indeed  x/c. 


which 

with 


The  example  shows  that  the  idea  of  a  development  cannot 
always  be  avoided,  even  if  one  would  include  the  specific  constant 
of  integration  in  the  formula  for  the  indefinite  integral. 

We  derive  a  number  of  recurring  limiting  expressions.  The 
order  of  the  error  in  terms  of  Aa  -  a  -  ir/2  is  shown 


(1 


-2  2 
sina)cos  a  -  (1/2)  +  0(Aa) 

cos  ^a-tga  »  0(Aa) 


(F.  1  ) 
(F.2) 
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Using  Eqs.  (133)  one  obtains 


U  -  -Y+0(Aa) 
V  -  X+0(Ao) 


(F.3) 


More  specifically, 

V-X  -  Ycosa-X( 1 -sina)  -  Ycosa- (1/2 )Xcos^a+0 ( Aa) ^  (F.4) 

U+Y  -  Xcosa+Y(l-sino)  -  Xcosa+(l/2) Ycos^a+0 ( Aa) ^  (F.5) 

U"^  -  -Y"^+U"^+Y“^  (F.6) 

U"^  -  -Y“^+(U+V)U"^Y"^  -  -Y"^-XcosaY“^+0( Aa)^  (F.7) 

VU~^  -  (Xsina+Ycosa) (Xcosa-Ysina)  ^ 

-  -XY"^  +(XY"’  +  (Xsina+Ycosa) (Xcosa-Ysino)"^ ) 

=  -XY"^ +(X^cosa-XYsina+XYsino+Y^cosa)  (P.8) 

VU"^  -  -XY"’-cosaR^Y“^+0(Aa)^  fv  o'. 


UY"^  -  -1+(U+Y)y"^  -  -1+cosaXY~’+(1/2)cos^a+0(Aa^) 
log((R-V)  (R-X)"b  -  log(1-(V-X)(R-X)"^  ) 

-  (-Ycosa+(l/2)Xcos^a) (R-X)"’  -  ( 1 /2 ) Y^cos^a( R-X)~^ +0 ( Aa^) 

log((R-V)  (R-X)“^  )  -  -C0SQt(RY“^+XY"S  -  ( 1 /2  )cos^a[  R^  y"^  +  RXY"^  ] +0  (  Aa) 


We  revert  to  the  original  notation  and  write  down  those 
expressions  which  contain  factors  tga  or  cosa 

IqqCo)  -  B^[arctg(v/u)  +  tgo  log(l<(R-X))  -  cos  \log(  (R-V) /U)  ] 

01  2  '■‘2  — P  2 

I^Q  -  B  [-sina  cos  a  U  log(kY)  +  cos  a  U  /y] 

IqI  -  -B^cos"’a  U  log  (KY) 

02  2  2  —  P  1 

-  B  {u  tg  a  log{(R-V)/U)  -  sina  cosa  log((R-X)/Y)  +  cos'  a(R/Y)} 

02  2  —  1 

-  B  lu  tga  log((R-V)/U)  -  cos  a  log(R-X)/Y} 

Because  of  the  definitions  of  X,  Y,  u,  and  v  one  has 
arctg(v/u)  »  -arctg(X/Y)  +  const.  This  term,  expressed  in  U  and  V, 
differs  from  the  corresponding  expression  in  terms  of  X  and  Y  by  a 
constant  only.  The  other  two  terms  of  Iqq  are  taken  together. 

Here  U  -  const.  Therefore,  changing  the  constant  of  integration  by 
cos~^a  log  kU,  one  considers 

tga  log(k(R-X))  -  cos  ^a  log(k(R-V)) 

-  (tga  -  cos  ^ a)log(k(R-X) )  -  cos  ^a  log( (R-V)/(R-X) 

The  first  term  is  0(Aa)  because  of  Eq.  (F.1).  The  second  term 
gives  RY  ^  +  XY  ^  +  0(Aa)  because  of  Eq.  (F.10). 

Consider  next 

I°Q(a)  +  “  B^U[cos"^a  UY'"*  +  cos'''a  RY'"'  +  tg^a  log(  ( R-V)  [  U  |  ) 

-  sina  cos  ^a  log(k(R-X))j 
cos  ^a  UY  ^  -  -  cos  ^a  +  cos  ^ a  XY  ^  +  (1/2)  +  0(Aa) 
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tg^a  log((R-V)|U|  -  sino  cos  log(k(R-X)) 

-  tg^a  log(CR-V)(R-X)~’ )  -  tg^o  log(klU|) 

+  (tg^a  -  sina  cos  ^a)log(k(R-X) ) 

-  -  cos'^a  (RY~^+XY~^)  -  (1/2)(R^y'^  +  RXY~^) 

-  (1 /2)log(k(R-X))  -  tg^a  log(k|Ul)  +  O(Aci) 

Therefore,  if  one  disregards  constants  of  integrations 

I°Q(ir/2)  -  I°J(ir/2)  I°Q(ir/2)  -  -a^Y[cos~^  aXY~^  +  cos'^o(RY~’  +  XY~b 

-  (1/2)(R^Y“^  +  RXY~^)  -  (1/2)log(k(R-X))  ] 

I^q(h/2)  -  (6^/2)[Y  log(k(R-X))  ♦  XRY'^  +  R^Y~^ ] 

but  R^y“^  -  Y  +  X^Y"^  and  for  a  -  it/2  -  Y  -  const. 

Therefore,  by  another  change  of  the  constant  of  integration, 
one  obtains 

I°Q  -  (a^/2)[Y  log(k(R-X)  +  (XR  +  X^)y“‘')] 

I°^(a)  -  B^U{tga  log((R-V)  |U|~b  -  cos~''a  log(k(R-X))} 

-  a^uKtga  -  cos"’a)log((R-V)  1U~^  I  )  +  cos'^a  log(  (R-V)  (R-X)"’ ) 

-  cos  ^  a  log(k |U 1 ) } 

-  B^UlO(Aa)  -  RY"^  -  XY“^  -  cos~\  log(k|U|)} 

1°  (Tr/2)  -  a^(R  +  X) 


Hence 


The  limiting  process  o+ir/2  is  immediately  obvious  in  I'.., 

212  ~  1 
I  and  ^22’’  other  words,  no  expression  cos  a  or  tga  are 

encountered  in  these  formulae.).  Limiting  processes  are,  however, 

needed  in  Qq.  One  has 

QqCo)  -  (V/U)(log(k(R-X))-1 )  -  tgo  log(k(R-X)) 

+  cos  log((R-V)|Ul  S 

Here 

-tga  log(k(R-X))  +  cos'^a  log( (R-V) | U | )~^ ) 

-  cos  ^a  log( (R-V) (R-X)  -  cos  ^a(log(k|U|)  +  0(Ao) 

-  -RY~^  +  y"^  -  cos“^a  log(k|U|)  +  0(Aa) 

Therefore,  disregarding  constants 

Qg(ir/2)  -  -XY'hog(k(R-X))  -  RY"’ 

-  (1/2) lV^U~^log(k(R-X))-tg^alog(k(R-X)+sinacos*^alog(R-V) |U|“^ ) 

-  (1 /2)V^U~^-tgaVU"^  -  cos"^a  RU"^ 

-  (1/2)  {X^Y"^(log(k(R-X))  -  (1/2))  +  tg^a  log(  (R-V)  (R-X)~S 

-2  2  -2 
+  (sino  cos  a  -  tg  a)log(k(R-V) )  -  sino  cos  a  log(k|Ul) 

-  tgo(-XY“^  -  cosa  R^y“^)  +  cos“\  R(Y~^  +  cosa  XY“^)  +  0(Aa)} 

-  (1/2)  {X^Y~^log(k(R-X))  -  (1/2))  +  tg^a(-co3o(RY~’  +  XY~^ ) 

-  (1 /2)cos^a(R^Y“^  ♦  RXY"^)]  +  ( 1 /2)log(k(R-X) )  +  cos'^a  XY~^ 

+  R^y“^  +  cos'^a  RY“^  +  RXY“^  -  sina  cos'^a  log(k|Ul)  +  0(Aa)} 


-  {1/2)lX^Y  ^log(k(R-X))  -  (1/2)  +  (1/2)log(k(R-X)) 

♦  (1/2)r2y~2  +  (1/2)RXY"2}  -  sina  cos"^a  log(k|U|)  ♦  0(Ao)} 

but  (1/2)R^Y~^  -  (1/2)  ♦  (1/2)X^Y~^. 

Therefore, 

-  (1 /2) {log(k(R-X))(X^Y‘^  ♦  (1/2))  +  (1/2)RXY~^  +  f(o)} 
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With  these  expressions  the  limiting  process  in  1 


0  0 


can  be  carried 


Figure  1.  Half-wing  with 

trapezoidal  elements. 


Figure  2.  Half-wing  with 

triangular  element 


Figure  4.  Subdivision  of  the  wing 
surface  in  which  no 
element  side  is  a  line 
n  =  const. 
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Figure  7.  Two  subdivisions  of  the  wing  surface  into  elements  which 
are  self-similar.  In  Figure  7a  element  boundaries  are 
lines  n  =  const.  In  Figure  7b  element  boundaries  of 
adjacent  elements  lie  on  one  straight  line. 
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trailing  edge. 
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Figure  9.  Subdivision  of  Cn-  and  xy-elements . 


Figure  10.  Numbering  of  the  corners  in  an  ^h-element 
with  one  side  parallel  to  the  n-axls. 
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